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Intermittent Behavior in Oscillators 


Ey W. A. EDSON 


Oscillators of all sorts may, for certain values of the parameters, show low 
frequency disturbances. Usually the disturbance takes the form of a low-fre 
quency interruption of the desired oscillation. By the method here presented 
it is possible to determine whether or not such intermittent behavior will occur 
in a given oscillator and what circuit modifications are required to promote 
stability. The intentional generation of a modulated wave by control of the low 
frequency behavior of an oscillator is also considered. Oscillators of the nega 
tive resistance type are not considered. 


I. INTRODUCTION 


T HAS been known for a long time that all kinds of oscillators are subject 
to the trouble variously referred to as intermittent oscillation, motor 
boating, or squegging. In conventional circuits such as the Hartley the 
phenomenon is most likely to be observed if the grid leak and grid condenser 
are abnormally large. It is found that the time constant of this combina- 
tion must be reduced as the frequency is raised and as the ( of the resonant 
circuit is decreased. At frequencies above a few hundred megacycles the 
problem of producing a practical circuit with suitable margin of stability 
is quite difficult. 

With the advent of the oscillator having automatic output control the 
problem assumed a new aspect.''* By application of an amplitied contro] 
circuit a high degree of constancy of output together with low harmonic 
output is obtained. Satisfactory operation is secured, however, only when 
suitable attention is given to the characteristics of the control circuit. 

The intentional generation of pulses by means of intermittent oscillations 
of relatively high frequency has been studied to some extent, and circuits | 
of this kind are employed in some television systems. Usually the high- 
frequency oscillation is limited to a small portion of the low-frequency cycle, 
the charge stored during this period being allowed to dissipate itself relatively 
slowly during the remainder of the cycle. 

In all of these circuits satisfactory performance depends upon a proper 
proportioning of elements not directly associated with the operating fre- 

1L. B. Argimbau, “An Oscillator Having a Linear Operating Characteristic,” Proc. 
LR.E., Vol. 21, p. 14, Jan. 1933. 


IF ? J. Groszkowski, “Oscillators with Automatic Control of the Threshold of Regenera 
SA tion,” Proc. 1.R.E., Vol. 22, p. 145, Feb. 1934. 
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quency. When continuous oscillation is necessary it is desirable to prov ide 





adequate margin against intermittent operation. When intermittent opera- 
tion is desired the opposite is true. In either case an understanding of the 
same general problem is necessary. 

The present analysis applies only to oscillators of the feedback type. No 
method of extending it to cover negative resistance oscillators such as the 
Dynatron and the Transitron has been found. Relaxation oscillators as 
such are not considered here inasmuch as they are seldom affected by inter- 
mittent operation. No specific frequency limits apply but it is sometimes 
difficult at very high frequencies to achieve desirable values of the constants. 
At very low frequencies oscillators employing automatic output control are 
relatively unsuitable because their performance tends to be unduly sluggish. 

The term linear oscillator is used to indicate an oscillator in which the 
range of operation is controlled within such limits that the harmonic content 
of the output is inappreciable. 

The general equation describing a simple amplitude-modulated wave is 


V = Vo(1 + m sin 2xft) sin 2xFt 


This may be taken as defining the modulation factor m, a complex number 
which is limited to magnitudes between zero and one. 


Il. GENERAL THEORY OF OSCILLATION 

It is found that three separate functions are necessary and sufficient for 
the operation of an oscillator of the feedback type.* These are indicated 
in the block diagram of Fig. 1. 

The amplifier must be provided to overcome the losses of the rest of the 
system. The power output, if any, depends upon the fact that the output 
of an amplifier is greater than the input. 

Selectivity must be provided to insure that the output has a definite 
frequency. Ordinarily a tuned circuit of relatively high ( is used although 
some excellent oscillators employ resistance-capacitance networks. The 
term filter is employed as being sufficiently general to include these extremes. 

A limiter of some form is necessary to establish the level at which sustained 
oscillations occur. In many circuits the functions of amplifier and limiter 
are performed simultaneously in the vacuum tube. In an important class 
of oscillators the limiter is a thermal device such as a tungsten lamp. In the 
Meacham circuit the functions of limiter and filter are combined in a bridge 
employing a tuned circuit and a tungsten lamp. 

To simplify the analysis it is convenient to assume that the amplifier of 
Fig. 1 is completely linear and operates with equal gain at all frequencies 


3’ This topic is discussed more fully in “Hyper and Ultra-High Frequency Engineering.” 
R. I. Sarbacher, and W. A. Edson, John Wiley & Sons, Inc., 1943. 
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from zero to infinity. Similarly the filter is assumed to consist of linear 
circuit elements and to have a detinite curve of loss versus frequency. Asso- 
ciated with this loss characteristic is some specitic phase characteristic.‘ 
The limiter is assumed to have a loss which is independent of frequency but 
which is explicitly related to the input (or output) voltage. 

Although amplifiers having the ideal performance indicated are not physi- 
cally realizable there are no new or unfamiliar concepts involved. Similarly 
the performance of passive networks, such as constitute the filter, has been 
extensively studied and is well understood. It is therefore appropriate to 
devote the following section to the third function. 
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Fig. 1—Functional block diagram of an oscillator. 























Ill. Types or Lrwirers 

The limiters which are now in common use may be separated into four 
relatively distinct groups. 

1. Vacuum tubes in which the gain is decreased by simple overload as the 
level of oscillation rises. This is the most common form of limiter. 

2. Varistors in which the impedance depends upon the instantaneous value 
of current. Copper oxide, thyrite, and electronic diodes are examples. 

3. Thermistors in which the resistance depends upon the rms value of 
current but does not vary appreciably during any one cycle. Carbon and 
tungsten filament lamps are the most common examples. 

4. Vacuum tubes in which the gain is reduced by application of a bias 
which depends upon the level of oscillation. Usually the bias is developed 
by rectifying a portion of the output. 

The limiters of the first two groups depend for their operation upon the 
generation of harmonic voltages and currents. The limiters of the second 


_ 4H. W. Bode, “Relations Between Attenuation and Phase in Feedback Amplifier De- 
sign,” Bell Sys. Tech. Jour., Vol. 19, pp. 421-457, July 1940. 
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two groups operate with very little harmonic distortion. The output of 
oscillators employing such limiters may, therefore, be made quite free fromm 
harmonic voltages. Oscillators of this sort are referred to as linear because 
the tube or tubes serve as simple Class A linear amplifiers. 


IV. CRITERION OF SELF MODULATION 


The block diagram of Fig. 1 is characterized by the fact that the separate 
elements are connected to each other in the form of an endless ring. The 
output may be assumed to come from any of the three junctions. It is this 
fact of closure which complicates the problem of oscillator study. For 
purposes of analysis it is convenient to open the loop as shown in Fig. 2. 
For this example it makes no difference where we choose to make the cut, 
but in actual circuits some caution must be exercised. This matter is dis- 
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Fig. 2— Test for self-modulation in an oscillator. 





cussed more fully later. It is also necessary to choose the impedances of the 
test generator and test detector so that the operation of the components of 
the original system is not disturbed. 

If a continuous wave of suitable voltage and frequency is supplied by the 
test generator it will be found that the terminal voltage of the test detector 
is identical in magnitude and phase with that of the generator. In this 
condition the requirements which are fundamental to oscillators are satis- 
fied. That is, the frequency and level at which oscillation should occur if 
the circuit were closed as in Fig. 1 have been established. The net phase 
shift of the system is zero and the net gain is zero. 


Whether the oscillations so produced would be stable or interrupted is now 
determined by adding amplitude modulation of relatively low frequency and 
very small magnitude to the test generator. It is clear that this modulation 
will be transmitted through the amplifier, filter, and limiter to the test 
detector and that the phase and percentage of the modulation may both be 
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modified. By examining the transmission of a lightly modulated wave for 
various frequencies of modulation it is possible to determine whether or not 
the normal oscillation will be self modulated when the loop is closed as in 
Fig. 1. 

The carrier is held constant at the frequency F and amplitude V for which 
the input and output are identical, and the frequency f of the modulation is 
varied from zero to infinity. In the following treatment it is assumed that 
the significant portion of the characteristic is observed for modulation fre- 
quencies small compared to F. The theory is simplified in this way without 
being seriously restricted in usefulness. The percentage of modulation must 
be held very low so as not to exceed the normal operating range of the limiter. 
The criterion is most conveniently stated in terms of the transmission of the 
modulation envelope which may be considered as a vector quantity. 





Fig. 3—Nyquist diagram showing magnitude and phase of loop transmission 
Legend: U is unstable 
C is conditionally stable 
S is absolutely stable 


A plot of the vector ratio of output to input modulation for various fre- 
quencies is prepared as in Fig. 3. The system characterized by curve (is 
unstable and will generate a self modulated rather than a continuous wave. 
The system characterized by curve S is unconditionally stable and will be 
free from self modulation. The system characterized by curve C is condi- 
tionally stable and may generate either a continuous or an interrupted wave 
depending upon the manner in which the oscillation is started and other 
factors. 

V. ANALOGY OF THE OSCILLATOR TO THE FEEDBACK AMPLIFIER 

The behavior of oscillators of the type here considered is entirely de- 
pendent upon feedback. It is therefore appropriate to review the funda- 
mental principles which apply to feedback in general. 

In the feedback amplitier, negative feedback is applied to improve the 
linearity, stability, impedance, or frequency characteristics. Considerable 
improvements in some or all of the properties may be secured if a consider- 
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able amount of negative feedback is applied and properly controlled. Pos 

tive feedback is sometimes used to increase gain or selectivity, but stabilit 
under such circumstances is poor. Any considerable amount of positive 
feedback results in oscillation. 

The criterion by which stable feedback systems are distinguished fron 
unstable ones has been presented by Nyquist and veritied by others.®:' 
A closed feedback system having input and output terminals is illustrated 
in Fig. 4. In Fig. 5 the loop is opened at some arbitrary point and a tes! 
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Fig. 4+—Typical feedback amplifier. 
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Fig. 5— Test for stability of feedback amplifier. 


oscillator and detector are connected. Here as in Fig. 2 certain precautions 
as to impedance are observed. The test generator must produce a pure 
sinusoidal wave of such small magnitude that no part of the tested system 
overloads and the vector ratio of the detector voltage to the generator voltage 
is observed for a large number of frequencies. The polar plot of Fig. 3 
applies directly to the feedback amplifier except that the radius vector 
represents the transmission of a simple wave rather than of an envelope. 


5H. Nyquist, “‘Regeneration Theory,” Bell Sys. Tech. Jour., Vol. 11, pp. 126-147, 
Jan., 1932. 

®E. Peterson, J. G. Kreer, & L. A. Ware, “Regeneration Theory and Experiment,” 
Proc. I.R.E?, Vol. 22, pp. 1191-1210, Oct., 1934. 





















Pos 


vilit: 


‘itive 
from 


ated 


test 


6 


6 


ions 
pure 
item 
tage 
2.3 
ctor 
ope. 


-147, 


ent,” 











INTERMITTENT BEHAVIOR IN OSCILLATORS 


The conditions of absolute and conditional stability and instability are 
exactly the same as those already given. 

It must be appreciated that Nyquist’s criterion supplies no information as 
to the type or frequency of oscillations which will be generated by an unstable 


system. This is true because the analysis is limited to linear systems. The 
only information imparted is that a very small oscillation of some frequency 
will increase exponentially with time until the amplitude is limited by the 
action of some non-linear device. A small or relatively large shift of fre- 
quency may occur and the oscillation may be regular or intermittent. The 
present work extends the usefulness of Nyquist’s criterion by using it in 
modified form to determine whether or not a particular unstable system 
(oscillator) has or lacks stability as to self-modulation. There is no apparent 
reason why a system lacking in both fundamental and envelope sta- 
bility might not be analyzed a third time for the stability of the 
self-modulation. 


VI. ANALYSIS OF AN OSCILLATOR HAVING AUTOMATIC OUTPUT CONTRO! 


Figure 6 presents a simple form of feedback oscillator having a separate 
rectifier as limiter. For small amplitudes of oscillation the tube operates 
in a linear fashion with cathode self-bias. No bias is produced by the diode 
rectifier until the peak voltage in the coil 13 exceeds that of the bias battery 
B. All voltage in excess of this value is rectified, smoothed by the condenser 
C, and applied to the resistor r as bias. It is seen that a small percentage 
change in the output level may result in a large change in the bias. Accord- 
ingly an output which is quite stable with respect to the tube condition and 
applied voltages, except that of B, is to be expected. 

The stability of this circuit with respect to self modulation is mos. con- 
veniently tested by opening the oscillatory loop at the plate of the tube. 
In so far as the plate resistance of the tube is high with respect to that of the 
associated circuit it is not necessary to control the impedances of the test 
generator and detector extremely accurately. A block diagram equivalent 
to Fig. 6 is presented in Fig. 7. The conditions which must exist for the 
test of stability are shown in Fig. 8. In both those figures it should be noted 
that the gain control is actuated by the input, not the output, of the ampli- 
fier. It is therefore possible for a marked decrease of output voltage to 
result from a small increase of input voltage. This behavior is very different 
from that of the conventional, back-acting, automatic-volume-contro| 
amplifier in which the output change is in the same direction as the input 
change but of reduced magnitude. It is this difference which is the basis 
of most difficulty with amplitude controlled oscillators. 
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Fig. 6— Oscillator having automatic output control. 
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Fig. 7—Block diagram of automatic output control oscillator. 
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Filter 

The filter of Fig. 8 consists of only a single tuned circuit. Its transmission 
is readily represented in terms of the circuit Q by the familiar universal] 
resonance curve. The transmission of a modulated wave through such a 
passive network is conveniently determined by separating the wave into 
its carrier and two sidebands. The carrier will be the frequency F corre- 
sponding to zero phase shift which, in this case, is also the frequency of 
maximum transmission. The sidebands will be shifted in phase by equa! 


1.0 














fy) 10 f¢ 2.0 ) 10 fQ 2.0 
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Fig. 9 Envelope transmission of a modulated wave through a single tuned circuit of 
selectivity Q. 
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Fig. 10—Data of Fig. 9 plotted in polar form. 


and opposite amounts and attenuated according to the frequency f by which 
they differ from the carrier. This behavior is interpreted in Fig. 9 as trans- 
mission and phase shift of the envelope. It is seen that the transmission 
approaches zero and the phase shift approaches 90° as the modulation 
frequency is indefinitely increased. The same data is presented in polar 
form in Fig. 10. Specifically Fig. 10 shows the vector ratio of the modula- 
tion factor m of the output wave to that of the input wave for all frequencies. 
In Fig. 9 the magnitude and phase angle of the ratio are shown separately. 


Limiter 


The limiting action of the tube and diode combination is determined by 
direct circuit analysis. For very low modulating frequencies the condenser 
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C of Fig. Oserves only as a high-frequency by-pass; the direct voltage across ; 
being the instantaneous difference between the peak voltage induced in / 
and that of the stabilizing battery B. For very high modulating frequencie- 
the modulation as well as the carrier is by-passed by C and no modulatio: 
voltage appears across r. Thus the bias is constant and the output wave is 
identical with the input wave. This corresponds to an envelope transmission 


of (1,0). For intermediate values of the modulating frequency the voltage 


developed across r varies in magnitude and phase approximately as if a 
constant current of the modulating frequency f were applied to r and C in 
parallel. 

The output of the amplifier depends not only upon the bias developed 
across r but also upon the input. For systems having a large amount of 
control the action of the bias is predominant. Thus for a low modulating 
frequency the variation of the bias overpowers the initial modulation, the 
phase of the modulation is reversed, and the percentage magnified by the 
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Fig. 11—Envelope transmission of a modulated wave through controlled amplifier 


action of the limiter. In Fig. 11 the envelope transmission is plotted in polar 
form for conditions of relatively large and relatively small amounts of control. 


Loop Transmission 


The separate diagrams of Figs. 10 and 11 are combined in Fig. 12 to de- 
termine the stability of the system. For any chosen frequency f the vector 
of Fig. 10 is multiplied by the vector of Fig. 11 corresponding to the same 
frequency to locate one point of Fig. 12. The resultant vector has an angle 
which is the sum of the two component angles and a magnitude which is the 
product of the two component magnitudes. 

It is seen that the loop may be made to cross the axis considerably to the 
left of the point (1,0) if the points A and A’ of the previous figures cor- 
respond to the same frequency. Similarly the loop may be made to come 
very close to the point (1,0) by increasing the size of C or lowering the 0 
of the tuned circuit so that the points B and B’ correspond to the same 
frequency. With the circuit elements drawn in Fig. 6 the stability margin 
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may be reducec. to zero, but actual looping of the point (1,0) is not indi- 
cated. Parasitic elements, not here considered, can readily affect the 
performance enough to produce instability. 


o LA 


f=0 f=co”’ (1,0) 








Fig. 12—Nyquist diagram applying to Fig. 6. 
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Fig. 13—Hartley circuit. 


VII. ANALYSIS OF THE HARTLEY OSCILLATOR 


The familiar Hartley Oscillator circuit is shown in Fig. 13. In this 
arrangement the tube serves as amplifier and limiter by the action of over 
loading. Harmonic voltages and currents are produced but if the selectivity 
of the tuned circuit is high the voltage returned to the grid of the tube is 
nearly sinusoidal. 

The stability of this circuit is tested in exactly the same way as was that 
of the previous circuit. The loop is opened at the plate of the tube to 
determine the transmission of a modulated signal. If, as is usually the case, 
the coupling of the coil is close, the filter reduces to a single tuned circuit. 
The limiting action results from bias produced by rectification at the grid. 
Accordingly the block diagram of Fig. 7 is directly applicable, and the 
behavior of the filter is correctly given by Fig. 9. 

Generally the circuit operates in class ““C” 


with high bias and large grid 
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voltage swings. If the time constant of the grid-leak-condenser combination 
is long in comparison to the period of a modulation cycle the bias will no‘ 
be able to follow the applied voltage and the modulation of the outpu! 
will be larger than that of the input. Moreover it is in phase with that o/ 
the input. When the modulating frequency is low the bias is able to folloy 
the level of modulation and the output modulation is very small. Thus the 
transmission of a modulated signal is greatest at high modulating frequen 
cies, and the modulation output is in phase with the input. Because of the 





ie) 
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Fig. 14 Envelope transmission of a modulated wave through a grid-leak-biased Class ¢ 
amplifier. 





Fig. 15-—Nyquist diagram applying to Fig. 13. 


action of the grid-leak-condenser network a phase shift at intermediate 
modulating frequencies occurs. This behavior is represented in ‘polar form 
in Fig. 14. 


The stability of the system is determined by combining in Fig. 15 the 


separate diagrams of Figs. 14 and 10. As in the previous system a 
thoroughly stable system results if the element values are such that the points 
A and A’ of Figs. 10 and 14 correspond to the same frequency. If on the 
other hand the elements are such that B and B’ correspond to the same fre- 
quency the curve loops (1,0) indicating instability. In general stability is 
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promoted by increase of the Q of the tuned circuit and by decrease of the 


time constant of the grid-leak-condenser combination. 


VUL. THe Lamp STasitizeD OSCILLATOR 
The circuit of Fig. 16 is of particular interest because the functions of 
amplifier; limiter, and filter are performed separately by units which are 
readily identified with their functions. The present method of analysis 
was developed in connection with this particular circuit. The output 
frequency and amplitude are beth quite stable and the harmonic content 
of the output is low. 


L 
; 


Under operating conditions the gain of the tuned amplifier, which is 
ordinarily in the order of 40 db, is equalled by the loss of the lamp bridge. 
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Fig. 16--Schematic diagram of lamp stabilized oscillator 


The lamps operate at such a temperature that their resistance is slightly 
less than that of the associated linear resistors. If the gain of the amplifiers 
is for any reason somewhat reduced, the current through the lamps decreases, 
the temperature and resistance of the lamps is reduced, and the loss through 
the bridge is reduced to the new value of amplifier gain. 

The d-c characteristic of a lamp bridge is shown in Fig. 17. A curve 
identical with Fig. 17 is observed if the measurement is made with an alter- 
nating current whose period is very short in comparison to the thermal 
time-constant of the filaments. Up to ZL the operation is nearly linear. In 
the region of M the output is essentially independent of the input. At V 
the bridge is nearly balanced and a small percentage change in the input 
voltage results in a large and opposite percentage change in the output. 
It is thus seen that an alternating current having a small superimposed 


modulation of low frequency will result in an output having a considerably 
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When the modulation 


larger percentage modulation in the opposite phase. 
frequency exceeds a few hundred cycles the lamps are unable to follow the 
individual cycles and the output wave is identical in form to the input. {1 


intermediate modulating frequencies the transmission of a modulated wave 
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Fig. 17--D-C characteristics of a lamp bridge. 
0 (1,0) 
Fig. 18 Envelope transmission of a modulated wave through a lamp bridge. 
ie} 0.5 (1,0) 
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Fig. 19—Envelope transmission of a modulated wave through two similar tuned circuits 


of selectivity Q. 


involves a phase shift. The behavior of a typical lamp bridge is presented 
in Fig. 18. 

If the Q of the grid and plate circuits are both relatively high the filter 
circuit may be taken as equivalent to two separate tuned circuits. The 
The combined transmission of the 

Because two tuned circuits are 


transmission of each is given by Fig. 9. 
pair is given in polar form in Fig. 19. 
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employed, the diagram of Fig. 19 differs markedly from that of Fig. 10. 
Specifically the phase shift corresponding to a given value of attenuation is 
vreatly increased. As in previous cases the curve of over-all loop transmis- 
sion may or may not loop the point (1,0) depending upon the relative 
frequency scales. Thus if the points 1 and A’ of Figs. 18 and 19 correspond 
to the same frequency the Nyquist diagram passes near the point (2, 0) 
indicating instability. If the points B and B’ correspond to the same 
frequency the loop passes very near to the point (1,0) and instability is 
likely. 

By making the tuned circuits very selective or by reducing the thermal 
time constant of the lamp circuit the points C and C’ may be made to cor- 
respond to the same frequency. In this case the loop passes to the left of 
the point (1,0) and the system is absolutely stable. The same result may 
be secured more easily by making one of the tuned circuits much more 
selective than the other. This is ordinarily accomplished by increasing the 
(and impedance level of the grid circuit while keeping the and impedance 
level of the plate circuit much lower so as to provide a suitable power output 
to operate the lamp bridge. 


IX. Tur VARISTOR STABILIZED OSCILLATOR 


A circuit which differs from that of Fig. 16 only in that the lamps are 
replaced by varistors is shown in Fig. 20. At low levels of oscillation the 
impedance of the varistors is relatively high, the loss of the limiter is low 
and the amplitude of oscillation rises. At some higher level the varistor 
impedance is reduced, the bridge approaches balance to the fundamental 
frequency, and a stable condition is reached. Because the initial un- 
balance of the bridge is opposite to that of Fig. 16 a reversal of phase is 
necessary to establish oscillation. 

The stable condition reached differs from that of the lamp stabilized 
oscillator in that the varistor goes through its entire range during each high- 
frequency cycle. The lamp resistance changes by only a small amount 
during any one cycle, its resistance depending on an integration of many 
previous cycles. Two important facts arise from this difference. Har- 
monics are produced in the bridge and, in so far as the varistors face react- 
ances of these harmonic frequencies, intermodulation may produce currents 
of fundamental frequency but shifted in phase with respect to the original. 
Thus the bridge may produce a phase shift which is a function of level of the 
oscillation frequency. A degradation of frequency stability results from 
such a condition. More important to the present problem is the fact that 
all modulation frequencies are transmitted alike. A small modulation is 
reversed in phase and magnified by an amount depending upon the bridge 
balance but not upon the modulation frequency. 
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Because the limiter introduces no phase shift it follows that the envelope 
loop transmission is merely an enlarged and reversed copy of that for (he 
filter. This can loop the (1,0) point only if there are at least three shunt 
elements in the filter section. That is, instability can result only if the phase 
shift of the filter system exceeds 180° for frequencies relatively near the 
operating frequency. This circuit is therefore much less likely to produce 
intermittent operation than any other circuit here considered. 
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Fig. 20--Schematic diagram of varistor stabilized oscillator. 


X. NEGATIVE FEEDBACK IN OSCILLATORS 

Because positive feedback is the necessary condition for the operation 
of an oscillator it is not obvious that the application of negative feedback is 
ever desirable. Actually it is frequently possible to introduce negative 
feedback into an oscillator with no loss of performance and under certain 
circumstances advantages are gained. 

The circuit of Fig. 16 serves as a convenient example. Removal of the 
cathode by-pass condenser is likely to reduce the amplifier gain by about 
6 db and to increase the stability of the gain with respect to applied voltages 
by a corresponding amount. Coincident with removal of the by-pass 
condenser the operating level drops a small amount, the bridge loss decreases 
6 db to reestablish equilibrium, and the stabilizing effect of the bridge is cut 
in half. Accordingly the over-all stability of the output with respect to 
applied voltages is unchanged. The advantages gained are that the loss 
which must be held in the bridge is reduced so that stray reactances are less 
likely to disturb the operation, and that the harmonic content of the output 
is reduced. 

Stated in a different way, the output stability of an oscillator using a non- 
feedback amplifier is limited in practice by the bridge balance which may 
be maintained. After this value of gain has been reached additional stability 
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may be secured by supplying increased inherent gain which is offset by 
direct negative feedback. 


XI. DESIGN OF A CONTROLLED OSCILLATOR 


To clarify the material already presented and to convey some additiona! 
concepts an oscillator having a large amount of control will be designed. 
The block diagram is to be that of Fig. 7 and the circuit is to be similar to 
that of Fig. 6. 

It may readily be seen that the gain contro! must satisfy two fundamental 
requirements. It must deliver a d-c bias which increases rapidly with 
increase of the level of oscillation and it must not return any appreciable 
voltage of oscillation frequency. Otherwise the frequency will be affected 
by the elements in the control circuit as well as those in the filter, and the 
performance will be generally poor. Because of its balance a push-pull 
rectifier is helpful in meeting the latter requirement. The principal require- 
ment is achieved by amplification and by the use of a constant counter emf 
or back bias. No bias is produced until the level of oscillation exceeds some 
threshold value. Above this threshold the bias increases approximately 
volt for volt with the peak value of the signal. The same amplifier which is 
used to increase the control may be used advantageously as a buffer so 
that appreciable power outputs may be produced without degrading the 
frequency or amplitude stability. 

It will be assumed that a Q of 100 is available in the coil and that a fre- 
quency of one megacycle is to be generated. The transmission of a modu- 
lated wave in terms of the sideband displacement through such a one-circuit 
filter is shown in Fig. 21. Because the cutoff occurs very slowly it will be 
convenient to incorporate a rapid cutoff in the auxiliary filter of the gain 
control, thus avoiding an excessive phase shift at any one frequency. 

The circuit features already discussed are shown in Fig. 22. A basic 
oscillator with a single tuned coil, a buffer amplifier having little selectivity 
and therefore contributing very little to the equivalent filter section, a source 
of biasing voltage, a balanced rectifier, and an auxiliary low-pass filter are 
shown. The condenser C is only large enough to allow the rectifier to be 
driven without serious loss at one megacycle. It has relatively little effect 
upon the modulation performance. 

It is assumed that the buffer-amplifier, rectifier, etc. are so chosen that a 
modulation of very low frequency of one part per million applied at the 
plate terminal of the oscillator will result in a modulation of one part in a 
thousand returned to that point. This is equivalent to saying that the 
envelope gain is 60 db at low frequencies, and corresponds to 60 db of 
negative feedback in a conventional amplifier. 

The auxiliary filter will be designed to approximate the attenuation and 
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Fig. 23—Characteristics of auxiliary filter. 
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phase characteristics shown in Fig. 23. The choice of this particular shape 
is best explained by reference to Fig. 24 which presents the over-all envelope 
loop transmission of the system. It is seen that the phase shift is relatively 
constant at 90° over a wide band of frequencies and that the gain falls off 
approximately linearly over the same band. In particular the gain becomes 
zero around 5000 cycles whereas the phase does not reach zero below 300,000 
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Fig. 24—-Overall envelope transmission of Fig. 22. 
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Fig. 25—Configuration of auxiliary filter. 


cycles. In terms of Nyquist’s criterion this represents a very stable system 
which is little disturbed by transient effects. A system having even greater 
stability could be achieved by beginning the cut-off at lower frequencies. 
It would then be found that the output was somewhat sluggish in reaching 
a new equilibrium after being disturbed. Such a behavior is not uncommon 
but is generally undesirable. 

Elements which give approximately the characteristics called for in Fig. 
23 are shown in Fig. 25. The peak of loss at one megacycle is contributed 
by the series resonant trap. The rest of the behavior is due to the 0.5 uf 
condenser in combination with the associated resistors. 
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XII. AuXImrary Controt oF THERMALLY LIMITED OSCILLATORS 


In the Meacham and certain other oscillator circuits a thermistor js 
associated with reactive elements in a bridge circuit which functions as both 
limiter and filter. In these circuits a large increase in the frequency stability 
This may sometimes be conveniently expressed as a magnitica- 
tion of the effective Q of the filter. 

The advantages of great frequency stability and good amplitude stability 
of these systems are accompanied by an undesirable tendency towari 
The thermal constants of the thermistor are not 


is observed. 


intermittent operation. 
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suitable envelope stability is likely to impair the frequency or amplitude 
stability for which the circuit is chosen. 

This dilemma may be resolved by the addition of an auxiliary network 
which does not affect the envelope transmission to very low frequencies but 
does modify the behavior at higher frequencies in such a way as to promote 
the stability of the system. 

A simple circuit illustrating the principle appears in Fig. 26. It will be 
noticed that the circuit is so arranged that the average bias applied to the 
tube is only that due to the cathode resistor. The steady voltage developed 
across C, by the rectifier is unable to affect the bias because of the blocking 
Accordingly the rectifier circuit does not affect the norma! 
operating condition, which is characterized by a bridge loss equal to the 


condenser Co. 


The added elements come into play only if there is a tend- 
Then displacement currents of modulation 


amplifier gain. 
ency toward self-modulation. 


frequency flow through C2 in such a magnitude and phase as to modify the 
tube gain and compensate the modulation returned from the bridge. 
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INTERMITTENT BEHAVIOR IN OSCILLATORS 


The exact nature of the control which must be added is best ascertained 
by opening the circuit at the plate of the tube. The loop transmission of a 
modulation envelope may then be determined, either experimentally or 


analytically. If instability is found an auxiliary circuit must be designed 
to produce an over-all system which is stable. In general the elements of 
the auxiliary circuit are to be chosen so that the loop transmission is con- 
siderably less than unity in the region of zero phase. This is ordinarily 
accomplished by increasing the final cutoff frequency at which the over-all 
loop envelope transmission is negligible. 
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Fig. 27—Self-modulating oscillator 


XIII. A Sete Moputatep OSCILLATOR 


The previous sections have been devoted primarily to the problem of 
preventing self-modulation in oscillators. Let us now consider an oscillator 
having envelope instability. The Nyquist diagram indicates that self- 
modulation will occur and tells the approximate frequency of the envelope 
wave. More detailed analysis of the circuit is necessary to determine the 
wave form of the envelope and the manner in which its amplitude is limited. 

If a circuit is to function well as an oscillator the Nyquist diagram for the 
operating frequency must loop the (1,0) point with considerable margin. 
This is necessary so that a small loss of gain will not stop oscillation. At the 
operating level the limiter reduces the loop transmission to unity. In the 
region of (1,0) amplitude stability is favored if the rate of change of gain 
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with respect to level is high. Similarly the frequency stability is favored {j 
the rate of change of phase with respect to frequency is high. 

If a circuit is to function well as a self-modulated oscillator, the above 
conditions must be met and in addition the Nyquist diagram for the envelope 
must meet similar requirements. That is, there must be a limiter and filter 
in addition to the effective amplifier in the envelope system. 

A circuit which meets these requirements is shown in Fig. 27. It is seen 
to be similar to that of Fig. 6 but to have a more complicated low-frequency 
path. The operation is best explained in terms of the relative size of the 
various elements. The by-pass condensers C; and C2 are comparative| 
small. The blocking condensers C3 and C; are quite large. The choke /, 
is large. Thus these elements serve as open or short circuits but do not 
enter into the setting of either of the frequencies. 

The stability tests are carried out by opening the mesh at the plate of the 
tube. At the operating frequency, as defined by the plate coil and condenser 
the loop gain is high at low levels. Thus the fundamental conditions for 
oscillation exist. 

The next step in the analysis is to supply a signal of suitable magnitude 
and frequency to reduce the loop transmission to (1,0). A small modula- 
tion of very low frequency is returned magnified and reversed in phase, as 
with previous systems. The phase of the envelope transmission changes 
with increase of modulating frequency until it is zero at the resonant fre- 
quency of 1, and Cs. At this frequency a considerable gain exists so that 
the Nyquist diagram for the envelope also loops the point (1, 0). 

The tungsten lamp in conjunction with the other impedances of the bridge 
serves to limit the degree of self-modulation. The operating frequency may 
be set by means of Cs in conjunction with a suitable value of Zs. The 
operating amplitude may be controlled by adjustment of the bias battery B. 
The frequency of the self-modulation is set by means of Cs in conjunction 
with 5. 

XIV. Conciusions 

A method of applying known feedback theory to the problem of self- 
modulation in oscillators has been presented. Although the discussion has 
been limited to electrical circuits it is clear that the analysis is applicable 
to other systems, such as electromechanical or mechanical oscillators. 

The analysis has been applied to several familiar oscillators to illustrate 
the method and to clarify some details of their operation. A sample design 
of a bias controlled oscillator is presented to show application to new designs. 

The application of bias control to thermistor stabilized oscillators is 
described. The design of a self-modulated oscillator is undertaken to show 
how intentional modulation may be introduced and controlled. 
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Evaluating the Relative Bending Strength of Crossarms 


i. By RICHARD C. EGGLESTON 
id tilter [3 Ovrr a million crossarms are produced annually in the United States. 
In the open wire lines of the Bell System alone there are now about 
is seen 20 million arms in use. It is natural, therefore, that public utility engineers 
quency should have an interest in the strength of such an important item of outside 
of the plant material; and, consequently, an interest in any tool or means of evalu- 
ratively } ating the strength of such material. It is believed that the moment diagram 
10ke 1, F% js a convenient and reasonably reliable tool for estimating the loads an 
do not F¥ arm will support, for measuring the effect of knots of various sizes and of 





pinhole locations on arm strength, and for answering similar questions 
e of the F relating to the bending strength of crossarms under vertical loads. 


ndenser Ff Two moment diagrams are shown in Fig. 1 for Bell System Type A cross- 
ons lor [# arms; and in the pages that follow are presented the method used in con- 


structing the diagrams and a discussion of their use. While the calculation 
gnitude fF) results apply particularly to the type and quality of arm referred to, they 
nodula- FF would also be of value as a time saving reference in future studies that may 


lase, as FP] be proposed relating to the strength of the same or other types of arms 
hanges involving different knot allowances. 

int tre- The resisting moment of a beam is the product of its section modulus by 
so that the unit stress on the remotest fiber of the beam. The section modulus of a 


beam of uniform cross-section is constant and readily determinable. The 


bridge section modulus, however, of a beam of nonuniform cross-section, such as a 
cy may crossarm, varies because of the different cross-sectional shapes and dimen- 
» The [% sions involved. 

tery B. In this study the following five different shapes were recognized: 

inction (1) Roofed section between pinholes 


(2) Roofed pinhole section 
(3) Roofed brace bolt hole section 
(4) Rectangular pole bolt hole section 
(5) Rectangular section without bolt holes 
The dimensions of the sections investigated were as follows: 


of self- 


ion has 





dlicable 
a Dimensions 
lustrate Section of Arm 

. Minimum Nominal 
» design 
lesigns. f s ; ; (Inches) (Inches 

’ . § @ Roofed section, except at end of arm.... 335 x43 33x43 
itors Is f Roofed section at end of arm..... lat eae ee 35 x4 31 x 49% 
‘o show | 4 Unroofed sections. 33% x 4,5; 31 
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Since there is little, if any, engineering interest in the strength of structura! 
members of maximum size, no investigations were made of sections of 


maximum dimensions. 
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Fig. 1—Moment diagram for Type A southern pine and Douglas fir crossarms per 
Specification AT-7075: 

Graph 1—Resisting moments of arms of nominal dimensions, straight grained and free 
from knots. (Fiber stress 5000 psi) 

Graph 2—Resisting moments of arms of minimum dimensions, having maximum slant 
grain (1” in 8”), and containing knots of the maximum sizes permitted (viz., 
sizes shown at bottom of arm sketch). (Fiber stress 3250 psi) 

Graph 3—Bending moments from a load of 50 pounds at each pin position. 


Section modulus calculations were made of each shape of minimum and 
nominal size, both with and without knots. Tests have shown that, be- 
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RELATIVE BENDING STRENGTH OF CROSSARMS 


cause of the distortion of the grain that occurs around them, knots are fully 
as injurious to the strength of structural timbers as knot holes.' Therefore, 
in dealing with sections containing knots, it was assumed for the purposes 


of this study that the knot extended across the section in the same manner 
as a hole having a diameter equal to the diameter of the knot. It was 
also assumed that the knot was located in, or reasonably close to, the most 
damaging position in the arm section. 

In the calculations of the section modulus of all roofed arm sections, it 
was necessary first to compute the moments of inertia of the whole or parts 
of the top segments of such sections (viz. nominal and minimum sections 
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Fig. 2—Brace bolt hole section containing a }{ inch knot located immediately below 
the top segment (knot and bolt hole shaded). 


between pinholes, and nominal and minimum pinhole sections). Accord- 
ingly, four such computations were made and the results used in calculating 
the section moduli of all the roofed sections investigated. The details of 
the four computations are shown in the Appendix. To insure uniformity 
in the results, the degree of precision used in these computations was con- 
siderably greater than is ordinarily employed in dealing with timber prod- 
ucts. All of the work, however, was done on a computing machine, and it 
was just about as easy to carry the operations to eight decimal places (which 
was the capacity of the machine used) as to a lesser number. As a matter 
of interest in this connection, it was found by actual trial in Computation I 
that absurd results would occur if fewer than five decimal places were used. 
For convenience, all of the section modulus calculations were made in 
tabular form. In such form the procedure employed would not be readily 


' Pg. 6 Dept. Circular 295, U.S. Dept. of Agriculture, “Basic Grading Rules and Work 
ing Stresses for Structural Timbers,” by J. A. Newlin and R. P. A. Johnson. 
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apparent. 
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Therefore, a sample calculation follows showing the metho: o/ 


finding the section modulus of the brace bolt hole section containing 


inch knot. 


Sam ple Calculation 


Referring to Fig. 2, it will be noted that the knot and bolt hole divide the section into 


three parts: the top segment (7 


) and two rectangular portions (R1 and R2). 


The moment 


of inertia (7) of such a compound section about its neutral axis (at a distance c from J/-\/ 

is equal to the sum of the moments of inertia (/7, 7R1 and /R2) of the component parts 
T, R1 and R2 about axes through their own centers of gravity, plus the areas of the con 

ponent parts multiplied by the squares of the distances of their own centers of gravity 


the neutral axis of the compound section. 
of course, by dividing its moment of inertia by the distance (y) from the neutral axis oi 


The section modulus (.8) of this section is found 


section to the most remote fiber. 


Dimensions: 


>> 
a ee ee 


> 


Areas: 


Moments about M — M: 
Tt 
Rirl 


R2r2 


‘ 
Moments of Inertia: 

iT 

IR 

TR2 

T(t — c}* 

Ri(r1 


R2(¢ — 2)? 


-— 


-\2 
—¢)? 


~ 


Section Modulus: 


= 3.1875” (Width of Section) 

= 0.7500” (Diameter of Knot) 

= 3.7625” (See Computation I in Appendix) 
= 0.7000” (d — 2.125” — 0.1875” — k) 

= 1.9375” (2.125” — 0.1875”) 

= 0.1330" (See Computation I) 

= 3.8955” (d+ g) 

= 2.6625” ($ Al + 2.3125”) 

= 0.96875” (3 h2) 


4.09375” (Depth of Section) 


0.7099 sq. ins. (See Computation T) 
2.2313 (bh1) 
6.1758 (bh2 


9.1170 sq. ins. 


2.7054 
5.9408 
5.9828 
14.6890 = 9.1170 c; and hence 
1.6112 


ll 


0.0053 (See Computation T) 


= (0.0911 (b415 + 12) 
= 1.9319 (bh2* + 12) 
= 3.7043 
= 2.4661 
= 2.5490 


= 10.7477 


= 2.48255 (D — ) 
= I = 4.3293 
: 


The same general procedure shown in this sample calculation was fol- 


lowed in dealing with the other cross-sectional shapes. 


For this reason, 


only the final results of the several calculations are presented; although, for 
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RELATIVE BENDING STRENGTH OF CROSSARMS 


TABLE 1.— Section Modulus of Roofed Sections between Pinholes 
ethan of 


‘ _ Knot Diameter—Inches 
hing a} & ‘ 





ie oe fled a ee Doak 2 
Calculation 1: 

ction into (Knots located at 
€ moment top ol section , 
om MM Section Size™: ; de oa. . es 
nent parts Minimum 6.86 5.08 3.57 2.33 1.35 0.04 
Pthe com. End. Min. Pack! 4.78 ; 2.13 a 0.53 
sit's ee a Nominal. baad 5.50 3.91 | 2.59 | 1.54 | 0.76 
iy a 
Segue : Calculation 2: 

4 (Knots located at 

| ~—_ hottom of section) 

cs Section Size*: ; 

£ Minimum 6.11 | 5.24 | 4.42 | 3.71 | 3.05 | 1.92 | 1.05 | 0.47 

B alculation 3: 

a (Knots located im- 

' mediately — below 

+ top segment) 

4 Section Size*: 

i Minimum 8.03 | 5.45 | 4.56 | 3.86 | 3.34 | 2.95 | 2.50 | 2.34 | 2.37 

End. Min. 7.65 3.65 
Nominal 8.60 4.16 


TABLE 2.—Section Modulus of Roofed Pinhole Sections 


Knot Diameter—Inches 


2 So Cella il 


No Knot } j : 1 1} 2 
Calculation 4: 
(Knots vertical) 
Section Size*: 
Minimum 4.50 3.84 a.21 2.63 2.20 
End. Min. 4.29 2.50 
Nominal Py | ae 2.88 
Calculation 5: 
(Knots horizontal) 
Section Size*: 
Minimum 3.63 2.96 2.40 1.97 1.41 1.11 
End. Min. 2.26 1.33 
Nominal 2.76 1.64 


* Section Sizes: 
Minimum = 3%” x 43°)” 
End. Min. = 335” x 4” (viz. minimum at end of arm) 
Nominal = 31” x 4;4” 


convenience, reference is made to the calculations by number in the pages 
that follow. These results are shown in Tables 1, 2, 3 and 4, and a brief 
discussion of the scope and use made of them follows. 


was {ol- 
Ss reason, 
ugh, for 
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TABLE 3.—Section Modulus of Bolt Hole Sections 
Knot Diameter—Inches 


No ; 


Knot . » } ‘ 1 MN 1 
Calculation 6: 
Brace bolt hole 
section 
Section Size*: 
Minimum. 7.97 | 6.47 5.28 4.33 | 3.58 2.62 
Nominal 8.55 4.71 2.78 


Calculation 7: 
Pole bolt hole 
section 23” Knot 3” ky 
Section Size*: 
Minimum. 9.25 7.42 5.63 3.24 1.51 7 
Nominal 9.74 6.05 3.61 1.66 JSS 


* Section Sizes: 


Minimum = 3;%” x 434” 
End. Min. = 37” x 4” (viz. minimum at end of arm) 
Nominal = 3}” x 4,5” 


TABLE 4.— Section Modulus of Rectangular Section without Bolt Holes 
(Calculation &) 


Section Size Knot Diameter Section Modulus 
Minimum (3,;%” x 4,4”) (No Knot) 9.32 
1 8.24 
5 ' By ?- 
3 6.28 
1 5.40 
1} 3.84 
2 2.54 
23 1.51 
3 75 
Nominal (3}” x 43”) (No Knot) 9.78 
1 8.67 
, 7.62 
3 6.64 
1 5.72 
13 4.10 
2 2.74 
23 1.66 
3 .85 


ROOFED SECTIONS BETWEEN PINHOLES 


As indicated in Table 1, three tabular calculations were made for roofed 
sections between pinholes. In Calculations 1,2 and 3 it was assumed that 
the knots present were located (1) at the top, (2) at the bottom, and (3 
immediately below the top segment of the section, respectively. The re- 
sults relating to the 34%” x 4;°)” section are plotted as Curves 1, 2 and 3 
respectively, in Fig. 3. 
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RELATIVE BENDING STRENGTH OF CROSSARMS 


With respect to the knot positions considéred, it is apparent from an exam- 
‘nation of the three curves (Fig. 3) that knots up to approximately 1)” in 
diameter are most damaging when located immediately below the roofed 
portion of the arm; and that the worst position for knots over 13” in diam- 
eter is at the bottom of the arm. However, since under usual loading 
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Fig. 3—Seetions between pinholes. Section modulus of crossarm sections containing knots 
of the sizes shown on the base line and located in the positions indicated. The data 
apply to sections of minimum size (37%%" x 435” 


conditions knots at the bottom of an arm section are in compression, and 


thus would have less influence on strength than they would have on the 
tension side, it was felt that the strength value shown by Curve 2 may be 


ignored; and that the values shown by a smooth curve, combining the values 


2On Page 69 of U. S. Dept. of Agriculture Tech. Bul. 479, “Strength and Related 
Properties of Woods Grown in the United States” by L. J. Markwardt and T.R.C. Wilson, 
is the following statement: “Knots have approximately one-half as much effect on com 
pressive as on tensile strength.” 
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of Curve 3 up to the 13” knot point with those of Curve 1 for 2” and larger 
knots, would be the practical minimum section moduli for roofed sections 
between pinholes. Accordingly, such a smooth curve was constructed and 
is shown as Curve 2 in Fig. 4. The results of Calculations 1 and 3 for nom. 


CURVE | FOR 34% 4% SECTIONS 
CURVE 2 FOR 37% * 4% SECTIONS 
CURVE 3 FOR 33% 4 SECTIONS 


w oa 


SECTION MODULUS -INCHES ? 
a 





1 2 
KNOT DIAMETER - INCHES 
Fig. 4—Sections between pinholes. Section modulus of crossarm sections containing knots 


of the sizes shown on the base line and located in damaging positions. 


inal and arm-end minimum sections were also plotted, and Curves 1 and 3 
drawn for those sections. 
ROOFED PINHOLE SECTIONS 


Two calculations were made for the pinhole sections: Calculation 4, in 
which the knots were assumed to be located adjacent to the pinhole in a 
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bess vertical position; and Calculation 5, in which the knots were assumed to be 
sections - ‘ ; ee pon 
Sad oni immediately below the top segment in a horizontal position. The results 
ted and z - ‘ a ° 
ia of these two calculations are shown in Table 2. It has heretofore been gen- 
or nom- Wee P a a 
7 erally assumed that in pinhole sections knots less than 1” in diameter were 
“| more damaging in a vertical position than in a horizontal position. The 
* results of Calculations 4 and 5, however, show that the horizontal knots 
+ immediately below the top segment are the more damaging. In order to 
compare the effect of knots so located with the effect of knots at the extreme 
“ ” 
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Fig. 5—Pinhole sections. Section modulus of crossarm sections containing knots of the 
sizes shown on the base line and located in damaging positions. 
top of the section, the following two computations assumed 1” and 2” hori- 
zontal knots at the latter location: 
1” Knot at Section Top: 
ing knots ares a 
; .02875 (3.09375) - 
1S = = 1.48156 
6 
land3 } : P 
S = 2.9631 
2” Knot at Section Top: 
.92875 (2.09375)° - 
on 4. in & 1S = = 67857 
és , ; 6 
ole ina = 


1.3571 
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As the section modulus (8) values for sections containing 1” and 2” hori- 
zontal knots located immediately below the top segment are 1.97 and 1.11, 
respectively, (Calculation 5, Table 2) it is apparent that in pinhole sections 
horizontal knots immediately below the top segment are the more dam- 
aging. The results of Calculation 5 were accordingly plotted in Fig. 5 and 
smooth curves drawn to show the section modulus for each of the three 
sections containing knots of any size. 


RooreD Brace Bott HOLE SECTION 


The worst position for knots in the brace bolt hole section was assumed | 
to be substantially the same as in the roofed sections between pinholes; 


and in Calculation 6, the results of which are shown in Table 3, knots up t 
15” in diameter were assumed to be so located, viz. immediately below the 
top segment. 

To check this assumption with respect to worst position, the following 
analysis was made of the minimum sections: 


Distance from top of section: 


To top of bolt hole 1.78” 
To bottom of bolt hole 2.16” 





Distance from bottom of top segment: 


To top of bolt hole 
To bottom of bolt hole 


B3 
1 

It is apparent that any knot ranging in diameter from 1.78” to 2.16”, 
when located at the top of the section, would enter the bolt hole. The 
section modulus of any section containing a knot within that size range 
would be the section modulus of the remaining portion of the section, or 


-) 


', where 6 is the width of the section and d the depth below the bolt hole. 
) 
Thus 


m= =2\2 
S (minimum arm) = = — = 1.9943 
It is also evident that any knot from 1.45” to 1.83” in diameter, when 
located immediately below the top segment, would likewise enter the bolt 
hole; and that the section modulus, on this basis of knot location, would be 
the same for any section containing a knot within the size range mentioned. 
Continuing the analysis the following tests were made: 


2” Knot: 





The distance between the top segment and the bottom of the bolt hole 
of a minimum section is 1.83”. Therefore, a 2” knot located immediately 
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below the top segment would extend beyond the hole; and its effect would 


be the same as in Calculation 3 (Table 1), where the section modulus of a 
section containing a 2” knot similarly located was found to be 2.50. On 
the other hand, since a 2” knot is within the limits 1.78” and 2.16”, the 
section modulus of a section containing such a knot located at its top 
would be 1.99, 


1.78” Knot: 


A knot of this size immediately below the top segment would enter the 
bolt hole since it is within the 1.45” and 1.83” limits, and the section 
modulus value associated with it would be the same as shown in the 
Calculation 6 results (Table 3) for a section containing a 14” knot, or S = 
2.62. But, as evident from previous discussion, the section modulus 
associated with this knot, if located at the top of the section, would be 
1.99, 


1.5” Knot: 


It can be shown that the section modulus of a section containing a knot 
of this size located at the top of the section would be 2.55; and that the 
section modulus associated with a similarly located 1” knot would be 4.55. 
The foregoing analysis for minimum sections may be summarized as 


follows: 


Section Modulus 
Knot Size 


Knot at Top Knot below Top Segment 
(Inches (Inches Inches* 
2.0 1.99 2.50 
1.78 1.99 2.62 
1.5 | 2.55 2.62 
1.0 4.55 3.58 


A study of this summary shows that knots 13” and over are more dam- 


aging when located at the section top; and that knots under 14” are more 
damaging when located immediately below the top segment. The section 
modulus values associated with 23” and 3” knots would be the same as 
shown in the Calculation 1 results (Table 1). 


By a similar analysis for arms of nominal size it can be shown: 

(1) That the more damaging position for knots 13” and under is imme- 
diately below the top segment; 

(2) That the more damaging position for any knot within the diameter 
range from 1.875” to 2.25” and all the larger knots is at the top of the 
section ; 
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(3) That the section modulus associated with 1.875” to 2.25” knots would 
3.25(1.9375)? 
be 


6 
(4) That the section modulus values associated with 23” and 3” knots 


= 2.0334; and 


would be the same as shown in the Calculation 1 results (Table 1). 


CURVE 1 FOR 34'x 4)" SECTIONS 
CURVE 2FOR 3% 49 SECTIONS __ 


oa 


SECTION MODULUS -INCHES? 





' 2 3 
KNOT DIAMETER - INCHES 


Fig. 6—Brace bolt hole sections. Section modulus of crossarm sections containing knots 
of the sizes shown on the base line and located in damaging positions. 


The results of Calculation 6 (Table 3), and of the foregoing analyses, 
together with the Calculation 1 results for 2}” and 3” knots, were plotted in 
Fig. 6 for both minimum and nominal sections. 
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The most damaging position for knots in the pole bolt hole section was 
3” knots assumed to be at the top of the section. They were so figured in Calcula- 


e 1). 
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Fig. 7—Pole bolt hole section. Section modulus of crossarm section containing knots 
of the sizes shown on the base line and located in damaging positions. 
lalyses, 


tted in i i ah ‘ _ ee 
ted tion 7, the results of which are shown in Table 3 and plotted in Fig. 7 for 


both minimum and nominal arms. 
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RECTANGULAR SECTIONS WITHOUT BoLt HOoLEs 


Here too the most damaging position for knots was assumed to be at the | 
top of the section. In Calculation 8 the section moduli of sections contai: 
ing knots from }{” to 3” in diameter were determined for both minimum « 


nominal sections. The results are shown in Table 4. As section modulus 


values for sections containing knots of other sizes than those shown may | 


"> 


; : , : i OE 

found so simply by the formula for rectangular sections, S = 6 7 no curves 
) 

of the results of this calculation were drawn. 


MomeENT DIAGRAMS 


From the results of this study as shown in Table 4 and in Figs. 4, 5, 6 and 
7, section modulus values for clear arms and for arms containing knots of 
various sizes may be read and multiplied by appropriate fiber stresses t 
determine the resisting moments throughout the length of such arms. For 
example, the section moduli of clear arms of nominal dimensions, and oi 
arms of minimum dimensions with the maximum knots permitted under the 
current Bell System crossarm specification (AT-7075) are as follows: 








Arms of nominal size Arms of minimum size wit J 
and free from knots maximum knots % 
Section of Arm 
Section Modulus aie. See emg gales 

Pole bolt hole 9.74 5.63 .” 
Brace bolt holes 8.55 4.33 :” 
Pole pinholes $.1] 3.28 ;” 
Other pinholes in middle section® 5.11 2.38 :” 
End pinholes a0] 1.33 13” 
Other pinholes in end sections® oes 1.41 So 
Unrooted part of middle section 9.78 3.84 1! 
Roofed part of middle section 8.60 2.95 i 
Solid part of brace bolt hole zones! 8.60 4.56 .” 
Between pinholes in end sections 8.60 vkt Zz 
Extreme ends 8.60 2.03 ” 


These section modulus values were used in preparing the moment dia- 
grams shown in Fig. 1. The clear arm of nominal dimensions was also 
assumed to be straight grained. The fiber stress factor used for it was 5000) 
psi, which is the ultimate fiber stress value that has been employed in the 
Bell System for many years for sawn southern pine and Douglas fir. The 


* For the purposes of specifying knot limitations, crossarms under Specification AT-7075 
are divided into a middle section (between brace bolt holes) and end sections (beyond brac« 
bolt holes). 

‘Where a brace bolt hole zone is less than four (4) inches from a pinhole zone, thes 
zones and the portion of the arm between them are considered as a single zone. 
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fiber stress factor used in computing the resisting moments for the arm of 
minimum size with maximum slant grein and maximum knots was 3250 
psi, Which is simply 5000 psi discounted 35°¢ to allow for slant grain of 1” 
in 8” 


count is, of course, unnecessary for the presence of knots, since allowance 


, which is the maximum permitted by Specification AT-7075. A dis- 
for their effect on strength was made in the section modulus values used. 

Since the 5000 psi value is an ultimate fiber stress and not a working 
stress, and since the arms were assumed to be made of clear, straight grained 
material, Graph 1 (Fig. 1) represents an idealized condition. The resisting 
moments shown are probably the maximum that may be expected from 
any commercial lots of southern pine or Douglas tir crossarms,°? notwith- 
standing the fact that the dimensions of some of the arms may exceed the 
nominal specified. With respect to Graph 2 (Fig. 1), the objection may be 
raised that 35°; is not a sufficient discount for a 1” to &” slant of grain and 
that the 3250 psi value makes no allowance for the effect of long continued 
loading. On the other hand, the graph assumes the simultaneous occur- 
rence of the maximum knot in a most damaging position in every section 
of an arm of minimum dimensions and having the maximum slant of grain 
allowed. Since the probability of such simultaneous occurrence of these 
defects and conditions is extremely small, it is felt that the resisting moments 
of Graph 2 represent the minimum strength of any arm of the two species 
concerned that may be furnished under Specification AT-7075. 

Under the assumptions made, Graphs 1 and 2 (Fig. 1) may be regarded 
as the upper and lower limits of the bending strength of specification cross- 
arms. On the same diagram may be plotted the graph or graphs of the 
moments resulting from any given load at each pin position, or any single 
load concentrated at any point on the arm. As an illustration, Graph 3, 
showing the bending moments from a load of 50 pounds per pin, is shown 
in the diagram (Fig. 1). A load of 50 pounds per pin is calculated to be the 
load of size 165 wire coated with ice having a radial thickness of } inch in 
span lengths of 235 feet, or of wire of the same size in 100 foot spans where 
the radial thickness of the ice coating is 5 inch. Since Graph 3 is wholly 
below Graph 2, even an arm of lowest specification quality would support 
the assumed loads with some margin of strength to spare. This margin or 
factor of safety, would, of course, be increased greatly if the quality of the 
arm under consideration approached the quality assumed in Graph 1. As 
previously indicated, the probability is extremely remote that any single 
arm will ever be furnished of a quality as low as assumed in Graph 2. It 

5 Graphs 1 and 2 (Fig. 1) are for southern pine and Douglas fir crossarms. _ It is estimated 


that the resisting moments of comparable graphs for the other woods included in Speciti 
cation AT-7075 should be about 20°; lower. 
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follows, therefore, that the average strength of any lots of southern pine or 
Douglas fir arms produced under Specification AT-7075 may be expected 
to lie well above the Graph 2 limit. 

Graph 2 and a bending moment graph for vertical loads at each pin 
position are of considerable value to the material design engineer, since th 
degree of parallelism between the two will show whether a_ consistent 
strength relationship exists throughout the length of the crossarm. As 4 
matter of interest in this connection, moment diagrams were used as a guide 
in setting the knot limitations shown in Specification AT-7075. 

Resisting and bending moment graphs may also be used to determine the 
location of the critical section of a crossarm by noting the point of coinci- 
dence between a maximum bending moment graph and the resisting moment 
graph fora cleararm. It can be shown by such graphs that this point in al! 
types of Bell System crossarms, designed for vertical loads, is located at the 
pole pinholes. If the comparison were made between a maximum bending 
moment graph and the resisting moment graph of an arm containing all o/ 
the maximum defects permitted, the location of the point of coincidence 
between the graphs might or might not fall at the pole pinholes, depending 
on the magnitude and location of the defects allowed. It should be noted, 
however, that for such arms the critical section locations so determined apply 
only when the arms are actually of the assumed minimum quality; and, 
since the probability of such being the case is so extremely remote, it is 
concluded that the maximum stress or critical section locations in arms of 
that quality are of academic interest only, and that for all practical purposes 
the critical section of any 34” x 44” x 10’ crossarm is located at the pole 
pinhole. 

This conclusion does not mean that every arm broken in service or under 
test will break at the pole pinhole; for, obviously, if some other section is rela- 
tively weaker because of some hidden defect which reduces its section 
modulus or its fiber strength, it will break at such section regardless of any 
mathematical determination of the break location. But the conclusion 
does mean that, generally speaking, when a crossarm breaks the break will 
occur at, or be closely related to, the pole pinholes. To check the accuracy 
of this conclusion, an examination was made of all available crossarm 
strength test data in which the break locations were recorded. The exam- 
ination revealed that, out of 258 arms tested, the breaks in 219, or 85 per 
cent, were either at, or directly related to, the pole pinholes. Six per cent 
of the breaks were located between the two pole pinholes, and 9 per cent at 
points outside the pole pinholes. 

As an illustration of another use to which such a moment diagram may be 
put, the following specific example is cited. Before the present standard 
Bell System specification for crossarms was drafted, it was decided to 
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n pine or ‘ include a new type (“W6") with 16 pin positions. It was felt that, if the 
expected additional pin holes in the type W6 did not unduly weaken the arm, it could 
@ pot only replace the old type “JW” arm with 8 pin positions but also be used 
‘n installations where greater flexibility in wire spacings might be required. 
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Fig. 8—Resisting moments and maximum bending moments for clear JW and W6 
r 85 per @ crossarms. 


per cent 

rcent at @ In order to obtain an estimate of the strength relationship between the two 
| types, strength tests were made of 10 matched arms of each type. The test 

jmay be @ arms were made of air-seasoned, clear Douglas fir. The dimensions of the 

standard crossarm blanks were 3!” x 41” x 20’. In selecting the 10 blanks from 

cided to @ which the test arms were made, only straight grained pieces free from 
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evidence of manufacturing and other defects were chosen. Each blank 
chosen was cut into two 10’ lengths, one of which was made into a JW arm 
and the other into a W6 arm, making 10 matched arms of each type. The > 
tests were made on an Amsler testing machine. The average breaking load | 7 
at the end pinholes was 1159 pounds for the JW arms and 1002 pounds for 
Wo arms. 

At the same time an estimate was made of the theoretical strength relu- 
tionship between the two types by means of the moment diagrams shown in 
Fig. 8. In this igure are shown the graphs of the resisting moments (fiber 
stress factor —5000 psi) of clear JW and clear W6 arms, together with the 
graphs of the bending moments due to the maximum loads these arms would 
withstand when the loads are concentrated at the end pinholes. These 
maximum loads were determined by dividing the moments at the points of 
coincidence between the graphs (critical pole pinhole sections) by the dis- 
tances to the end pinholes. The maximum loads, so determined, are 608 
pounds for the JW arm and 532 pounds for the W6 arm. The fact that 
these loads are low as compared with the actual breaking loads shows, of 
course, that the average ultimate fiber stress developed by these selected 
arms was considerably greater than 5000 psi, which is not surprising in view 





of their exceptionally high quality. However, so far as the information 





sought is concerned—-namely, to determine not the actual strength but the 
strength relationship between the two types—the result would be the same 
regardless of the fiber stress factor used in the moment diagram. 

The ratio of the strength of the W6 arm to that of the JW arm as shown 
both by the actual strength tests and by the moment diagrams was as fol- 


lows: 
Strength Ratio W6 to JW 
Per cent) 

Actual st th test ine x 100 86.5 

‘tual streng ‘sts — - = < 
Actual strength tes 1159 ) 
M li oe x 100 87.5 
I ‘nt diagrams — — = —— 

ome qgiagrams 608 ‘ 


These ratios show a remarkably close agreement between theory and ac- 
tuality and justify the belief that the crossarm moment diagram may be 
employed to obtain reasonably accurate estimates of relative bending 
strength. 


SUMMARY 


The results of this study may be summarized as follows: 
1. The moment diagram is a useful guide in setting specification 
limitations on defects. 
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2. It is shown that the critical section of a crossarm is located at the pole 


pinholes. The practical value of this observation is that it emphasizes the 
need for keeping the pole pinhole sections and the portion of the arm be- 
tween them reasonably free from strength reducing defects. 

3. Only by breaking tests can the actual bending strength of crossarms 
be determined. The relative bending strengths, however, of two or more 
arms of different types or quality may be estimated with sufficient accuracy 
by means of the moment diagram, regardless of the fiber stress used in its 
construction. 

4. If the fiber stress factor employed is dependable, the moment diagram 
may be used to estimate the minimum factor of safety that would obtain 
for an arm of any type or any assumed quality. In this connection, it is 
believed that the strength of Bell System crossarms is well above the mini- 
mum required to support the loads ordinarily carried. 

5. The section modulus curves of Figs. 4, 5, 6 and 7 will simplify the con- 
struction of moment diagrams for arms of the same sizes shown in the figures 
but differing with respect to type and quality. 

The uses listed lead to the general conclusion that the crossarm moment 
diagram is a convenient and reasonably reliable engineering tool. 


APPENDIX 


” 


Computation IT. Moment of Inertia of Top Segment of Minimum (3;'," x 
4335") Section between Pinholes: 
The moment of inertia (/7) of a segment (7) with respect to anaxis 
through its center of gravity and parallel to its base may be found by the 
formula 


IT = I pp = Ax* 


where J gz is the moment of inertia of the segment about the axis BB, A 
the area of the segment and x the distance between the two axes. The 
values Jy", A and x are given by: 


2 sin*® a cos a 

. > a ‘ 

Ios = ‘Ar 1+ : (1) 
a — sin acosa 


A = Wf (2a — sin 2a) (2) 
,r sin’ a 
x= 3 { (3) 
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Fig. 9 
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—Crossarm section between pinholes, 


The significance of r and a in these formulae, and of the other symbols used 
in the computations that follow will be clear from a glance at Fig. 9. 


D = 4.09375” 


b = 3.1875” 

1b = = 1.59375” 
(=a 

r?> = 16.000000 
(1/2 b)? = 2.540039 


p? = 13.459961 
p = 3.068782” 


d=pt+(D-r= 
0.7099 sq. ins. [Area of T by Formula (2) | 
‘ 3.8018” [By Formula (3) | 
eta a: 


~ 


Sin a 


a 

ad 

2a 

Sin? a 

Sin 2a 

Cos a 

Sin a Cos a 
3.7625” 


0.1330” 


Tnx = 10.2054 [By Formula (1)] 


it 


II 


10.2601 


0.0053 


sf = ().39843750 
23° 28’ 49.93” 
0.40981266 radians 
46° 57’ 39.86” 
0.063252925 
0.73089017 
0.91719548 
0.36544507 


(Note: While the results of this and the following computations are shown 
to four decimal places, the actual work was done by machine and carried 
to eight decimal places as mentioned in the text.) 


Since the width of the section in this computation and the radius of its 
3” 


roof is the same as for the minimum 3,3;” x 4” section at the end of the arm, 
the top segments of the two are identical, and the only value that will differ 
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will be the depth (d) of the rectangular portion of the section, which for the 
smaller will be p + (D — r), or 


3.6688 + (4 — 4) = 3.6688” 
Computation IIT. Moment of Inertia of Top Segment of Nominal (3\” x 
4335”) Section between Pinholes: 


As this computation was made in exactly the same manner as Computa- 
tion I, only the results are here shown: 
d = 3.8593” 
g = 0.1317” 
A = 0.7168 sq. ins 
IT 0.0053 


Il 


Computation IIT. Moment of Inertia of Top Segment of Minimum (3 
4335”) Pinhole Section: 


3 . 
i6|6|6St 


It will be noted in Fig. 10 that the top segment is divided into four parts: 
the small segment (7) at the top of the pinhole, the rectangular portion 
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Fig. 10—Crossarm pinhole section. 


Rj, with a width of 5; and a depth of d, and two portions designated Tc. 
The purpose of this computation is to determine the moment of inertia of 
one of the Tc portions with respect to its gravity axis parallel to its base. 
The moment of inertia of the two Tc portions about the axis BB may be 
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found by deducting the moments of inertia of 7, and R, about this axis from 
the moment of inertia of the entire top segment about the same axis. 


1/26, 


D = 4.09375” Sin a = ah 0.16625 
b, = 1.33” a = 9° 34’ 11.49” 

3b; = 0.6605” a = 0.16702554 radians 
r = 4.00” 2a = 19° 8’ 22.98” 

r? = 16.000000 Sin’ a = 0.0045949941 

(1/2 6;)? = 0.442225 Sin 2 a = 0.32787285 

pi = 15.557775 Cos a = 0.98608364 

pi = 3.9443345” Sin a Cos a = 0.16393640 


d (Computation 7) = 3.7625” 
d; = py + (D— 1) — d = 0.2756” 
ry = pi — 1/2d, = 3.8065” 
Area R; = bid; = 0.3665 sq. ins. 
A, = 0.0494 sq. ins. [Area of 7, by Formula (2) | 
x; = 3.9666” [By Formula (3) | 
By Computation J, 77g” = 10.2654 
TT, »»% [Formula (1)] = 0.7777 
Rie ;= sl + Ry? = 5.3126 


6.0903 
21T¢ wn = 4.1751 
The moment of inertia of the 2 7c areas with respect to the axis through 
their own centers of gravity is given by 
21Tc = 2Tcas — 27¢ 2 
where 
2 Tc is the area of the two Tc portions of the top segment and is given by 
in which -1 is the area of the entire top segment as shown in Computation 
I; and 


where, by the principle of moments, 


7 Tx - T1% = Ryn 


rf 2Tc 
in which Tx, T,x; and Ryn; are the moments of the areas of T, T; and Rj, 
respectively, about the axis BB. (7x = Ax of Computation IL.) 
Thus 


2Tc = 0.2940 sq. ins. 
z = 3.7680” 
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As previously shown, 2/T¢g% = 4.1751 
2Tc 2? = 4.1738 


21T¢c = 0.0013 
ITc = 0.0007 
D — r = 0.09375” 
z = 3.7680” 
‘ 3.8618” 
ad = 3.1625” 


g for Tc = 0.0993” 


II 


The results of this computation apply also to the minimum 33,” x 4” pin- 


hole section at the ends of the arm. The depth (d) of the rectangular por- 
tion of the end pinhole sections will be the same as at the extreme ends of the 
arm, viz. 3.0088”. 


’ 3 7 


Computation IV. Moment of Inertia of Top Segment of Nominal (34” x 4;'5”) 
Pinhole Section: 
Since this computation was made in the same manner as Computation 


III, only the results are here shown: 
d = 3.8593” 


g = 0.1019” 
Tc = 0.1630 sq. ins. 
ITc 


II 


0.0008 

















Mathematical Analysis of Random Noise 
BY S. O. RICE 


(Concluded from July 1944 issue) 


PART III 


STATISTICAL PROPERTIES OF RANDOM, NOISE CURRENTS 


3.0 INTRODUCTION 


In this section we use the representations of the noise currents given in 
section 2.8 to derive some statistical properties of /(¢). The first six sec- 
tions are concerned with the probability distribution of 7(¢) and of its zeros 
and maxima. Sections 3.7 and 3.8 are concerned with the statistical prop- 
erties of the envelope of /(¢). Fluctuations of integrals involving /°(1 
are discussed in section 3.9. The probability distribution of a sine wave 
plus a noise current is given in 3.10 and in 3.11 an alternative method of 
deriving the results of Part IIT is mentioned. Prof. Uhlenbeck has pointed 
out that much of the material in this Part is closely connected with the 
theory of Markoff processes. Also S. Chandrasekhar has written a review 
of a class of physical problems which is related, in a general way, to the 
present subject.” 


3.1. Tue Distripution oF THE Noise CuRRENT™ 


In section 1.4 it has been shown that the distribution of a shot effect 
current approaches a normal law as the expected number of events per 
second, v, increases without limit. 

In line with the spirit of this Part, Part III, we shall use the representation 


. 
I(t) = D> (a, cos wat + by sin wy!) (2.8-1) 
n=l 
to show that /(¢) is distributed according toa normal law. This is obtained 
at once when the procedure outlined in section 2.8 is followed. Since a, 
and 6, are distributed normally, so are a, cos w,f and 0, sin w,f when ¢ is 
regarded as fixed. /(t) is thus the sum of 2.V independent normal variates 
and consequently is itself distributed normally. 
2 Stochastic Problems in Physics and Astronomy, Rev. of Mod. Phys., Vol. 15, pp. 
1-89 (1943). 


%3 An interesting discussion of this subject by V. D. Landon and K. A. Norton is given 
in the J.R.E. Proc., 30 (Sept. 1942) pp. 425-429. 
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The average value of J(t) as given by (2.8-1) is zero since d, = b, = 0: 
I(t) =0 (3.1-1) 


The mean square value of /(¢) is 
N ee 


P(t) = ps (a5, COS’ wat + 0; sin’ wt) 


n=l 
v 


= » w(f,)Af (3.1-2) 


— | w(f) df = (0) = Wo 


In writing down (3.1-2) we have made use of the fact that all the a’s and b’s 
are independent and consequently the average of any cross product is zero. 
We have also made use of 


° < 
a, = 6 


n= wfddf, fa = ndf, tn = 2afn 


which were given in 2.8. y(r) is the correlation function of J(¢) and is 
related to w(f) by 


il 


@ 
vy; ¥(r) = [ w(f) cos 2xfr df (2.1-6) 
0 
as is explained in section 2.1. In this part we shall write the argument of 
y(7) as a subscript in order to save space. 

Since we know that /(¢) is normal and since we also know that its average 
is zero and its mean square value is Yo , we may write down its probability 
density function at once. Thus, the probability of /(/) being in the 
range I, J + dl is 

dl 2 
-12/2 
Ga ee (3.1-3) 
V 2m 

‘This 1s the probability ot finding the current between J and J + d/ ata 

time selected at random. Another way of saying the same thing is to state 


that (3.1—3) is the traction of time the current spends in the range /, 7 + d/. 
In many cases it is more convenient to use the representation (2.8-6) 


I(t) => Cn COS (Wal — ¢n), c, = 2w(f,)Af (2.8-6) 
n=l 


in which ¢; , «++ ¢n are independent random phase angles. In order to 
deduce the normal distribution from this representation we first observe 
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that (2.8-6) expresses /(¢) as the sum of a large number of independent ran. 
dom variables 


T(t) xy o- Xe ote eee 4 Xn 


Xn = Gn COS (Wat — Gp) 


and hence that as V — x J/(t) becomes distributed according to a norma 
law. In order to make the limiting process definite we first choose .V an 
Af such that VAf = F where 


:: alec [ w(f) df 


where € is some arbitrarily chosen small positive quantity. We now let 
VN — « and Af — 0 in such a way that VAf remains equal to PF. The: 


ll 
-M- 


2w(f,.)Af cos? (wnt — on) 
(3.1-4) 


II 


Supp fo w(s) ay 


. 
B=|mP +--+ + lave > (Qw( fis)?” | cos (wrt — gn) 2 


< 4a! [bois af 


where the bars denote averages with respect to the ¢’s, ¢ being held constant 
‘ é ; —3/2 : 

If we assume that the integrals are proper, the ratio BA —O0Oas V— =, 

and consequently the central limit theorem* may be used if w(f) = 0 for 


f > F. Since we may make F as large as we please by choosing € sma!! 


enough, we may cover as large a frequency range as we wish. For this 
reason we write © in place of F. 

Now that the central limit theorem has told us that the distribution of 
I(t), as given by (2.8-6), approaches a normal law, there remains only the 
problem of finding the average and the standard deviation: 


T(t) = D0 cn cos (wat — gn) = 0 
i 
N sa 0 Nesaeanaioats . oi 

P(t) = z. Cn COS? (Wnt — Gn) (3.1-5) 
T 


= | wit) df = vo 


* Section 2.10. 
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This gives the probability density (3.1-3). Hence the two representations 


lead to the same result in this case. Evidently, they will continue to lead 

to identical results as long as the central limit theorem may be used. In the 

future use of the representation (2.8-6) we shall merely assume that the 

central limit theorem may be applied to show that a normal distribution 

is approached. We shall omit the work corresponding to equations (3.14). 
The characteristic function for the distribution of /(f) is 


julit) 0.2 2 
ave,é" = exp — . ul (3.1-6) 


3.2. Tue Distripution or J (¢) ann J (¢ + 7) 


We require the two dimensional distribution in which the first variable 
is the noise current /(¢) and the second variable is its value /(f + 7) at some 
later time 7. It turns out that this distribution is normal”, as we might 
expect from the analogy with section 3.1. The second moments of this 
distribution are 


ee er ee | w(f) df 
re (3.2-1) 
we = I(t)T(t + 7) 


= 7, 


The expression for uy is in line with our definition (2.1-+) for the correla- 
tion function: 


T 
fy, ble) thee ot 1()1(t + 1) dt (2.1-4) 
0 


T-2 


In order to get the distribution from the representation (2.8-6) we write 


I = I(t) = XS & cos (wat — gn) 
1 
I, = I(t+1r) = y Cn COS (ant — gp “> WaT) 
1 


* It seems that the first person to obtain this distribution in connection with noise was 


H. Thiede, Elec. Nachr. Tek. 13 (1936), 84-95. 
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From the central limit theorem for two dimensions it follows that J; and /, 
are distributed normally. As in (3.1) 


m=H= Darl wa =v 
1 0 


72 2 
2 = Is = Ti; = Yo (3.1-2) 
N 
we = Kh, = eR Cn AVE, [COS (Wnt — gn) COS (Wnt — gn + wnt)} 
1 
Now the quantity within the parenthesis is 
2 . ‘ . 
COS” (Wal — On) COS WaT — COS (Wal — gn) SIN (Wal — gn) SIN WaT 
and when we take the average with respect to ¢, the second term drops 
out, giving 


wr = >, Cn} COST [ w(f) cos 2afrdf = y,  (3.2-3) 
I 0 


where we have used w, = 2f, and the relation (2.1-6) between w(f) and (+). 
The probability density function for 7; and J; may be stated. From the 
discussion of the normal law in 2.9 it is 
oa: & Ii = vols + | 3.2-4) 
- 2(yo — vs) i 
For a band pass filter whose range extends from f, to fe we have 


fo 
/ Wo Cos 2nfr df 
f 


vr 





sin wT — SIN wet ° 
Wo 5 = (3.2-5) 
™T 


= = sin r(fs — fa) cos wr fo + fa) 
Yo = wife — fa) 
where w» is the constant value of w(f) in the pass band and 
w, = 2nfr (3.2-6) 


We = 2tfa 


According to our formula (3.2-4), 7; and J, are independent when {, 
is zero. For the 7’s which make y, zero, a knowledge of J; does not add to 
our knowledge of J;. For example, suppose we have a narrow filter. Then 


v, = Owhen 7 = [2(fo + fa)] 
¥- is nearly — Yo when 7 = [fe + fol 
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zand J, | @ For the first value of 7, all we know is that /2 is distributed about zero with 

4 7; = Yo. For the second value of 7 /2 is likely to be near —/,. This is 

® in line with the idea that the noise current through a narrow filter behaves 
like a sine wave of frequency }(f, + f.) (and, incidentally, whose amplitude 
fluctuates with an irregular frequency of the order of 3(f, — fa)). The first 

(3.1-2) ® value of + corresponds to a quarter-period of such a wave and the second 

? value toa half-period. By drawing a sine wave and looking at points sepa- 

rated by quarter and half periods, the reader will see how the ideas agree. 
The characteristic function for the distribution of J; and Js is 


ave, Mitr = exp | -* (w+) — da | (3.2-7) 


7 





m. drops ‘The three dimensional distribution in which 





\ i= T(t) 
tg 
(3.2-3) | Ip = I(t + 7) 
ndy(r). : Tz = I(t + 11 + 72) 
rom the where 7; and 72 are given and ¢ is chosen at random is, as we might expect, 
normal in three dimensions. The moments, from which the distribution 
(3.2-4 may be obtained by the method of Section 2.9, are 
| Mu = Me = p33 = Yo 
Miz = Vr, 
Me23 = Vr. 
e wig = O71 + 72) = Yr ire 
(3.2-3) — 
The characteristic function for I; , Is , [3 is 
ave. etl tiselotizs!s 
Yo, 2 P P (3.2-8) 
= exp| -* (21 + 22 + 23) — miet1 22 — wes2223 — mis 2123 
(3.2-6) 3.3. EXPECTED NUMBER OF ZEROS PER SECOND 
We shall use the following result. Let y be given by 
when y, 
t add to y = F(a, ,a2, +++ dy 5%), (3.3-1) 
-. Then 


and let the a’s be random variables. For a given set of a’s, this equation 
gives a curve of y versus x. Since the a’s are random variables we shall call 
this curve a random curve. Let us select a short interval 1, 1 + dx, 














a 
Nm 
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and then draw a batch of a’s. The probability that the curve obtained |), 
putting these a’s in (3.3-1) will have a zero in x, 41 + dx is 


dx [ in| p00, 9; 1) dn (3.3-2 


and the expected number of zeros in the interval (11 , x2) is 


[ ax [ || p(O, n; x) dn (3.3-3) 


Je, 
In these expressions p(£, n; x) is the probability density function for the 
variables 
F(a, t+ dy} x) 
OF (3.3-4) 
Ox 
Since the a’s are random variables so are — and », and their distribution 
will contain x as a parameter. This is indicated by the notation p(&, 7; 4 
These results may be proved in much the same manner as are similar 
results for the distribution of the maxima of a random curve. This method 
of proof suffers from the restriction that the a’s are required to be bounded. 


Ser 
II 


7 = 


Results equivalent to (3.3-2) and (3.3-3) have been obtained independent! 
by M. Kac.” His method of proof has the advantage of not requiring the 
a’s to be bounded. 

Here we shall sketch the derivation of a closely related result: The prob- 
ability that y will pass through zero in x; , +1 + dx with positive slope is 


SP oe Ota er 
: tan 


Ciacielieibicticiiseg a a 
= 
= 





dx | np(O, n; x1) dn (3.3-5) 


We choose dx so small that the portions of all but a negligible fraction 
of the possible random curves lying in the strip (a1, 41 + dx) may be re- 
garded as straight lines. If vy = Eat .1., and passes through zero for x, <x < 


x1 + dx, its intercept on y = Ois x, — ~ where 7 is the slope. Thus é and 


3 [See 


must be of opposite sign and 


Ser 


a <a —- 2 <a + dx 


*S.O. Rice, Amer. Jour. Math. Vol. 61, pp. 409-416 (1939). However, L. A. MacCo! 
has pointed out to me that a set of sufficient conditions for (3.3-5) to hold is: (a) p(&, 7; 
is continuous with respect to (g, 7) throughout the &y-plane; and (b) that the integral 


[ P(an, n; x1) dn 
0 


converges uniformly with respect to a in some interval —a, < a < a2, where a; and 1 
are positive. These conditions are satisfied in all the applications we shall make use oi 
(3.3-5), 

26M. Kac, Bull. Amer. Math. Soc. Vol. 49, pp. 314-320 (1943). 
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According to the statement of our problem, we are interested only in positive 
values of n, and we therefore write our inequality &s 


’ 


—ndx <i <0 


For a given random curve i.e. for a given set of a’s £ and 9 have the values 
given by 


= F(a, eee dy ; %) 


” E | 
. Ox r=2; 


If these values of € and » satisfy our inequality, the curve goes through zero 
in 41, 41 + dx. The probability of this happening is” 


wre 


@ 0 @ 
[ an [ dé p(é, 3 11) = [ {0 — (—n dx)|p(0, n; x1) dy 
0 ndz 0 


where we have made use of the fact that dx is so very small that € is effec- 
... = . ; - 
‘tively zero. The last expression is the same as (3.3-5). 


In the same way it may be shown that the probability of y passing through 
zero in 4, ¥, + dx with a negative slope is 


0 
—dx [ np(O, n; x1) dn (3.3-6) 


Expression (3.3-2) is obtained by adding (3.3-5) and (3.3-6). 
We are now ready to apply our formulas. We let ¢, /(¢) and yg, play the 
roles of x, y, and a, , respectively, and use 
N 
I(t) = D> ex cos (ant — ¢n); c, = 2w(f)af (2.8-6) 
n= 1 
*7 MacColl has remarked that the step from the double integral on the left hand side 


of this equation to the final result (3.3-5) may be made as follows: 
It is easily seen that the probability density we are seeking is 


Fs se j 
d » x)d 
Ee “f. P(E, 0; x ‘| 


Proceeding formally, without regard to conditions validating the analytical operations 
(for such conditions see the footnote on page 52), we have 


d @ 0 am 
—— [ an [ p(é, 3 x) dé = [ np(—nAx, n; x) dn 
dAx Jo giz 0 


and hence the required probability density is 


[ np(0, n; x) dn 
0 
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The first step is to find the probability density function of the two random 
variables 


N 
t= > Cn COS (wnt — on) 
n=l 


(3.3-7) 
7 = I'(t,) = a = Cr Wn sin (wn ti — Yn) 


where the prime denotes differentiation with respect ¢. From section 2.16) 
un = & = yo 


N " — 
—_ i oe ee es 
be = 7° = Cn Wn SiN (Wn ty — Gn) 


n=1 


= D» (2mfx)?w( fr)Af 


+ 4g? I ful) df = —y 


N — r a Sy AE 
— 2 . 
be = tn = -> Cn Wn COS (wn ty = $n) sin (wn fh =a $n) 
n=l 
= 0 

The expression for p22 arises from (2.1-6) by differentiation. In this expres- 
° ” : . ‘ . 
sion Yo denotes the second derivative of ¥(7) with respect to 7 at 7 = 0: 


W"(r) = —4r” [ f’w(f) cos 2nfr df (3.3-8) 
0 
Hence the probability density is 
[—vovo!” esa 
5) = — exp] —- 3.3-9) 
P(e 951) = — ew] at gn] G8 


WP « e . . . 
where Y is negative. It will be observed that the expression on the right 
is independent of ¢. Hence the probability of having a zero in ¢; , ty + di, 


+00 [—Yovo]” sayy dt ae 
uf -e7'? dn = - : 3-10 
+ 7 FL VO itil 
which follows from (3.3-3), is independent of ¢. 


The expected number of zeros per second, which may be obtained from 
(3.3-3) by integrating (3.3-10) over an interval of one second, is 


moe | [tena] 
1f wo] _ afc: 
| YO ‘| ig w(f) af | 


(3.3-11) 
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For an ideal band pass filter whose pass band extends from /, to f, the 
' expected number of zeros per second is 


os) 

+ 
_ 
to 


2pfinaey 
3 fo <— da 


When f, is zero this becomes 1.155 f, and when /, is very nearly equal to 


fy it approaches f, + fa. 


2 . . , : 
In a recent paper M. Kac™ has given a result which, after a slight gene- 


ralization, leads to 
a 1 rm /2 
q Hite Sa as (3.3-13) 
2r Yu 


for the probability that the noise current will pass through the value / 
with positive slope during the interval ¢,¢ + dt. The expected number of 
such passages per second is 
—72/2y ° } ve 
e “°° X |} the expected number of zeros per second] (3.3-14) 
The expression (3.3-13) may also be derived from analogue of (3.3-5) 
obtained by replacing the zero in p(O, 9; x1) by y. 
In some cases the integral 


a 


ee | fw(f) df 


does not converge. 

An example occurs when we apply a broad band noise voltage to a re- 
sistance and condenser in series. The power spectrum of the voltage across 
the condenser is of the form 

oo = 1 5 

w(f) Pre f 
Although yo is infinite, yo is finite and equal to #/2a. A straightforward 
substitution in our formula (3.3-11) gives infinity as the expected number 
of zeros per second. 

Some light is thrown on this breakdown of our formula when we consider 
a noise current consisting of two bands of noise. One band is confined to 
relatively low frequencies, and its power spectrum will be cenoted by 
w(f). The other band is very narrow and is centered at the relatively high 
frequency fo. The complete power spectrum of our noise is then 


w(f) = wif) + A°6(f — fr) 


*°On the Distribution of Values of Trigonometric Sums with Linearly Independent 
Frequencies, Amer. Jour. Math., Vol. LXV, pp 609-615, (1943). 
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where the unit impulse function 6 is used to represent the very narrow ban¢ 
The power spectrum of the narrow band is approximately the same as tha: 
of the wave Av/2 cos 2xfet. 

The integrals occurring in our formula are 


[wong= [anges 


=W+ A 


Lie CAE To See oe 


eee 








[ wary = [ fon a + af 
=U + A’fs 
We suppose that A and f, are such that 
W > 4° 
UK A'pfe. 
Then our formula (3.3-11) gives us the expected number of zeros 


Afe 
2 wie 


We may give a qualitative explanation of this formula if we regard our 
noise current as composed of a small component 


Tz = 2'"A4 cos 2nfot 


due to the narrow band superposed cn a large, slowly varying component 
due to the lower band. Since the r.m.s. value of the second component is 
r1/2 : : : : ° ° . 
W' we may assign it a representative frequency f; and write it approxi- 
mately as 
11/2 
I, = (2W)”™ cos 2nfit 
The zeros of the noise current are clustered around the zeros of the secon( 
wave. Near such a zero 
r\1/2 
I, = &(2W)2nf At 
where Af is the distance from the zero. The oscillations of J; produce zeros 
when | /; | is less than the amplitude of 72 or when 
71/2 | 
A> W''2rfi | At | 
and the interval over which zeros are produced is given by 
1/2 
AW 
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The number of zeros is this multiplied by 2f.. Since there are 2f; such 
intervals per second the number of zeros per second is 


4 1/2 
* AWS, 


This differs from the result given by our formula by a factor of 2/7. 
This discrepancy is due to our representing the two bands by the sine waves 
I, and Ip. 

From this example we obtain the picture that when the integral for Yo 
converges corresponding to .1 — 0, while at the same time the integral for 


” 


y diverges, corresponding to f, > * in such a way that Af.— «, the 


© noise current behaves something like a continuous function which has no 


derivative. It seems that for physical systems the integrals will always 
converge since parasitic effects will have the effect of making w(f) tend to 
zero rapidly enough. The frequency which represents the region where 
this occurs is of the order of the frequency of the microscopic wiggles. 


So far we have been considering the formulas of this section in the most 


' favorable light possible. There are experiments which indicate the possi- 
bility of the formulas breaking down in some cases. Prof. Uhlenbeck has 


pointed out that if a very broad band fluctuation current be forced” to flow 
through a circuit consisting of a condenser, C, in parallel with a series com- 
bination of inductance, L, and resistance, R, equation (3.3-11) says that the 
expected number of zeros per second of the current, /, flowing through R 


1/2 


(and LZ) is independent of R. It is simply —(LC) The differential 
T 


equation for J is the same as that which governs the Brownian motion of a 
mirror suspended in a gas’’, the gas pressure playing the role of R. Curves 
are available for this motion and it is seen that their character depends 
greatly upon the pressure”. Unfortunately, it is difficult to tell from the 
curves whether the expected number of zeros is independent of the pressure. 
The differences between the curves for various pressures indicates that there 
may be some dependence*. 


3.4 THe DistRipuTion oF ZEROS 
The problem of determining the distribution function for the distance 
between two successive zeros seems to be quite difficult and apparently 


** For example, by putting the circuit in series with a diode. 

*® This problem in Brownian motion is discussed by G. E. Uhlenbeck and S. Goudsmit, 
Phys., Rev., 34 (1929), 145-151. 

E. Kappler, Annalen d. Phys., 11 (1931) 233-256. 

* Since this was written M. Kac and H. Hurwitz have studied the problem of the ex 
pected number of zeros using quite a different method of approach which employs the 
“shot-effect” representation (Sec. 3.11). Their results confirm the correctness of (3.3-11) 
when the integrals converge. When the integrals diverge the average number of elec 
trons, per sec. producing the shot effect must be considered. 
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nobody has as yet given a satisfactory solution. Here we shall give some 
results which are related to the general problem and which give an idea oj 
the form of the distribution for the region of small spacings between {hy 
ZeTOS. 

We shall show (in the work starting with equation (3.4-12)) that the 
probability of the noise current, 7, passing through zero in the interva 
r,7 + dr with a negative slope, when it is known that / passes through zer 
at r = 0 with a positive slope, is 


dr o LT? TM > €<~3/2 4 
| ¥ | Fy (Wo — v7) 11 + Heot (—H)| (3.441 


27 L—vYo 
where Wo. and Mo3 are the cofactors of p22 = —o and po3 = —y, in the 
matrix 
~~ e, * % 
lien: iv ae a a 
[\¥; —vtr —vo 0 (3.4-2 
, 


vr —, 0 Yo ’ 
H = Ma{[Mi — M3,)"”. 


We choose 0 < cot '(—H) < a, the value m being taken at 7 = 0, and the 
value r/2 being approached as 7— &. It should be remembered that we 
are writing the arguments of the correlation functions as subscripts, e.¢., 
—y; is really 


—y’"(r) = ax’ | fwlf) cos 2nfr df (3.3-8 


As 7 becomes larger and larger the behavior of J at 7 is influenced less 
and less by the fact that it goes through zero with a positive slope at r = 0). 
Hence (3.4-1) should approach the probability that, for any interval oi 
length dz chosen at random, I will go through zero with a negative slope. 
Because of symmetry, this is half the probability that it will go through 
zero. Thus (3.4-1) should approach, from (3.3-10), 


mmyi/2 
= | | (3.4-3 


as r— *. It actually does this since M approaches a diagonal -matrix 

. fs 2," . 
and both Me; and # approach zero with Mo3/H — Mx2— —Yoo. For a 
low pass filter cutting off at fy (3.43) is 


drfs” (3.44) 


The behavior of (3.4-1) as r > Ois quite a bit more difficult to work out. 
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4 2 2 10 
Woo and Mo3 go to zero as 7, Moy — Mo; as 7 , and consequently H goes 


to infinity as 7 ‘The final result is that (3.4-1) approaches 


(4) 2 
dr* as — vo ] (3.4-5) 
8 —yYovo 


as 70, assuming y exists. Here the superscript (4) indicates the fourth 
derivative at r = 0, 


yi = 168" [ ful f) df (3.4-6) 
0 
For a low pass filter cutting off at f, (3.4-5) is 
dr ‘3 , (2m) (3.4-7) 


When (3.4-1) is applied to a low pass filter, it turns out that instead of + 
the variable 


¢ = Infor, dg = xf, dr (3.4-8) 


is more convenient to handle. Thus, if we write (3.4-1) as pl¢) dy, it fol- 
lows from (3.4-4) and (3.4-7) that 


1 
ple) ~ deV/3 = 0919 as po «x 
(3.4-9) 


p(¢) has been computed and plotted on Fig. 1 as a function of ¢ for the 
range 0 to 9. From the curve and the theory it is evident that beyond 
9 p(y) oscillates about 0.0919 with ever decreasing amplitude. 

We may take p(y) dy to be the probability that 7 goes through zero in 
¢,¢ + dg, when it is known that J goes through zero at ¢ = 0 with a slope 
opposite to that at ¢. p(y) dg exceeds the probability that J goes through 
zero at ¢ = Oand ing, ¢ + dg with no zeros in between. This is because 
p(¢) dg includes all curves of the latter class and in addition those which 
may have an even number of zeros between 0 and g. From this it follows 
that the curve giving the probability density of the intervals between zeros 
must be underneath the curve of p(y). 

A partial check on the curve for p(y) may be obtained by comparing it 
with a probability density function obtained experimentally by M. E. 
Campbell for the intervals between 754 successive zeros. He passed thermal 
noise through a band pass filter, the lower cutoff being around 200 cps and 
the upper cutoff being around 3000 cps. The upper cutoff was rather grad- 
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ual and it is difficult to assign a representative value. The crosses on figuy: 
1 are obtained from his data when we assume that his filter behaves |ike 

low pass filter with a cutoff at f, = 2850, this choice being made in orde; 
to make the maximum of his curve coincide with that of p(y). 

It is seen that some of the crosses lie above p(y). This is probably dy 
to the fact that the actual filter differs somewhat from the assumed low jyas. 
filter. 

On Fig. 1 there is also plotted a function closely related to (3.4-1). |; 
is the low pass filter form of the following: The probability of [passing 


ona (638 





Fig. 1—Distribution of intervals between zeros—low-pass filter 
vag is probability of a zero in Ag when a zero is at origin. 
vag is probability of a zero in Ag when a zero is at origin and slopes at zeros ar 
opposite signs. 
Va = ple), fp = filter cutoff, 7 = time between zeros. 


through zero in 7, 7 + dr when it is known that J passes through zero a! 
r = Ois 


dr Yo po M2 2 2\—3/2 1 I ~- H 3.4-1 
” i t ib ’ = v 
a a y | ¥;) +H tan” Hl] ( 


- 
a 


* ° . T wp 
where the notation is the same as in (3.4-1) and — ; < tan ‘H< 


This curve should always lie above p(y) and the small difference betwee! 
the curves out to g = 4 indicates that 'the true distribution of zeros is give! 
closely by p(y) out to this point. 

When (3.4-1) is applied to a relatively narrow band pass filter or some 
similar device we may make some approximations and obtain an expression 
somewhat simpler than (3.4-1). As a guide we consider our usual idea 
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band pass filter whose range extends from f, to f,. The correlation function 
is given by (3.2-5). 


W a e . 
Y, = —sin mr( fy — fa) cos wr(fo + fa) 
IT (3.2-5) 


Yo = wol fo — fa) 

From physical considerations we know that in a narrow filter most of the 
distances between zeros will be nearly equal to 
yA 1 

Jo ‘i fa 
i.e., nearly equal to the distance between the zeros of a sine wave having 
the mid-band frequency. We therefore expect (3.4-1) to have a peak very 
close to 7,. Wealso expect peaks at 37 , 57, etc. but we shall not consider 
these. We wish to examine the behavior of (3.4-1) near 7, . 

It turns out that Mo; is nearly equal to Mo so that H/ is large and (3.4-1 


becomes approximately 
dr | Yo | ; Moz 
2L—-vo} lvo — vel" 


In order to see that My; is nearly equal to Mo. we use the expressions 


where 7 is near 7). 


Mx = —Wo (Ws — v3) — vo? 
Mos = ¥7 (Ho — 0) te 
Mx + Mos = (Yo — ¥)[(Yo + IY) — Wo) — WE 
(Yo — ¥,)[B + C] 
Mn — Mss = (Yo + ¥)[(Yo — ve)(— vs — Wo) — HF] 
(Wo +y¥,)[- B+ C} 
B= Wb; — vedo 
C = —Wo +o -— ve 


From (3.2-5) it is seen that ¥, may be written as 


¥, = A cos Br, 8 = (fot fa) 


where 87; = mand A is a function of + which varies slowly in comparison 
with cos Br. We see that near 7, ¥; is nearly equal to —Yo. Likewise 
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y. hovers around zero and y, is nearly equal to —W. Differentiating wit! 
respect to 7 gives 

¥, = A’ cos Br — AB sin Br 

v, = (A” — AB’) cos Br — 24’B sin Br 

vo = Ao — AB’, = Yo = Ad 


” . . . 
where -19 and Ap are the values of A and its second derivative at 7 equa! 
to zero. These lead to 


B = (AgA" — AAG) cos Br — 2AgA’B sin Br 
C = (AA” — A") cos* Br — Apdo + (AR — AVS? 


We wish to show that C + Band C — Bare of the same order of magni- 
tude. If we can do this, it follows that Mos — Mo; is much smaller than 
Mo. + Mo; since Yo — Yr, is approximately 2p» while Yo + Yr is quite smal! 
Consequently we will have shown that Mo; is nearly equal to Mx . 

So far we have made no approximations. We now express the slowly 
varying function A as a power series in r. Since Wo and vo. must be zero 
for the type of functions we consider, it follows that 


2 
A= Ay + 5 Ao + - 
A' =1Ay +>> 

2 
A” =A +; Ag? + «+ 


where we neglect all powers higher than the second. Multiplication and 
squaring gives 


2 2 2 ” 
A* — Ao = TAoAo 


"2 


2 
AA” — A” = AyAo +5 (Ase? = Ay°) 


II 


AoAc + F 


2 


AyA” — AAy = 5 (AyAM — At?) = F 


Since, for small r, A and A” are nearly equal to Ao and Aj, respectively 
. . . , . 
we see that the difference on the left is small relative to Ao Aq, i.e., 


|F| << |AoAo | 
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Our expression for B and C become approximately 


B,.= F cos Br — 2Ay Ay Br sin Br 


C = Fos’ Br — AoA sin? Br — Ap Ag Br’ 


When 7 is near 7; , 8t isapproximately r. Hence both C + BandC — B 
s i ° re ° ‘ 
are approximately — A» Ay mw and are of the same order of magnitude. Con- 
sequently Mss and M2; are both nearly equal and 
Moz = WiC + BI 
241 2 
= —Aj Agr 


When this expression for Mo; is used our approximation to (3.4-1) gives 
us the result: If the correlation function is of the form 


¥, = A cos Br 


where A is a slowly varying function of 7, the probability that the distance 
between two successive zeros lies between 7 and 7 + dr is approximately 
dr a 
2 {1 + a@(r — 71)2|3/2 


where a is positive and 


‘ Ay? x 
a= We 2y5 Ti =- 
—ApgT) pb 


For our ideal band pass filter with the pass band f, — fa , 


9 
> (fs + fa 1 
¢e= V3 ema 1S = 
ho =, Ta to + Ja 
‘ ‘ 1 wan 
and the average value of |} 7 — 7! isa =. Thus 
ave.|r — 71 | 1 f—-fa 1 band width 


7 ar} VJ3 (fo + fa) 24/3 mid-frequency 


When the correlation function cannot be put in the form assumed above 
but still behaves like a sinusoidal wave with slowly varying amplitude we 
may use our first approximation to (3.4-1). Thus, the probability that the 
distance between two successive zeros lies between 7 and 7 + dr is approxi- 
mately 

bdr 
re 23/2 
lyo ,< ¥71 
when 7 lies near 7; where 7; is the smallest value of + which makes y, 


approximately equal to —Yo. This probability is supposed to approach 
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zero rapidly as 7 departs from 7; , and 6 is chosen so that the integral over 
the effective region around 7; is unity. 

It seems to be especially difficult to get an expression for the distribution 
of zeros for large spacing. One method, suggested by Prof. Goudsmit, js 
to amend the conditions leading to (3.4-1) by adding conditions that / he 
positive at equally spaced points along the time axis between 0 and _ ;, 
This leads to integrals which are hard to evaluate. For one point between 
0 and 7 the integral is of the form (3.5-7). 

Another method of approach is to use the method of “in and exclusion” 
of zeros between 0 and 7. Consider the class of curves of J having a zero 
at 7 = 0. Then, in theory, our methods will allow us to compute the fun 
tions polt), pil, 7), pol, Ss, 7), associated with this class where 

po(z) dr is sechialisier of curve having zero in dt 

pilr, 7) dr dr is probability of curve having zeros in dz and dr 

por, s, t) dr dr ds is probability of curve having zeros in dr, dr, and ds 
In fact po(r) dr is expression (3.4-10). The method of in and exclusion 
then leads to an expression for Po(7) dr, the probability of having a zero 
at Oand a zero in 7, 7 + dr but none between O and 7. It is 


Po(r) = polr) — if pilr, tr) dr + = | [f po(r, s, r) drds 


-if [ [ ps(r, s, t, 7) dr ds dt + 


Here again we run into difficult integrals. Incidentally, (3.4-11) may be 
checked for events occurring independently at random. Thus if v dz is 


(3.4-11 


the probability of an event happening in d7, then, if v is a constant and the 

° : 2 3 
events are independent, we have po, pi, pe, -*:* given by »,v, v,-:: 
wr 


From (3.4-11) we obtain the known result Po(r) = ve 
We shall now derive (3.4-1). The work is based upon a generalization of 


(3.3-5): If y is a random curve described by (3.3-1), the probability that y 


will pass through zero in x, x4, + dx; with a positive slope and through 
zero in x2 , X2 + dx. with a negative slope is 


+ 0 
—dx, ave [ dm / dn 1 N2 pO, Mm. X15 0, N25 X2) (3.4-12) 
0 2 


where p(&, m1, 1 3 &, ™, X2) is the probability density function for the 


four random variables 


& = F(a), a2, +++, aw; %) 


OF , 
n= ax 7, ‘= 1, 2. 
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fhe x, and x2 play the role of parameters in (3.4-12). This result may be 
established in much the same way as (3.3-5). 

When we identify F with one of our representations, (2.8-1) or (2.8-6), 
of the noise current /(f) it is seen that p is normal in four dimensions. We 
may obtain the second moments directly from this representation, as has 


been done in the equations just below (3.3-7). The same results may be 
obtained from the definition of ¥(+7), and for the sake of variety we choose 
this second method. We set x) = 4, x» = 4; + 7. Then 


g#=8=7O =v 
ho = T(O)I(t+7) =% (3.4-13) 


a1\ (al ae aw am 
™ N2 — (Gr)(S) —_— we a I (t _ r)I (t) dt 


where primes denote differentiation with respect to the arguments. Inte- 
grating by parts: 


T T 
| t+ r)di(t) = (V(t + IO" - I I(t + r)I(t) dt 
0 0 


We assume that / and its derivative remains finite so that the integrated 
portion vanishes, when divided by 7, in the limit. Since 

rv a 

PCO+1) « its) 

T 


we have 
a ” 
nM s , ¥(7) = —y, 
Or" 


Setting 7 = 0 gives 


2 = ” 
m= mm = —Wo 
in agreement with the value of 2» obtained from (3.3-7). In the same 
way 
os a 
fim. = Limit 5 I(t + r)I(t) dt = — Y(r) 
T-2 7 0 Or 
=y, 
; 1 T 
fm = Limit I'(t)I(t + 1) dt 
T-2 T 0 


“ec dia 1 4 , 
= ( ef I'(t + r)I(t) dt 


, 


= —y7, 





66 BELL SYSTEM TECHNICAL JOURNAL 


where we have integrated by parts in getting &,. Setting r = Oand using 

/ . 
Yo = 0 gives 

fim = £m = 0 

In order to obtain the matrix M of the second moments u;, in a form 
fairly symmetrical about its center we choose the 1, 2, 3, 4 order of our 
variables to be £, m, m, &. From equations (3.4-13) etc. it is seen that 
this choice leads to the expression (3.4—2) fer M. 

When we put & and & equal to zero, we obtain for the probability density 
function in (3.4-12) the expression 


My |? 1 , . 
te? exp | -; va (Moni + 2Mosmine + Msi) | 


Because of the symmetry of M, Ma is equal to M33. When, in the integra! 
(3.4-12) we make the change of variable 


Mo Me 2 
~~ soley te “a 


| ag \3/2 pe @ 
dx, dx, | M } | ‘ ax [ dy ve r2~y?2+2(Mo3/Meoo)ry 
rT Mx Jo = yp 


we obtain 


The double integral may be evaluated by (3.5-4). Let 


e = cos” (- ue) =cot'(-H),  H = MalMz — M3)” 
Mx 


where // is the same as that given in (3.4-2). Our expression now becomes 


dx, dx, _|M #? 1 
ta? Mi, — Mi," + Hoot (-H)] 

From a property of determinants 
MxM3 — M3; = |M| (yi — ¥3) 


Using this to eliminate | M | and dividing by 


én fv)" 
dr Yo 


which, from (3.3-10), is the probability of going through zero in x41, «1+ dx 
with positive slope, gives the probability of going through zero in dx with 
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negative slope when it is known that J goes through zero at x; with positive 
slope: 
1/2 
= a [Mis — Mil!?(y8 — y3)701 + A cot (—)) 
2x L—Yo 
This is the same as (3.4-1). 

The expression (3.4-10) is the same as the probability of J going through 
zero in dr when it is known that J goes through zero at the origin with posi- 
tive slope. This second probability may be obtained from (3.4-1) by add- 
ing the probability that J goes through dr with positive slope when it is 
known to go through zero with positive slope. Thus we must add the ex- 
pression containing the integral in which the integration in both m and m» 
runfrom0to «. In terms of x and y this integral is 


a oO 
—r?—y?—2( Mo3/Mo2)r1 
| rdx | dy ye" * os 
0 0 


This is equivalent to a change in the sign of Mo; and hence of 17. After 
this addition we must consider 


1+ Hcot'(—H) +1— H cot H 
2+ H cot’ (—H) — cot™ H] 
2+ lr — 2 cot” H) 

21 + H tan‘ 7] 

and this leads to (3.4-10). 


3.5 MULTIPLE INTEGRALS 


II 


We wish to evaluate integrals of the form 
@ Vs) ae 3 
} r?—2ar,r9—1- — 
J =| ax, [ wes (3.5-1) 
0 0 


Our method of procedure is to first reduce the exponent to the sum of 
squares by a suitable linear change of variable and then change to polar 
coordinates. This method appears to work also for triple integrals of the 
same sort, but when it is applied to a four-fold integral, the last integration 
apparently cannot be put in closed form. 

The reduction of the exponent to the sum of squares is based upon the 
transformation: If* 


hy v1 + hy Do V2 + hs D3, Vs +--+ +h, Daryn 
O + h2Dav2+ +: +thaDazvn (3. 


*1 


wv 
| 
tN 


Xe 


mi =0 +0 oe FO tiem 


* T. Fort, Am. Math. Monthly, 43 (1936), pp. 477-481. See also Scott and Mathews, 
Theory of Determinants, Cambridge (1904), Prob. 63, p. 276. 
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where Dp = 1, Dy = an, Dr» = Dr», and D,, is the cofactor of a., ‘or 
of a,, because they are equal) in D, : 


| ay Gy2°** Gir 
—1/2 
D, = | @y2 de ’ h, = [Dy—1D,| ’ 
| Qi °° Orr | 


then, if none of the D,’s is zero, 
n 
2 2 2 
Do Gre Xe Xe we Dy + ye + +++ + Yn 
1 
From (3.5-2); the Jacobian A(x, +++ n)/A(M1, °** Ya) is equal to D> *. 
Applying our transformation to the exponent: 
x1 = y1 — aDz" "ye 
O + Dz" 5. 
Dz =1- a’ 


ll 


Xe 


Since x2 runs from 0 to « so must v2. The expression for x; shows that y, 
. —1/2 aie : 
runs from a Dy‘“y2 to «. The integral is therefore 


J = D;'” [ dys / ei dy, 
0 an, '!* ye 


We now change to polar coordinates: 


VM p cos 6 


dy, dy. = pdp d@ 


yo = psin 6 


yo > OgivesO SO< 4 


v1 > aDz"* ye gives cot @ > aDz"” 


and obtain 


: cot~!ap;! 2 ) . 
J= p:** | a | pe” dp 
0 


an 1p;'" cot (aDz"”) 


where the arc-cotangent lies between 0 and x. This may be written in the 
simpler form 


— a 
JI = 31 — a) cos’ a = csc ¢ 


oom 
6 


where 


= cos ¢, 


a 
it being understood that 0 < ¢ < x. 
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a ‘or — Other integrals may be obtained by differentiation. Thus from 


Y ) rd) 

Bi —rx?2—y2—2 o ” 

: | ax [ dye * VY ro? = hecsce (3. 
0 0 


we obtain 


ae 

nm 
| 

ws 


ww 
4 
| 

+ 


@ o 
—z2—y2—2ry cos ¢ 2 
I dx | dyxye* Yue? = esc o(1 — gcoty¢ (3.3 
0 
By using the same transformation we may obtain 


I ax | dy al ain = . xa (3.5-5) 


dD, 


Of course, we may expand part of the exponential in a power series and 
jntegrate termwise but this leads to a series which has to be summed in each 
particular case: 


that vy [ dx : dy “" y" et tary 


“Ee retest tt) 


r=() r a 


If we take —1 < R(m) < —}, —1 < Rim) < —4}, the series may be 
summed when a = 1. The result stated just below equation (3.8-9) is ob- 
tained by continuing m and w analytically. 

The same methods will work when the limits are + *. We obtain, when 
m and n are integers, 


+00 +00 
[ dx dy x" y” env? —tey cos ¢ 





0, n + m odd 
r(™ ott) (3.5-6) 
n a 2 
a c \( Vx (sin grt 
In the a ao i 
F{ —n, —m; . - 5 gi : . 3 n + m even 


The hypergeometric function may also be written as 
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By transformations of this we are led to the following expression for ‘he 
integral 7° 
0, n + m odd, 


fmt+t1\.f/n+1 : 
, ? : ? n m 1 
2 tI e(- 


(in gel 3° 9°93 cos: :) » m,n both even, 


ia 


; m 

r(is4)r (1+ ere cata 

ot os -~ cos? yo), 
(sin g)"t™t! wre a ei i 


m,n odd 








As was mentioned earlier, the method used to evaluate the double inte- 
grals may also be applied to similar triple integrals. Here we state two 
results obtained in this way. 


os iy ay [ dz exp[—x° — y — ¥ — exy — Ibex — ays] 


1| 
ae lat+tB+y—-a] 


[ dx [ dy | dz yz exp [—a"” — y’ — 3 — Lexy — 2bsx — 2ayz] 
0 0 0 


= Me[ttertos 2) ~ a 
4 6=—6hh he pit Sa@t+e8+y—7) (3.5-7) 





i 
where 8 and y are obtained by cyclic permutation of a, }, ¢ from % 
cos a sin” Ds a ; 
a= cos - —— bs 
(1 — cy) an b?) 1/2 (1 a cya titan 6?) : 
i a — be ss 
= cot 7 
oy 
where a, 3, y all lie in the range 0, + and where 
tS Bi 
Ds='¢ 1 a =1+2ab—-a -—B— 
b « I 


For reference we state the integrals which arise from the definition of the 
normal distribution given in section (2.9) 


+00 +2 n 1/2 
i dx, ++: » exp | - > Ars Xy x | = | 
. dx, -° ra dX, X¢Xy exp > Ars Xy =| = ! as oe 


(3.5-8)} 
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where the quadratic form is positive definite and a@_ is its determinant. 


{., is the cofactor of a,,. Incidentally, these may be regarded as special 


cases of 


F dx, -°° / aeot(X Oakes *) F (= ix.) 
L- 26 i 20 1 1 


n—171/2 p+ 0 
/ dx | dyy" “f(x + y) 
-x (3.5-9 


) 


(esta 


) 
e 4 bbs | 
F 4x - - Ts r : 
\ |a| 


which is a generalization of a result given by Schlémilch.* 


| 

\ 

> 

(> 
/ 


3.6 DISTRIBUTION OF MAXIMA OF NOISE CURRENT 
Here we shall use a result similar to those used in sections 3.3and 3.4. Let 
y be a random curve given by (3.3-1), 


y = Fla, +++ ay 3%). (3.3-1) 


If suitable conditions are satisfied, the probability that v has a maximum in 

the rectangle (x1, 41 + dx, vi, ¥1 + dy), dx, and dy, being of the same 
‘ . + 32 

order of magnitude, is 


0 
—dx; dy, [ p(n ’ 0, fare dé (3.6-1 ) 


and the expected number of maxima of vy in a < x < 6 is obtained by in- 
tegrating this expression over the range —* < vy < *¥,a oun Sb. 
P(g, 0, ©) is the probability density function for the random variables 


é = F(a, eee, ay : x1) 


OF 
ait 3.6-2 
(7). one 
oF 
Ox? r=2} 


* Héheren Analysis, Braunschweig (1879), Vol. 2, p. 494, equ. (29). 

® 4m. Jour. Math., Vol. 61 (1939) 409-416. A similar problem has been studied by 
E. L. Dodd, The Length of the Cycles Which Result From the Graduation of Chance 
Elements, Ann. Math. Stat., Vol. 10 (1939) 254-264. He gives a number of references 
to the literature dealing with the fluctuations of time series. 
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In our application of this result we replace x and y by ¢ and J as before 
Then 


t =I = > cos (wnt — ¢n) 
1 

g= 1 

t=” 


where the primes denote differentiation with respect to?¢. According to the 
central limit theorem the distribution of £ 7, ¢ approaches a normal law. 
The second moments defining this law may be obtained either from the 
above definitions of £, », ¢, or may be obtained from the correlation/funct ion 
as was done in the work following equation (3.4-13). 


= ” ied 
#2 = yo, v= —o,; iy = 0 


T 
ne = OI"® = Limit ; [ roroa 
To 0 


ne ee es 
Limit pp (T) — I*(0)} = 0 


T 
&¢ = Limit | I(t)I'"(t) dt 
T-2 T 0 


. . Or) " 
Limit —“~ = yy 
ma ae 


Il 


T 
Limit 7 [ I'"(t)1""(t) dt 
T-0 T 


“> | 
to 
Il 


Y 
Limit 5 [ 1 (t)1(t) dt 
T Jo 


T—0 


(4) 
0 


where the superscript (4) represents the fourth derivative. The matrix Vf 
of the moments is thus 


ur 


Yo .. Yo 
M = 10 —yYo 0 


an" 
Vo 0 yw 
The determinant | VM | and the cofactors of interest are 


412 


|M| = —vo(os® — yo”) (3.6-3) 


Mu = —ovs, Mis= Yo", Mas = —vovo 





tes 
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The probability density function in (3.6-1) is 
is before 


p(T, 0, $) = (2n)-*? | |? exp 





(3.6-4) 


1 > 2 
| -; MM (Mul + Mao + 2M Te) | 


and when this is put in (3.6-1) and the integration with respect to ¢ per- 
formed we get 


ng to the | —3/2 ite Soe 
‘mal lay, || @/ dt ie [ M |""e — 
— the ‘is wT (3.6-5) 
function PF M r( ul ) a be ie asf M13 ) 

c ve Ve) + ert Mi Mul 


for the probability of a maximum occurring in the rectangle d/ dt. As is 
'} mentioned just below expression (3.6-1), the expected number of maxima 
1 in the interval ¢, , 42 may be obtained by integrating (3.6-1) from f; to fs 
after replacing x by ¢,and / from — « to + ~ after replacing yvby 7. When 
we use (3.6-4) it is easier to integrate with respect to J first. The expected 
number is then 


‘ wi" [ ? ( =) 
sm tO ee ee < 
ty ” 2 bal ’ a 2 | M! ™ My q 


C (4) V2 
9 bs -_ ‘ 
= (4 — h) an My” =~ * “ | ve »| 








‘a Yo 

| Hence the expected number of maxima per second is 

4 x pe 

i (4) 4/2 | fw(f) df | 

5 1 | Yo _ | 0 

a me ” = ) (3.6-6) 
a dr —wYo 2 

{ | | fen ar| 


atrix M ) Fora band pass filter, the expected number of maxima per second is 


3h — sp" ma 
Bead — 


where f, and f, are the cut-off frequencies. Putting f, = 0 so as to get a 
low pass filter, 


(3.6-8) 
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From (3.6-8) and (3.6-5) we may obtain the probability density function 
for the maxima in the case of a low pass filter. Thus the probability that 
a maximum selected at random from the universe of maxima will lie in 


I, I + dl is 


dl ) 92/8 Sr i y2/2 7 5 1/2 
3¥V 2m = * iy ye l + er! y 8 (3.6-9) 


where 





I = OuTPUT NOISE CURRENT 
W.= RMS VALUE OF I 
I 


OW 


\ 


' 2 3 a 





Fig. 2—Distribution of maxima of noise current. Noise through ideal low-pass filter 


Pr\y) = . , : 
,— dl = probability that a maximum of / selected at random lies between / and J +d/ 
V Yo 


When y is large and positive (3.6-9) is given asymptotically by 
aT V5 rev 
Vino 3° 
If we write (3.6-9) as p,(y) dy, the probability density p,;(y) of vy may be 
plotted as a function of y. This plot is shown in Fig. 2. The distribution 
function PUmsx < vV bo) defined by 


Plex < Vd) = [ply dy 


and which gives the probability that a maximum selected at random is 
less than a specified y~/y) = J, is one of the four curves plotted in Fig. 4. 

If J is large and positive we may obtain an approximation from (3.6-5). 
We observe that 


M (4) 
= - os ie > ¥ 
|M| poo — vo yo 
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function 
lity that P@ so that when J is large and positive 


ll lie in FR eo Muri Ml ee QI? 


@ Also, in these circumstances the 1 + erf is nearly equal to two. Thus re- 
(3.6-9) FF taining only the important terms and using the definitions of the M’s gives 


the approximation to (3.6-5): 


dl dt —Yo % -17/2% 
ral i | Te (3.6-10) 


From this it follows that the expected number of maxima per second lying 
‘ P . 33 P 
above the line 7 = /; is approximately” when /; is large, 


ur q 
= =f! Par *¥0 
2rL Yo (3.6-11) 


—17/2¥9 


=e} X }3|the expected number of zeros of J per second| 


It is interesting to note that the approximation (3.6-11) for the expected 
number of maxima above J, is the same as the exact expression (3.3-14) for 
the expected number of times J will pass through J; with positive slope. 


3.7 RESULTS ON THE ENVELOPE OF THE NOISE CURRENT 
The noise current flowing in the output of a relatively narrow band pass 
filter has the character of a sine wave of, roughly, the midband frequency 
} whose amplitude fluctuates irregularly, the rapidity of fluctuation being 
ass filter. [7§ of the order of the band width. Here we study the fluctuations of the 


Se 
a 





d/+di. | envelope of such a wave. 
First we define the envelope. Let f,, be a representative midband fre- 
quency. Then if 
om = thn (3.7-1) 
the noise current may be represented, see (2.8-6), by 
may be N 
‘ibution | I Cn COS (Wnt — mt — on + wml) 
n=l 


where the components J, and J, are 


dom is | N 
Fig. 4. | Te = D> Cn COS (wnt — wmt — on) 
(3.6-5). § 45 (3.7-3) 


N 
I, = >. en sin (wat — emt — vn) 


n=l 


% This expression agrees with an estimate made by V. D. Landon, Proc. I. R. FE. 29 
(1941), 50-55. He discusses the number of crests exceeding four times the r.m.s. value 
of I. This corresponds to [7 = 16y,. 
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The envelope, R, is a function of ¢ defined by 
R = [Ie + 1" (3.7-4) 


It follows from the central limit theorem and the definitions (3.7-3) of /, 
and J, that these are two normally distributed random variables. They are 
independent since J7,J, = 0. They both have the same standard deviation, 
namely the square root of 


afefe« I wif) df = Yo (3.7-3) 


° wl 


I 


Consequently, the probability that the point (7., 7;) lies within the ele- 
mentary rectangle d/.d/, is 


tell gy EE) ies 
Im Wo “ 


In much of the following work it is convenient to introduce another ran- 
dem variable @ where 


I. R cos 6 
I, = Rsin 6 


(3.7-7) 


Since 7, and J, are random variables so are R and @. The differentials are 
related by 


dIdI, = Rd6dR (3.7-8 


and the distribution function for R and @ is obtainable from (3.7-6) when 
the change of variables is made: 
d0 RdR —R2/2y 
———€.° (3.7-9) 
2r Yo 
Since this may be expressed as a product of terms involving R only and @ 
only, R and @ are independent random variables, @ being uniformly dis- 
, ‘ waa +, 34 
tributed over the range 0 to 27 and R having the probability density 


R —R*/2y 
mi /2¥0 

¥o 
Expression (3.7-10) gives the probability density for the value of the en- 
velope. Like the normal law for the instantaneous value of J, it depends 
only upon the average total power 


(3.7-10) 


w= [wng 


% See V. D. Landon and K. A. Norton, /.R.E. Proc., 30 (1942), 425-429. 





¢ 
3 


Beet Sian 





Ei eke ss a ee 





aa 
ras 


Sa 














aie ac bee nino Ha SH ik 









(3.7-4) 
-3) of /, 


‘hey are 
Viation, 


the ele- 


(3.7-6) 


her ran- 


(3.7-7) 
tials are 
(3.7-8 
5) when 
(3.7-9) 


y and @ 


ily dis- 
3 


3.7-10) 


the en- 
Jepends 








Eo ee ea TO 





a ae 


BS acre Pk oa 5H 


MATHEMATICAL ANALYSIS OF RANDOM NOISE 


We now study the correlation between R at time ¢ and its value at some 
later time ¢ + 7. Let the subscripts 1 and 2 refer to the times ¢ and ¢ + 1, 
respectively. Then from (3.7-3) and the central limit theorem it follows 
that the four random variables 7 , Ja, J, 7.2 have a four dimensional 
normal distribution. This distribution is determined by the second mo- 


ments 
Ti; a Ti = Tes = Ii = Yo = un 
lala = ale = 0 
1 No 
Tale - Ti Ts = ) ps Cn COS (wn T — Wm T) 
& n=) 
“ 3.7-11 
| wf) cos 2x(f — fa)r af = w13 (3.7 ) 
0 
1 N B 
Tale = Tele = 7 bs Cn sin (w, 7 =e: Wm T) 
& n=l 


= | w(f) sin 2x(f — far df = uu 


The moment matrix for the variables in the order 7-1 , Te, Zo, [2 is 


¢ 0 Migs) Mas ] 

0 a 

v= Yo Mis M13 
pe — p14 Yo O 


Mas h413 O yo 


and from this it follows that the cofactors of the determinant | M ! are 


9 


My = Mn = My = My = WoW — wis — pis) ; 
% ‘“ ‘ (3.7-12) 
= po, A = fo — wis — ws 
My. a Mog = 0 
My: = Moy = —p13A 
My = —Mex3 = —uuAd 
| M| = A? 


The probability density of the four random variables is therefore 


1 ae 1 2 2 2 2 
4A exp vA lYo(Ti + 2 + T3 + Ih) 


— us(TiT3 + I2Ts) — Quu(il, — I2]s)) 
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where we have written J; , J>, 7/3, 7s for Ja, Ja, Ta, Ten. 3~We now make 


the transformation 
I; 
Ts = R, sin ™; I, = Rosin 4, 


Ry cos A, T3 = R cos Ao 


and average the resulting probability density over 6; and @ in order to vet 
the probability that R; and R» lie in dR; and dR:. It is 


R, R. dR, dR, [*" - 
: 3 ‘ - ‘ d6, | dO. exp 


1 9 ° : 
abe v4 lo Rj + Wo Rs — 2u13 R,; Re cos (02 — 6;) sie Quis R, Rs sin (05 oe 


’ 
Since the integrand is a periodic function of 4 we may integrate from 
A = 6; to 0. = 6 + 2 instead of from 0 to 2x. This integration gives the 
Bessel function, /o , of the first kind with imaginary argument. The result- 


ing probability density for R; and R, is 


R, R» Ip (® R; 


Re (ER ais + ail”) ep — YR + RD) B73 


where, from (3.7—12), 


2 2 2 
A= Yo —™ Mi3 — M14 


a3 and yyy are given by (3.7-11). Of course, Ri and R, are always positive. 


For an ideal band pass filter with cut-offs at f, and f, we set 


me + fa 


9 ,’ 


~ 


fu w(f) = Wo for fa < D < to 


and obtain 


Yo 


II 


wo( fo — fa) 


tb 4 f 
wo sin w( fo — fa)r 
Mi = / wo cos 2n(f — fn)t df = i fo) 


Sa TT 


tb 
by = | Wo sin 2xr(f = tm) df 0 


a 


The J term in (3.7-13), which furnishes the correlation between R; and R:, 


becomes 
sin x 
Ry R: xv 
a inlay 
Yo sin? 4 
1 ——— 
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where xis t(fy — fa)t. When xisa multiple of x, R, and Ry are independent 
random variables. When x is zero R; and Rs are equal. Hence we may 
say, roughly, that the period of fluctuation of R is the time it takes x to in- 
crease from 0 to or (fs — fa) '. This is related to the result given in the 
next section, namely that the expected number of maxima of the envelope 
is 641 (fp — fa) per second. 


3.8 MAXIMA oF R 


Here we wish to study the distribution of the maxima of R.* Our work 
is based upon the expression, cf. (3.6-1), 


0 
—dR dt [ p(R, 0, R”’)R” dR" (3.8-1) 


for the probability that a maximum of R falls within the elementary rec- 
tangle dR dt. p(R, R’, R’”) is the probability density for the three dimen- 
sional distribution of R, R’, R’’ where the primes denote differentiation with 
respect to ?. 

We shall determine p(R, R’, R’’) from the probability density of J. , £. 
ee Fs I! , which we shall denote by x1, %2,°:: 4%. The interchange 
of I, and I) is suggested by the later work. It is convenient to introduce 
the notation 

@ 
b, = (2x) | wT — fa)" df 
0 (3.8-2) 
bo = Yo 
where f, is the mid-band frequency, i.e., the frequency chosen in the defini- 
tion of the envelope R. 6, is seen to be analogous to the derivatives of 
¥(7r) at 7 = 0. 
From the definitions (3.7—3) of 7. and 7, we obtain the second moments 
= Yo = bo 
= by 


ai =/ 
I 


= 
oa 


N 
42 


1? = > w(fr)Af4e' (fan — fn)” = bs 
1 


Shel 


| se 


= 4 
Gro) wr! orn 


bs 


= * oe Ds 


* Incidentally, most of the analysis of this section was originally developed in a study 
of the stability of repeaters in a loaded telephone transmission line. The envelope, &, 
was associated with the “returned current” produced by reflections from line irregularities. 
However, the study fell short of its object and the only results which seemed worth sal 
vaging at the time were given in reference” cited in Section 3.3. 
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D wl fa) Af2a( fn — fn) = bi 














X1 Xe = ea = 
I 
4X5 = 1,1. = —b, 
N 
mts = ToDo = —Liw(f)afar'(fa — fu)” = —de 
xX = 1,1, = —by 
mx = 1,1, = —bs 
x5x—e = 1-1, = bs 
All of the other second moments are zero. The moment matrix M is thus 
[ bp fb be «(OO 0 0 
by bo —by Q 0 0 
_| —b —b by O 0) 0) 
vals 0 0 0 bo —b —by 
Q 0 0 —b, bo bs 
| 0 0 0 —b, dbs bs | 
The adjoint matrix is 
’ Bo B, — By i) 0 0 1 
B, Boy — B; 0 0 Q 
0 0 0 Bo —B, —B, 
0 0 o -2 BR Bh 
ee 0 o -— A By | 
Bo = — (baby — 53) B By = (bob, — 03)B 
B, = — (bb4 — b2b3)B Bs; = — (bobs — bib.) B 
By = (bibs — b3)B By = = (bob, — 5i)B (3.8-3 
B= bobobs + 2 by bobs 
— b} — bobs — babi 


|M|=B } 
where B is the determinant of the third order matrices in the upper left and 
lower right corners of M. 
As in the earlier work, the distribution of x, - 
dimensions. The exponent is — [2| M | |‘ times 
Bolxi + xi) bh 2By(x1X2 — 1%) — 2Bo(x1x3 + X4X¢6) 
+ Boo(x2 + x5) 


* , % 1S normal in six 


— 2B3(x2x3 — X5x6) (3.8-4) 


+ Bix + 26) 


' 
L 








is thus 


t and 


in. Six 


38-4) 
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In line with the earlier work we set 


x, = 1, = Rceosé xy = I, = Rsiné 
vw = 1, = R’ sind + R cos 66 
xs = I, = R’ cos@ — Rsin 06’ 
x3 = IT, = R” cos 6 — 2R’ sin 06’ 
— Ros 00° — Rsin 60” 
xg = I = R" sin 6 + 2R’ cos 06’ 
— Rsin 00° + R cos 60” 


The angle 6 varies from 0 to 27 and @’ and 6” vary from —* to+x. By 
forming the Jacobian it may be shown that 


dx, dx, +++ dxg = R* dR dR’ dR" dé d6’ do” 
Also, the quantities in (3.84) are 
ata =R 4143 + 44X68 = RR” — RO” 
X1X2 — 4X, = RO’ x +a; = R°? + R’6” 
X2%3 — Xx = RRO’ — 2R”6’ — R'RO” — RO” 
a3 +46 = R’” — 2RR"0” + 4R°0” + 4RRO'O” 
+ Re" + Re” 


The expression for p(R, 0, R’’) is obtained when we set these values of the 
x’s in (3.8-4) and integrate the resulting probability density over the ranges 
of 0, 6’, 0”: 


R*® '¥ 'S im a 
0, R”) = F de” 8-: 
o(R,0,R") = sop | do a8 ; (3.8-5) 


exp — [BoR” + 2B, R°0’ — 2B.(RR” — R’@”) 


+ By R°6” — 2B3R0'(R” — Ré”) 
+ BR” — 2RR"0” + Ro + R’0’”)| 


The integrations with respect to and 6” may be performed at once leaving 
p(R, 0, R’’) expressed as a single integral which, unfortunately, appears to 
be difficult to handle. For this reason we assume that w(f) is symmetrical 
about the mid-band frequency f,. From (3.8-2), 6; and 63 are zero and 
from (3.8-3), B, and Bs; are zero. 
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With this assumption (3.8-5) yields 
+00 
p(R, 0, R”) = R°(2r)*? By” do’ (3.8-6 
1 m i ae 
exp — 5 pz BoR + R(|Bs + 2B.)R0™ — 2B. R") + BR” — RO”)| 


The probability that a maximum occurs in the elementary rectangle ¢R 
dt is, from (3.8-1), p(t, R) dR dt where 


0 
pe eee [ p(R, 0, R")R” dR” (3.8-7 








We put (3.8-6) in this expression and make the following change of variables 


“— By” >g? oe cig By” ” ba 
eee fe * V2B : 
on ee =o (3.8-8 

V 2B, B Vv 2By 
(By + 2B) _ [3 4 ; : ‘ 
=_l = = — z = 3(3 — a) 4 
2B bi L2 283 : 


a = Bo 2B, = bobs 
2B of BB 


where we have used the expressions for the B’s obtained by setting 4; and 
43 to zero in (3.8-3). Thus 


3/2 p@ © 
p(t, R) = ae (#*) y ay | x? dx (3.89) 
02 


exp [—a’s? + 2bex + zy — (x + y)* 





As was to be expected, this expression shows that p(t, R) is independent of /. 
A series for p(t, R) may be obtained by expanding exp 22(y + bx) and 
then integrating termwise. We use 


[ dy a dx x" ye tv® = Va My + OPH + 1) 
0 ~ 0 2 


ae” fetes) 
7 








which may be evaluated by setting 


9 9 « . 3 
x = p cos ¢g, vy=p sin ¢ 





angle dR 


(3.8-7 





variables 
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The double integral in (3.8-9) becomes 


n!b”™ Tim + 3)T(n —m+ 2 


m=0 mM! (n — m)! se 4 


(2s)" 


where 49 = 1 and 


(3); 
> © 


m=0 m! 


1) 


SCs in + 1)b", 


Jess (m— 


O<n  (3.8-10) 


b 


) 


(n + 1)(1 — 8)? — — (1 — 8)*". 


n large 


The term corresponding to m = 0 in (3.8-10) is 7 + 1. 
We thus obtain 


e a*2z? (B:)° 2 2 


(t, R) = 
P 4bobs War 1-0 r( 


(3.8-11) 


We are interested in the expected number, .V, of maxima per second. 
From the similar work for /, it follows that .V is the coefficient of dt when 


(3.8-1) is integrated with respect to Rfrom Oto x. 
dR = ~/2B,b2" dz = (2b) B)'* by*” dz 
[2bo(a° — 1)]'" dz 


Thus from (3.8—7) and 


we find 


= [ p(t, R) dR 
0 


is “ i T 
(a*® — 1° /&\" } (; 
~ (2a)? () > 


a? A 
r(3 


Equations (3.8-11) and (3.8-12) have been derived on the assumption 
that w(f) is symmetrical about /,, , i.e. the band pass filter attenuation is 





84 BELL SYSTEM TECHNICAL JOURNAL 


symmetrical about the mid-band frequency. We now goa step further an 
assume an ideal band pass filter: 


wf)=w fax<f<fo 
w(f) = 0 otherwise 


2fm = fa + fo 


Putting these in (3.8-2) we obtain zero for 6; and 63 and also 


bo = wol(fr — fa) = Yo 


9 
T Wo 


by = —— (fr — fa)” 


nr Wo 


= (fr — fay 


9 


3(3 ‘scan a’) = 
[2bo(a” — 1)]"?z = [$yo]"”?s 


leh 


6.775 

8.333 
3: 9.9002 
38 7 ~=11.4736 
~ 1.5811 2 + .3953 


From (3.8-12) we find that the expected number of maxima per secon 
of the envelope is 


64110 (fy — fa) (3.8-15 

assuming an ideal band pass filter. . 
The distribution of the maxima of R for an ideal band pass filter may be | 
obtained by placing the results of (3.8-14) in (3.8-11). This gives 


dR —f. T vont 
p(t, R) dR = -. # er 











rther and 


(3.8-14 


(3.8-15 


r may be 
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1 It is convenient to define y as the ratio 


R R 


~ rms. 1) yi? 


= (§)'"2 


where R is understood to correspond to a maximum of the envelope. Since 
the value of R corresponding to a maximum of the envelope selected at 
random is a random variable, y is also a random variable. Its probability 
density is pr(y), where 


p(t, R) dR 


Pal») dy = 0 64410(f — fa) 


pre(y) has been computed and is plotted as a function of vin Fig. 3. 


ENVELOPE OF OUTPUT 


NOISE CURRENT 


Pe L — 
20 y 2.5 3.0 35 40 
Fig. 3—Distribution of maxima of envelope of noise current. Noise through ideal band- 
pass filter. 
Pr) 
/ 


R + dR. 


dR = probability that a maximum of R selected at random lies between R and 


The distribution function P(Rnax < wv yo) defined by 


= y 
P(Resox s yo) _e [ pry) dy 


and which gives the probability that a maximum of the envelope selected 
at random is less than a specified value vy = R, is plotted in Fig. 4 to- 
gether with other curves of the same nature. 
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[he asymptotic expression for px(v) may be obtained from the integral 
3.8-9) for p(t, R). Indeed, replacing the variables of integration x, y in 
3.8-9) by 

x 


y=auty, 


t B 


'@ integrating a portion of the y’ integral by parts, and assuming 6 < 1 
£ g s \ 


(a > 1, by Schwarz’s inequality, so that 6 < 1 always) leads to 


b, 5 e R2 Wo R’ 
um ~ (BP (E-1) 
P 2a Yo Yo 
when R is large. 


If, instead of an ideal band pass filter, we assume that w{f) is given by 


1 ( 2/92 
ttt fn > 6 (3.8-16) 


wf) _ ov/ Ie e ; 


we find that 
1 
= dao 
16r°- 30° 
3,6=0 
(2 + 1) 
Some rough work indicates that the sum of the series in (3.8-12) is near 
3.97. This gives the expected number of maxima of the envelope as 
N = 2.52¢ (3.8-17) 
per second. 


The pass band is determined by ¢. It appears difficult to compare this 
with an ideal band pass filter. If we use the fact that the filter given by 


‘iia = Wo CX le owt = P 
wf) vexp| 7(i-*) | 


passes the same average amount of power as does an ideal band pass filter 
whose pass band is fs — fa, we have 


fo “ee Je = ov/ 2 
and the expression for V becomes 1.006 (fs — fa). 


3.9 ENERGY FLUCTUATION 


Some information regarding the statistical behavior of the random vari- 
able 
t1+T 


E= I°(t) dt (3.9-1) 
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where /(¢) is a noise current and ¢; is chosen at random, has been given in 4 
recent article.” Here we study this behavior from a somewhat differen: 
point of view. 
If we agree to use the representations (2.8-1) or (2.8-6) we may write, as 
in the paper, the random variable E as 
T/2 
E= [ P(t) dt (3.9-2) 
7/2 
where the randomness on the right is due either to the a,’s and 6,’s if (2.8-1 
is used or to the ¢,’s if (2.8-6) is used. 
The average value of E is mr where, from (3.1—2), 


_ T/2 ae 7/2 
E =m, = / [2(t) dt = [ ¥(0) dt = Tyo 
7/2 T/2 


wo (3.9-3) 
=f [ w(f) df 
The second moment of £ is 
= T/2 T/2 
B= [dof dm P@Pe (3.04) 
T/2 T/2 


If, for the time being, we set ¢. equal to ¢,; + 7, it is seen from section 3.2 
that we have an expression for the probability density of 7(¢,) and I(t; + + 
and hence we may obtain the required average: 


+00 +20 
ii = aq | an, [ dl 1713 exp 


1 . 3.9-5 
(- — (voli + WoT) - hI) ( 
2A? 
A=p—-ve, L=T(h), I= Mh + 7) = M(t) 
The integral may be evaluated by (3.5-6) when we set 
2 2 
I, = Ax —, I,=A = 
1 x 2. y 3. 
(3.9-6) 


A =yYsing 


* “Filtered Thermal Noise—Fluctuation of Energy as a Function of Interval Length”, 
Jour. Acous. Soc. Am., 14 (1943), 216-227. 
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Thus 
Lil; = (1 + 2 cos® ¢) 
=o + 2, 
Incidentally, this gives an expression for the correlation function of J°(t). 
Replacing 7 by its value of tf, — t; and returning to (3.9-4), 


(3.9-7) 


r/2 T/2 
B= Ty + 2 | at, [ dtzy"(t, — t) (3.9-8) 
7/2 7/2 


When we introduce o7 , the standard deviation of E, and use 
2 = 2 
or =F — my, 


we obtain 


: : : T/2 T/2 ’ 
oy = (E — EY 2| dt, dt (t2 — ty) 
T/2 


T/2 


T 
i/ (T — x)¥*(«) dx 


where the second line may be obtained from the first either by changing the 
variables of integration, as in (3.9-27), or by the method used below in 
dealing with £3. I am indebted to Prof. Kac for pointing out the advantage 
obtained by reducing the double integral to a single integral. It should be 
noted that the limits of integration — 7/2, T/2 in the double integral may 
be replaced by 0, T by making the change of variable t = t’ — T/2 for both 
ty and to. 
When we use 


¥(r) = I w(f) cos 2xfr df (2.1-6) 


we obtain the result stated in the paper, namely, 


oe oes . we ; sin’ (fi + fo)T a 
=| w(fr) a: | w( fs) afe| rt ere (3.9-9) 


sin® mf — fr 
m(fi — fr)? 

If this formula is applied to a relatively narrow band-pass filter and if 

T(fo — fa) >> 1 the contribution of the f, + fe term may be neglected and 


we have the approximation 


oe . Po ii . sin” wf, — fe)T 
em [mah [weds Ga 


(3.9-10) 


wo T (fs — fa) 


Wo Mr ' 
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where, from (3.9-3) 
mr = wWol (fo — fa) (3.9-11) 
The third moment £3 may be computed in the same way. However, in 
this case it pays to introduce the characteristic function for the distribution 
of I(t), [(t2), T(tz). Since this distribution is normal its characteristic 
function is 


Average exp [iz) I; + ite IT, + i23 Js] 


Yo /.2 2 2 a 
= exp -|% (21 + 22 + 23) + V(tz — t)21 2% (3.9-12) 


+ ¥(ts — A)aszs +V(b — ha | 


From the definition of the characteristic function it follows that 


9 


LLL = —coeff. of r: 


omeee 


in ch. f. 


vi + WoW + vii + vb) elie 


+ 8Ya1 Wai 32 


where we have written Ya for ¥(t2 — 4), etc. When (3.9-13) is multiplied 
by dt; dt, dts , the variables integrated from 0 to T, and the above double 


integral expression for a; used, we find 
. T T ¥ 
(E — ky = 212 | dt, I dts [dist bards ° 
0 0 0 
Denoting the triple integral on the right by J and differentiating, 


T id 
OF ns 3 | dt, | dboWty — tV(T — t)W(T — tr) 
d7 0 0 


= 3[ ax | dy)(x — y)¥(x)v(y) 


6[ ax [ dy (x — y)v(x)y(y) 


In going from the first line to the second ¢, and ¢ were replaced by T — x and 
T — y, respectively. In going from the second to the third use was made oi 


the relations symbolized by 


a T T Ir T y 
| dx dy = | dx dy + dx / dy 
0 0 0 0 0 x 
T r T y 
= | av [ dy + | dy [ dx 
0 0 0 0 
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and of the fact that the integrand is symmetrical in x and y. Integrating 
dJ dT with respect to T from 0 to 7), using the formula 


T; T T; 
[ ar [ f(x) dx = [ (T; — x)f(x) dx, 
0 0 0 


noting that J is zero when T is zero, and dropping the subscript on 7; finally 


wives 
: T r 
(E — E)’ = as | ax [ dy(T — x)¥(xW(v)y(x — vy). 
0 0 


E4 may be treated in a similar way. It is found that 


T T T T 
(- BY - 38 — B= 30 [ dt [de [dts [ drdadnveve 
0 0 0 0 


which may be reduced to the sum of two triple integrals. It is interesting 
to note that the expression on the left is the fourth semi-invariant of the 
random variable E and gives us a measure of the peakedness of the dis- 
tribution (kurtosis). Likewise, the second and third moments about the 
mean are the second and third semi-invariants of E. This suggests that 
possibly the higher semi-invariants may also be expressed as similar multiple 
integrals. 

So far, in this section, we have been speaking of the statistical constants 
of E. The determination of an exact expression for the probability density 
of E, in which T occurs as a parameter, seems to be quite difficult. 

When T is very small E is approximately /°(t)T. The probability that 
E lies in dE is the probability that the current lies in —7, —7 —d/ plus the 
probability that the current lies in J, 7 + d/: 

2dI rP 
V 2m ll 2Po 


where E is positive, 


O\ 1/2 a 
I = (5) ; dimes (ET)? dE 


 /. 


’ E 
= (2m ET)” exp — ol dE (3.9-14) 


and T is assumed to be so small that /(¢) does not change appreciably during 
an interval of length T. 

When T is very large we may divide it into a number of intervals, say », 
each of length T/n. Let E, be the contribution of the r th interval. The 
energy E for the entire interval is then 


E = Ey + Fo + on + Ey 
If the sub-intervals are large enough the £,’s are substantially independent 
random variables. If in addition m is large enough E is distributed nor- 
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mally, approximately. Hence when T is very large the probability that £ 
lies in dE is 


dE } (E — mr)* 3 0-15 
nf weit 
where 
Mr = r | w(f) df 
0 
(3.9-16) 


7=T I : w'(f) df 


Q 
— 
ll 


the second relation being obtained by letting T — © in (3.9-9),. The 
analogy with Campbell’s theorem, section 1.2, is evident. When we dea! 
with a band pass filter we may use (3.9-10) and (3.9-11). 

Consider a relatively narrow band pass filter such that we may find a 7 
for which Tf, >> 24 but T(fs — fa) << .64. Thus several cycles of fre- 
quency f, are contained in T but, from (3.8-15), the envelope does not change 
appreciably during this interval. Thus throughout this interval J(¢) may 
be considered to be a sine wave of amplitude R. The corresponding value 
of E is approximately 

R’ 
Ra T= 
2 
where the distribution of the envelope R is given by (3.7-10). From this 
it follows that the probability of EZ lying in dE is 
dE E dE . . 
—_*«—s (3.9-17) 
Yor mr 
when £ is small but not too small. 
When we look at (3.9-14) and (3.9-17) we observe that they are of the 


form 


a”*'k" 


—ak . > 

sete @ OE 3.9-18) 

T'(n + 1) ( 

Moreover, the normal law (3.9-15), may be obtained from this by letting » 

become large. This suggests that an approximate expression for the dis- 

tribution of E is given by (3.9-18) wher. a and m are selected so as to give 
the values of my and o7 obtained from (3.9-3) and (3.9-9). This gives 





ic n+1=™ (3.9-19) 
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and if we drop the subscript T and substitute the value of a in (3.9-18) we 


get 


( y" 
oe mE mE me 

mnt si age 3.9-20 
I'(n + 1) i ( o ) . ( o ) ‘ . o va ’ 


An idea of how this distribution behaves may be obtained from the 
following table: 


X95 X75 

n T(fo — fa) X25 X50 X35 — a 

0 0 .29 .695 1.39 415 2.00 

1 1.45 .96 1.68 2.69 oa 1.60 

2 2.4 1.73 2.67 3.94 647 1.47 

3 3.4 2.54 3.67 S42 .692 1.39 

5 5.4 4.22 5.67 7.42 .744 1.31 

10 10.5 8.63 10.67 13.02 . 808 1.22 
24 25 21.47 24.67 28.17 .870 1.14 
48 50 44.1 48.7 53.5 .905 1.10 


where is the exponent in (3.9-20). The column 7( /, — f,) holds only for a 
narrow band pass filter and was obtained by reading the curve vy, in Fig. 1 
of the above mentioned paper. The figures in this column are not very 
accurate. The next three columns give the points which divide the dis- 
tribution into four intervals of equal probability: 


mE 95 . ee ee eee 
Xo = , E25 = energy exceeded 75°% of time 


2 
id 


mE 5 simtes ° 
Xs = —. , Eso = energy exceeded 50° of time 
o 
mE 75 © 4s 
X73 =-~, Ez = energy exceeded 25% of time 
a 
The values in these columns were obtained from Pearson’s table of the in- 
complete gamma function. The last two columns show how the distribu- 
tion clusters around the average value as the normal law is approached. 
For the larger values of we expected the normal law (3.9-15) to be 
approached. Since, for this law the 25, 50, and 75 per cent points are at 
m — .6750, m, and m + .6750 we have to a first approximation 


‘ 
— 5 = (n +1) = T(fo — fa) 


m s -- y— 3.9-2 
ta = 3? (m — 6750) = X.59 — 675+ x29 ( 1) 


X75 = X90 + .675>/ x. 


This agrees with the table. 
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Thiede” has studied the mean square value of the fluctuations of the 
integral 
t 
At) = [ Poet? dr (3.9-22 
p— I 
The reading of a hot wire ammeter through which a current J is passing js 


proportional to A(¢). @is a constant of the meter. Here we study A(¢) by 










675 





-|—-|-y. —- — P —- oc 
os ae Ye T (f= fa) 
0.403 ATO 
PROBABILITY DENSITY 
(3.9-20) 1s ASSUMED 
: kk = 0 20 30 ao. 5¢ 
TE f.) 


Fig. 5*—Filtered thermal noise—spread of energy fluctuation 


t,+T 
E= | I(t) dt, t; random, T is noise current. 


"1 
Se E.75/E.s0 , 7 E..25/E.s0 . 


ty — fa = band width of filter. 


first obtaining its correlation function. This method of approach enables 


us to extend Thiede’s results 
The distributed portion of the power spectrum of A(f) is given by (3.9 
30). When the power spectrum w(f) of J(¢) is zero except over the band 
fa < f < fy where it is wo , the power spectrum of A(t) is 
Wwf fe-f) , 
—___—..—— for 0<f<f- 
a + 4r°f? f<h—Sfe 
and is zero from f, — fa up to 2f,. The spectrum from 2f, to 2f, is not zero, 
and may be obtained from (3.9-34). The mean square fluctuation of -1\/ 
is given, in the general case, by (3.9-28) and (3.9-32). For the band pass 


case, when (fp — fa)/a@ is large, 


r.m.s A(t) — A o | = ¥ 
<a ad A 2(fo — fa) 


88 Elec. Nachr. Tek., 13 (1936), 84-95. This is an excellent article. 
* Note added in proof. The value of y2 at 0 should be .415 instead of .403. 
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We start by setting r = ¢ — uw which transforms the integral for A(t) into 
aD 
Alé) = [ I(t — ue *" du (3.9-23) 
0 


In order to obtain the correlation function ¥(7) for 1(¢) we multiply -1(¢ 
by A(t + 7) and average over all the possible currents 


W(r) = AACE + 7) 


Il 


@ @o 
3 [ ” a au | e' dvave. F(t — WI (t + 7 — 2) 
0 0 


Just as in (3.9-4) the average in the integrand is the correlation function of 
4 T(t), the argument being¢ + 77 -—27—-t+u=7r+u—-%. From (3.9-7 
it is seen that this is 


Yo + Wr + u — 2) 
where (7) is the correlation function of 7(t). Hence 


W(r) = Yo 2 | du | dve™ "V(r t+u—rz) (3.9-24) 
a 0 0 


4 From the integral (3.9-23) for A(t) it is seen that the average value of 
A(t) is 








em ( . 
es (3.9-25) 
Qa Q@ 
where we have used 
enables 3G ve , ’ 
q w=vo=/[ whag=P 
0 
yy (3.9 
e band Using this result again, only this time applying it to A(¢), gives 
A*(t) = ¥(0) 
oe « @ ‘ (3.9-26) 
=A + 2[ au | dve  “V'i(u — 2) 
0 0 
ot zero, ' 
of Ai The double integrals may be transformed by means of the change of 
nd pass variable w+ 7 = x,u —v = y. Then (3.9-24) becomes 


V(r) = A” +{[° dy a dx + [ dy [ ix eu W(r + y) 


a (3.9-27) 
=A 44 [ coy ty) + ¥% - lay 
a 40 
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When we make use of the fact that ¥(y) is an even function of y we see, from 
(3.9-26), that the mean square fluctuation of A(¢) is 


Te ae 7 a 
(4) —AP = a) - A? = | eM y) dy (3.9-28 
Qa 0 


(7) may be expressed in terms of integrals involving the power spectrum 
wif) of Z(t). The work starts with (3.9-24) and is much the same as in 
going from (3.9-8) to (3.9-9). The result is 

@ @ 
W(r) = A? + | df, I dfyw( fr) w( fr) 
| cos 2r(fi + fo)r + cos 2r( fi — fe)r 
a? + [2r(fit fe) a + [24(f — fr) 


It is convenient to detine w(—f) for negative frequencies to be equal to 
w(f). The integration with respect to f2 may then be taken from — « to 
+x and we get 


-2 ger tes cos 2r(fi — fe)r 
= d 1 dfoi 1)W\ fo = . 3.9-2! 
v= 4+ [oi f dhwinwn So pp OO 


The power spectrum W(f) of A(t) may be obtained by integrating V(7): 
Wf) = s| W(r) cos 2nfr dr 
0 


Let us concern ourselves with the fluctuating portion A(t) — A of AU 
Its power spectrum W’,(f) is 


Wf) =4 [ (W(r) — A’) cos 2xfr dr 
The integration is simplified by using Fourier’s integral formula in the form 
[ dr * dfe F( fe) cos 2x(u — fo)r = 3F(u) 
We get 


1 @ 
Wf) : [ dflw( fw f + fr) + w(fdw(—f + fl 


_ o + 4a’ f” 


1 +00 
~ at + dn’ f [ w(fiw if — fi)dfi 


(3.9-30 


The simplicity of this result suggests that a simpler derivation may be 


found. If we attempt to use the result 
2S E 02.5-3 
T 


a(f) = Limit 


T—2 
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e. from where S(f) is given by (2.1—2) we find that we need the result 


D) yy T ae: a P 
Limit f ar, | KPO FOF) 
T—0 0 0 
3. 9-28 te (3.9-31) 
=[  wAw(f - Aah 


ectrum 
€ asin Fe where f > 0 and /(t) is a noise current with w(f) as its power spectrum. 
This may be proved by using (3.9-7) and 
ao +2 
s | ¥(r) cos 2nfr dr = / w(x)w(f — x) dx 
0 h— 20 
fo)t which is given by equation (4C-6) in Appendix 4C. 
- fo) |? An expression for the mean square fluctuation of A(f) in terms of w(f) may 
yual t be obtained by setting 7 equal to zero in (3.9-29) 
—x to (A(t) — 4)? = VO) — A 
[afa4>" (3.9-32) 
= ay Off 2 5 4-%e 42 
3. 9-9 0 ~ a 2 dr (fi Ts) 


The same result may be obtained by integrating W(f), (3.9-30), from 0 
v(r): fe tow: 
: .« 2:2. 
i oe a anf? * dfiw(fiw(f — fi) (3.9-33) 
f A 5 Although this differs in appearance from (3.9-32) it may be transformed 
into that expression by making use of w(—/f) = w/(f). 
Suppose that /(¢) is the current through an ideal band pass filter so that 








) w(f) is zero except in the band f, < f < fy where it is wo. Then, if 3f. > fi, 
ie form (4 A=— (fo — fa) (3.9-34) 
E a 
4 mik~j,~f) C<fsh~y 
: | w(x)w(f — x) dx = wo(f — 2fa) Fas fsfotfa 


i t | wo(2h — f) fotfasf<e xn 


and is zero outside these ranges. The power spectrum HW’.(f) may be ob- 
3.9-30) P P ‘ ° er . .~ 2-2 
tained immediately from (3.9-30) by dividing these values by a + 4n°f’. 


From (3.9-33) 


a® + 41° f* 


fotSa (f — 2f,) > fre (%, —~ f) 
e «ja ] Wo | = _ 
Pe a 


may be 4 (A(t) = A)’ = 2u’g 


2 : (fa fo fF 
0 


df 


forte a + 42 /? 





(2.5-3) | + wo 
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Ii an exact answer is desired the integrations may be performed. When we 
assume that fy — fa << fy + fa we may obtain approximations for the last 
two integrals. 
el | fuw~ fe 1 2x( fs — fa) 
tan 
TH 


(A(t) — A)” = w 
Qa 


oe ON 


4dr - ( 


a” a* + 4 (fs + fa) 


a. (f, — fa) | 


Furthermore, if 2a(fp — fa)/a@ is large we have 


(A(t) = Ay = Wp 


and the relative r.m.s. fluctuation is 


ee Ee my | a i 
.m.s. 0 = = 
A 2( fo — fa) 


This result may also be obtained from (3.9-10) and (3.9-11) by assuming 
a so small that the integral for A(t) may be broken into a great many in- 
tegrals each extending over an interval 7. a7 is assumed so small that 


au 


is substantially constant over each interval. 
3.10 DistTRIBUTION OF NorIsE PLuUs SINE WAVE 
Suppose we have a steady sinusoidal current 

I, = I,(t) = P cos (wot — ¢p) (3.10-1 
We pick times f; , f2, «++ at random and note the corresponding values of 
the current. How are these values distributed? Picking the times at ran- 
dom in (3.10-1) is the same, statistically, as holding ¢ constant and picking 
the phase angles ¢, at random from the range 0 to 2x. If J, be regarded as 
a random variable defined by the random variable ¢,, its characteristic 
function is 


ave ete ee bf eg ene cos (wyt—¢) de 
dr Jo (3.10-2) 
- J (Ps) 
and its probability density is 
toe 1 2-1/2 
: e ''» Jo(Ps) dz = [La T,) in| <P (3.10-3) 
ogi: ( 0 \I,|>P 


In this case it is simpler to obtain the probability density directly from 
(3.10-1) instead of from the characteristic function. 
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Now suppose that we have a noise current /y plus a sine wave. By com- 
bining our representation (2.8-6) for Jy with the idea of ¢, being random 
mentioned above we are led to the representation 


It) =Il=1,+ 1s 
VW 
= P cos (wot — ¢p) + Zz, Cn COS (wal — gn),  (3.10-4) 
1 


cy = Qw(fn)af 


where ¢, and ¢1, +++ gw are independent random angles. 

If we note J at the random times ¢; , fo --+ how are the observed values 
distributed? Since 7, and Jy may be regarded as independent random 
variables and since the characteristic function for the sum of two such vari- 
ables is the product of their characteristic functions we have from (3.1-6) 


and (3.10-2) 


iz i2(J,+ 
ave. e’ = ave, 777’ 
— yor" (3.10-5) 
= Jo(Pz) exp 4 
which gives the characteristic function of 7. The probability density of J 
is”. 

1 +2 1 r 
~tz Yoz2/2 . I—P cos 6)2/2y > \ 
has Jo( Ps) ds = I € cov" do (3.10-6) 

at Leo TV lr 0 


In the same way the two-dimensional probability density of (J; , J), 
where J; = /(t) is a sine wave plus noise (3.10-4) and J; = I(t + 7+) is its 
value at a constant interval 7 later, may be shown to be 


2 2-1/2 ale 
(Yi ni i d@ exp ee | (3.10-7) 
where 
B(6) = ¥ol(i — P cos 0)” + (I, — P cos (8 + w,7))’] 
— WI, — P cos 6)(I. — P cos (8 + w,7)) 


The characteristic function for J; and Jo is 





july+ivlg a) 9 
ave, eM? = Jy Py/i? + o + uv cos wp) 


b Si ? (3.10-8) 
x ep| -4 (w +e) v0 


7 A different derivation of this expression is given by W. R. Bennett, Jour. Acous. Soc. 
Amer., Vol. 15, p. 165 (Jan. 1944); B.S.T.J., Vol. 23, p. 97 (Jan. 1944 
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Sometimes the distribution of the envelope of 
I = P cos pt + In (3.10--9 


is of interest. Here we have replaced w, by p and have set ¢, to zero. by 
the envelope we mean R(t) given by 


R(t)=R=(P+/)/4+ 6 (3.10-10 


where J, is the component of Jy “in phase” with cos pt and J, is the com- 
ponent ‘in phase” with sin pf: 


I. ys Cn COS [(wn — p)t — gal 
I; ms Cn Sin [(wn — p)t — gal 
Iy = I, cos pt — TI, sin pt 


r = r = r = Yo 
Since 7, and /, are distributed normally about zero with a variance of 
Yo, the probability densities of the variables 
x= fp + ie 
y= 1, 
are 


—1/2 _ (« — P)’ 
(2rYo) "exp The 





(2mo) 7"? exp — 2 


240 
respectively. Setting 
x = Rcosé 
y = Rsin@ 
and using these distributions shows that the probability of a point (x, y) 
lying in the ring R, R + dR is 
RdR [* 


am J exp | - 3, (R° ob P” — 2RP cos 0| do 


_RdR, [ R+P RP 
mi ep| Wo ]n(*?) (3.10-11) 


where Jo is the Bessel function with imaginary argument. 
CJ 2n 


Io(z) = 


nao 22" n!n! 
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and is a tabulated function. Thus (3.10-11) gives the probability density 
of the envelope R. 

The average value of R" may be obtained by multiplying (3.10-11) by R" 
and integrating from 0 to x. Expansion of the Bessel function and term- 
wise integration gives 


7 niraW —p? 2yo . P 
R" = (2p) I (; + 1)e * AS + 1;1; a 


= (2yo)” *¥ (° + I) (-%, 1 Ir -*) (3.10—12) 


+ . = . 38 . - ~ . 
where ,F; is a hypergeometric function.” In going from the first line to 
the second we have used Kummer’s first transformation of this function. 
A special case is 


R? = P? + Ao (3.10-13) 
When only noise is present, P = 0 and 


1/2 
(2¥o)'* P(g) = (“0”) 


R? = WY 


R 


(3.10—14) 


Before going further with (3.10-11) it is convenient to make the following 
change of notation 
R dR 7 
v = 1/2? dv = 1/2? a - 1/2 
Yo 0 Yo 

is the ratio (sine wave amplitude )/(r.m.s. noise current). 
Instead of the random variable R we now have the random variable 7 whose 
probability density is 


(3.10-15) 


oe? 
a 


p(v) = vexp | - . : | Io(av) (3.10-16) 


Curves of p(v) versus v are plotted in Fig. 6 for the values 0, 1, 2,3, 5 of a. 
Curves showing the probability that v is less than a stated amount, i.e., dis- 
tribution curves for 2, are given in Fig. 7. These curves were obtained by 
integrating p(v) numerically. The following useful expression for this 
probability has been given by W. R. Bennett in some unpublished work. 


I p(u) du = exp| —* = | a () T,(av) (3.10-17) 


« n=l 


% Curves of this function are given in ““Tables of Functions”, Jahnke and Emde (1938), 
p. 275, and some of its properties are stated in Appendix 4C. 
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a 
This is obtained by integration by parts using 


fw T,-1\(au) du = u"T,(au)/a 


When az >> 1 but 1 << a — 2, Bennett has shown that (3.10—17 


leads to 
Poe x “2: ae. v— a) 
I p(u) du = a ee 


A : (3.10-18) 
(1 aa 3(a +)" _ 4u" : 


Sava — v)° 








Fig. 6—Probability density of envelope R of I(t) = P cos pt + Ty 


This formula may also be obtained by putting the asymptotic expansion 
(3.10-19) for p(v) in (3.10-17), integrating by parts twice, and neglecting 
higher order terms. 

When av becomes large we may replace /o(av) by its asymptotic expres- 
sion. The expression for p(zv) is then 


' 1 » \? (v — a) ; 
v) ~ = 2X _ eee 3.10-19) 
P(r) (1 + te (=) exp| 5 | (3.10-1 


Thus when either a becomes large or 7 is far out on the tail of the probability 
density curve, the distribution behaves like a normal law. In terms of the 
original quantities, the normal law has an average of P and a standard devia- 
tion of Yj °. This standard deviation is the same as the standard deviation 
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of the instantaneous values of Jy. When av >1anda> »—a we may 
expand the coefficient of the exponential term in (3.10-19) in powers of 








Fig. 7—Distribution function of envelope R of I(t) = P cos pt + Ty 


(v—a)/a. Integrating this expansion termwise gives, when terms of magni- 
tude less than a ” are neglected, 


1 v-—a 


‘ 1 
| p(u) du = 5 + 3 erf V3 


we ' g—a i+ @¢=— @} (v — a)? 
lan 2e fa bg 8a” aa 2 
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When J consists of two sine waves plus noise 

T= Pcospt+OQOsingt+TIy, (3.10—21) 
where the radian frequencies p and g are incommensurable, the probability 
density of the envelope R is 

ao 

R | rJ (Rr) Jo( Pr) Jo(Onye °° dr (3.10—21 

0 
where Yo is J). When Q is zero the integral may be evaluated to give 
(3.10-11). When both P and Q are zero the probability density for R 
when only noise is present is obtained. If there are three sine waves instead 
of two then another Bessel function must be placed in the integrand, and 
soon. To define R it is convenient to think of the noise as being confined 
to a relatively narrow band and the frequencies of the sine waves lying 
within, or close to, this band. As in equations (3.7—2) to (3.7-4), we refer 
all terms to a representative mid-band frequency f,, = w»/2« by using 
equations of the type 

cos pl = cos \(p — wnt + wnt | 
= COS (P — Wm)t COS Wnt — Sin (Pp — wm)t SIN wl. 
In this way we obtain 
V = Acosw,t — Bsinw,¢ = R cos (wt + 6) (3.10—22 

where 4 and B are relatively slowly varying functions of ¢ given by 
A = P cos (p — w»)t + Q cos (g — Wy )t 


oa z Cr COS (wnt — wmt — on) 
: (3.10-23 
B = P sin (p — om»)t + Q sin (¢ — wade 


a : i Cr SiN (wrt — wml — gn) 


and 


R® = A? + B* R>0O 
(3.10-24 
tan @ = B/A 


As might be expected, (3.10-21) is closely associated with the problem 
‘ < . > 3 : 
of random flights and may be obtained from Kluyver’s result® by assuming 


38 G. N. Watson, “‘Theory of Bessel Functions” (Cambridge, 1922), p. 420. 
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the noise to correspond to a very large number of very small random dis- 
placements, 

Another way of deriving (3.10-21) is to assume (Pp — wall, (g — wy it, 
C1, ¢2 .-- tare independent random angles. The characteristic function 


’ 


of A, Bis 
juA+ivB 9 ° > . 
ave, 478 = J PV 8 + 2) J(OV 102 + ve 


The probability density of A, B is 


1 2 +O ~x 

—inA-—tvB uArivl 
(+) / au | Ce eae 
ais L— 26 — 30 


When the change of variables 
A = Reos@ u=rcose 
B= Rsin@ v =rsing 


is made the integration with respect to ¢ may be performed. The double 
integral becomes 


of rJo( Pr) Jo(Or)Jo( Re #8!” de 
Qn 9 


This leads directly to (3.10-21) when we observe that d1dB = RdRdé. 
Incidentally, if 
IT = Q(1 + k cos pt) cos gt + Ty 
in which p << gq, similar considerations show that the probability density 


of Ris 


a" ae 





R Qn @ 
— I da | rJo(Rr)JolQOr(1 + k& cos a) 
2m Jo 0 


when wm is taken to be g. The integration with respect to r may be per- 
formed. This relation is closely connected with (3.10-11). 

Returning now to the case in which / is the sum of two sine waves plus 
noise, we may show from (3.10-21) and 





ir (1 - *) 
@ 4 i * 
[ jp) ah « —— 
0 


(3) 
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that the average value of R" is, when —2 < re (n) < — 3, 
+ n 
~ I (1 ee ) p®@ 
R= ——— | r" "Jo Pr)J(Orye *"* dr 
i nN 0 3 
od a. 
) 
i k \n > - 
<ietixs eie (-") (—x)"(—y)" — (3.10-25 
n/2 i , kim 
= (2y)""r +1 Ss ee! Las io 
v (5 ip ns kik!m!m! 
n 
= a iy = x)* ; 
psn (* 4a a ee 
= iy) EF + ae k = 
y k=0 Rik! t-, 
It appears very probable that this result could be extended, by analyti 
continuation, to positive integer values of 7. We have used the notatio; 
(a)o = 1, (ay = afa +1) +++ (a +k — 1) 
| ‘Ow (3.10-26 
C= ‘ y= s 
2 2Wo 
and have denoted the Legendre polynomial by P;(2). The series converg: 
for all values of P, QO, and Yoand terminate when 7 is an even positive integer. 
When x or y, or both, are large in comparison with unity we may use the 
| integral for R" to obtain the asymptotic expansion, assuming Q < P so 
that.) <9; 
n n 
Rw py NAAN th (pt pity 4 (3.10-27 
ssh ‘ ' : 24") Sees pia ’ } 3. "Ll 
imo k! x Z Z x 
When # is an even positive integer this series terminates and gives the same 
expression as (3.10-25). When wis an odd integer the o/*; may be expresse« 
in terms of the complete elliptic functions E and A of modulus y' *x7' 
| ' 
y iy ae OE we 
oF -1,-11;2) = E- 1—-jK 
+ 7 T x 
(3.10-28 
: v = 
(4, 33 1;2) =-K 
x T 
The higher terms may be computed from 
a(1l — 3) oF (a +1,a+1;1;s) = (2a — 1)(1 + s)oFi(a, a; 1:2) 
+ (1 — a)oF\(a — 1,a — 13132) (3.10—29 
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which is a special case of 


ably + 1) — z “Fila + 1,6 + 13¢;2) = AoFi(a, bc: 2 
—(y—1)(c—a)(c—)b F(a — 1,6 —1:6¢:2 3.10-—30 
where y = ¢ — a@ — band 


A=(¥ -— Dy t 1 — aly — Dle — b)(b— 1) + (y + Dale —a - 1 
Although this expression does not show it, A is really symmetrical in a 
and 6. Asymmetrical form may be obtained by using the expression ob- 
tained by putting s = 0 in (3.10-30). 


3.11 SHot Errect REPRESENTATION 
In most of the work in this part the representations (2.8-1) or (2.8-6 
have been used as a starting point. Here we point out that the shot effect 
representation used in Part I may also be used as a starting point. 
For example, suppose we wish to find the two dimensional distribution of 
I(t) and [(t + 7) discussed in Section 3.2. This is a special case of the distri- 
bution of the two variables 


+” 
> Fit — t) 


k=—2o 


I(t) 
a (3.11-1) 
Jit) = > Gt-& 
h=—oo 


where we now assume 


| F(t) dt = | G(t) dt = 0 (3.11-2 


in order that the average values of J and J may be zero. In fact, to get 
I(t + 7) from J(t) we set G(¢) equal to F(t + 7). 

The distribution of J and J may be obtained in much the same manner 
as was the distribution of J alone in section 1.4. The characteristic func- 
tion of the distribution is 


- bu +ivd 
t(u, v) = ave. ” ill 


cas 
= exp v / fern Perreea) _ 1] dt 
I 30 


where v is the expected number of events (electron arrivals in the shot effect 
per second. The probability density of J and J is 


1 +2 +20 : 
4 :[ au | dve “"*™ f(u, ») (3.11-4 
T Jax — 20 
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The semi-invariants X,,,, are given by the generating function 
“ 

aie cate . in penuh 

™" (iu)™(iv)" + of(iu)*, (iv)*] 


log f(u, v) = 
aia mnai mn! 


and are 
— [ F"(t)G"(t) dt (3.11-5) 


Asy— « the distribution of J and J approaches a two dimensional norma! 
law. The approximation to this normal law may be obtained in much the 
same manner as in section 1.6. From our assumption (3.11-2) it follows 
that Ay and Aq, are zero. From the relation between the second moments 
and semi-invariants A we have 


+2 


Ao» + Ajo = rf F°(t) dt 


2 


Mil 


+2 


= Aun + AwAn = V / F(t)G(t) dt (3.11-6 


A) 


= 
| 


wr = No + An = » | G*(t) dt 


where the notation in the subscripts of the y’s differs from that of the 2’s, 
the change being made to bring it in line with sections 2.9 and 2.10 so that 
we may write down the normal distribution at once. 

The formulas (3.11-6) are closely related to Rowland’s generalization of 
Campbell’s theorem mentioned just below equation (1.5-9). 
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PART IV 
NOISE THROUGH NON-LINEAR DEVICES 


4.0 INTRODUCTION 


We shall consider two problems which concern noise passing through 
detectors or other non-linear devices. The first deals with the statistical 
properties of the output of a non-linear device, that is, with its average 
value, its fluctuation about this average and so on. The second problem 
may be stated more definitely: Given a non-linear device and an input 
consisting of noise alone, or of noise plus a signal. What is the power 
spectrum of the output? 

There does not seem to be much published material on the first problem. 
However, from conversation with other people, I have learned that it has 
been studied independently by several investigators. The same is probably 
true of the second problem although here the published material is somewhat 
more plentiful. This makes it difhicult to assign credit where credit is due. 
Much of the material given here had its origin in discussions with friends, 
especially with W. R. Bennett, J. H. Van Vleck, and David Middleton. 
Help was obtained from the recent paper” by Bennett, and also from the 
manuscript of a forthcoming paper by Middleton.” 


4.1 Low FREQUENCY OvUT?UT OF A SQUARE LAW Devict 
Let the output current J of the device be related to the input voltage 1 by 


where @ is a constant. When the power spectrum of V is confined to a 
relatively narrow band, the power spectrum of 7 consists of two portions. 
One portion clusters around twice the mid-band frequency of Vand the 
other around zero frequency. We are interested in the low frequency 
portion. The current corresponding to this portion will be denoted by 
1¢, and is the current which would flow if a low pass filter were inserted 
in the output to remove the upper portion of the spectrum. It is convenient 
to divide /;¢ into two components: 


Tut = Tae + Its (4.1-2) 


37 Loc. cit. (Section 3.10). 

40 Cruft Laboratory and the Research Laboratory of Physics, Harvard University, 
Cambridge, Mass. In the following sections references to Bennett’s paper and Middle- 
ton’s manuscript are made by simply giving the authors’ names. 
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where the subscripts stand for ‘total low” frequency, “direct curre: 
and “low frequency,”’ respectively. 


Tac = average 7¢= [yt 


2 72 

Mean Square /¢; average ([,¢ — Tg.)° = [i¢—- 
Probably the simplest method of obtaining J/g, is to square the given ex- 

pression for V and pick out the terms independent of time. 


P cos pl + © cos gt + Vy 


It; may also be obtained by picking out the low frequency terms. 
ever, here we wish to use the square law device, and the linear rectifier in the 
next section, to illustrate a general method of dealing with the statistica 
properties of the output of a non-linear device when the input voltage is 
restricted to a relatively narrow band. 


If none of the low frequency spectrum is removed by filters, 


where R is the envelope of 1. The probability density and the statistical] 
properties of Z,¢ may be derived from this relation when the distribution 
° . - . 41 
function cf R is known. 


Before discussing these properties we shal! 
establish (4.1-6). 
Equation (4.1-6) is a special case of a more general result established 


in Section 4.3. However, its truth may be seen by taking the example 


= P cos pl a a) COs qt + Vy 
where f, = p, 2m and f, = gq, 27 lie within, or close to, the band of the noise 
By using formulas of the type 


cos [(p — wm)t + wnt | 


b 
i 
fe 
ef 
4 
ie 
3g 
ey 
i 
| 
: 
h 


| 


COS (PP — Wm)t COS Wnt — sin (Pp — wm )t Sin wnt 





4 When part of the low-frequency spectrum is removed, the problem becomes much 
more difficult. Jy. may be obtained as above, but to get I, it is necessary to first deter- 
mine the power spectrum of J (Section 4.5) and then integrate over the appropriate por- 

Concerning the distribution of I¢; , our present knowledge tells us only that it 
lies between the one given by (4.1-6) and the normal law which it approaches when only 


ci nS 


a narrow portion of the low frequency spectrum is passed by the audio frequency filter 
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ve may refer all terms to the mid-band frequency f,, = w, 22, as is done 
in equations (3.7—2) to (3.7-4). 
In this way we obtain 
V = A cos wat — Bsin wat = Roos (wat + 8), 41-8 


where A and B are relatively slowly varying functiens of ¢ given by 
A = P cos (p — am)t + QV cos (¢ — wn)t + Du Cn COS (wnt — wnt — on), 
B = Psin (p — wm)t + Q sin (¢ — wnt + 3 C, SIN (w,t — onl —¢ 


and 


tan @ = B/S A. 


This definition of R has also been given in equations (3.10 
The envelope of V is R and the output current is 


I = aR : a : cos (2w,f + 2 | (4.1-10) 
Since R is a slowly varying function of time, so is R°. The power spectrum 
of R° is confined to frequencies much lower than 2/,, and consequently the 
power spectrum of R° cos (2w»f + 26) is clustered around 2f,,. Thus the 
only term in J contributing to the low frequency output is aR” 2 which is 
what we wished to show. 

We now return to the statistical properties of /.¢. First, consider the 
case in which V' consists of noise only, 1) = Vy, so that the probability 


density of the envelope R is 


TP atin (3.7-10) 
Yo 
where 
Yo = [rms Vx’ = V4 (4.1-11) 
Hence 
= ak? 
Fie bis Let —_ > 
@o 2 
ak R y 
«| Sora 
4 2 y 
= avo 
— ~~ aw oR’ ee ; 2 
y= Tete = | S oP aR — te 
“0 4Yo (4.1-12) 
= ayo 
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Second, consider the case in which 





Vv = Vy + P cos pt (4.1-13 


where p/ 27 lies near the noise band of Vy. The probability density of the 


envelope R is 


i 
} 


Fool ]n(8) ems 


From this and equations (3.10-12), (3.10-13), we find 
ak? , oP” 
Lac = > = ayo 7° (4.1-14) 


rs <= 


2 . a P' 
y= 7 R=a | 2 + 2P*yo + | 


Ty = Tt — Tie = alo + P’lo (4.1-15 


In (4.1-14) Yo is the mean square value of Vy and P*/2 is the mean 
square value of the signal. These two equations show that Ja. and the 
rms value of J; are independent of the distribution of the noise power 
spectrum in ly as long as the input V is confined to a relatively narrow band. 
In other words, although this distribution does affect the power spectrum 
of the output, it does not affect the d.c. and rms J7, when Yoand P are given. 
That the same is also true for a large class of non-linear devices was first 
pointed out by Middleton (see end of Section 4.9), 

When the voltage is” 


V = Vy + P cos pt + Q cos qt, (4.1-4) 


p ~ q, we obtain from equation (3.10—25) 


2 ? 

lee = % R = «(vo + 4 +) 

Tit = aR (4.1-16) 
Ii, = at vi + Po + Oyo + rey 


® These results are special cases, obtained by assuming no audio frequency filter, of 
formulas given by F. C. Williams, Jour. Inst. of E. E., 80 (1937), 218-226. Williams also 
discusses the response of a linear rectifier to (4.1-4) when P >> Q + Vy. An account 
of Williams’ work is given by E. B. Moullin, “Spontaneous Fluctuations of Voltage,” 
Oxford (1938), Chap. 7. 
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4.2 Low Frequency Output oF A LINEAR RECTIFIER 


13 In the case of the linear rectifier 
the 0, : 8 
Bas (4.2-1) | 
lal V>0 
i 
1] the low frequency output current, assuming no audio frequency filter, is 
aR 
lu = — (4.2-2) 
7 
4) This formula, like its analogue (4.1-6) for the square law device, assumes 


that the applied signal and noise lie within a relatively narrow band. It 
may be used to compute the probability density and statistical properties 
of 7,¢ when the corresponding information regarding the envelope R of the 
applied voltage is known. 


5 med ‘ ee 
‘ The truth of (4.2-2) may be seen by considering the output 7. It con- 
an | sists of the positive halves of the oscillations of al’. The envelope of J is 
he the same as that of al’. However, the area under the loops of 7 is only about 
ver [E 1/7 of the area under aR, this being the ratio of the area under a loop of 
1d. » sin x to the area of a rectangle of unit height and length 27. From the 
im low frequency point of view these loops of J merge into a current which 
oy. (varies as aR/7z. 
rst ; When V is a sine wave plus noise, 
V = Vy + P cos pt (4.1-13) 
4 the average value of 7,4 is® 
1/2 2 
ao Yo " 1 P 
Tic = R=al\- iF, ;1;- 
ar 2 2Wo 
sia (4.2-3) 
Wo —") x x 
=al— e° |A1+ x10 + xl, 
2 2 2 
6) > where Jo, 7; are Bessel functions of imaginary argument and 
r ave. sine wave power (4.2-4 
x= = .2-4) 
‘ 2Wo ave. noise power 
4 
j * This result was discovered independently by several investigators, among whom we 
of » may mention W. R. Bennett and D. O. North. The latter has applied it to noise measure 
Iso > ment work. He has found that the diode detector, when adapted to noise metering, is a 
int great improvement over the thermocouple, and has used noise meters of this type satis- 
i tactorily since 1940. See D. O. North, “The Modification of Noise by Certain Non 
Linear Devices”, Paper read before I.R.E., Jan. 28, 1944. 
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Yo being the average value of Vy. Equation (4.2-3) follows from 
formulas (3.10-12) and (4B-9). When wx is large the asymptotic expans 
(4B-3) of the ,/; gives 


= a a P+ Yo ob Yo 4 | (4.2 : 


- 2P * SP 


Similarly, the mean square value of /,¢ is 


Ie = vet R’ aa an (P* + 2Wo) (4.2 6 
rT 


7 


and the mean square value of the low frequency current J¢; , excluding t) 
d.c., is given by 


fii, = Fe 2. 


When x is large we have 


"es ee _a@ al = 
Ti; = |v dp? | =— x Wo E Ax | (4.2 i 


and when x = QO, 
—— ¢ 2 + ’ 
lt; = = vo(2 i 4 (4.2 


Curves for Ja. are given in Figures 1, 2 and 3 of Bennett’s paper. IH 


s 


also gives curves, in Fig. 4, showing //; versus x. These show that t! 
effect of the higher order modulation terms is small when J/¢; is comput: 
by adding low frequency modulation products. 

When V consists of two sine waves plus noise, 


V = Vy + P cos pt + Q cos qt, (4.1-4 


the average value of J,¢ is, from (3.10-25), a sort of double :F; functio: 


e 1/2 @ 2 aA 
Tac = - R-a (#) Zz. z ( 2)k-+m (—x)*(—y)” 


2a] 0 mao kiki mim! a 
nd 
_  (w\'* S (—de jt p (*% +9 
i a(¥ p> i di liad, x—y 
where 
odie 4 , y= c , P,(s) = Legendre polynomial (4.210 
2Yo 2fo 


If x is large and vy < x, we have from (3.10-27) the asymptotic expressio 


bed poet pee | 
i = (—3)e(—3)e oF (: —i,k—}; 1;2) (4.2-11 


us k=0 ki xk 
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m the JB The of; may be expressed in terms of the complete elliptic functions E and 
' > 1/2, —1/2 an i 
uns | AK of modulus yx “~. Thus 
1 4. 2 3 i 
(1, -t1;2)=f2-2( -*\k, 
— i i \ ou | 
4.2-5 (3.10-28) 
\ Senge 
of ,, A | if = K j 
x 7 | 
| 
42 ind the higher terms may be computed from the recurrence relation | 
2-6 E., Ps : ae 
b (3.10-29). The first term, k = 0, in (4.2-11) gives Jz. when the noise is 
44 i 
absent. 
ng ; , ; 
The mean square value of /;¢ is 
ee ae ar 
he= R= —2b+tP +] (4.2-14) 
From this expression and our expression for Jg-, the rms value of the low 
42-7 frequency current, /¢; , excluding the d.c., may be computed. For example, 
rh when the noise is small, 
2 a » 2 > 1 1 y\\ | 
lig~ =| P +O — ( PoPil —3, —3; 1; 
. 7 x 
(4.2-§ (4.2-15) 
; v\ A 
+ 2Yol 1 — oft me, 95 15: ) | 
r. He & t} * 
ahd The term independent of Yo gives the mean square low frequency current 
| | 
mpule: F in the absence of noise. As Q goes to zero (4.2-15) approaches the leading 
term in (4.2-7), as it should. When P = Q our formula breaks down and 
. . 45 
it appears that we need the asymptotic behavior of 
(4.1-4 
= 1/2 \ 
, I Yo > (—3 (2k)! ( . k 
tio! de = a =a) 
nctl ? — [RI y* 
In view of the questionable nature of the derivation given in Section 3.10 
(4.2-9 f— of equations (4.2-9) and (4.2-11) it was thought that a numerical check on 
their equivalence would be worth while. Accordingly, the values x 1, 
y = 3 were used in the second series of (4.2-9), It was found that the 
| largest term (about 130) in the summation occurred at k = 11. In all, 24 
' terms were taken. The result obtained was 
4,2-10 R 
2.5502 


pressio V 2 


“See W. R. Bennett, B.S.T.J., Vol. 12 (1933), 228-243. 
*® This may be done by the method given by W. B. Ford, Asymptotic Developments, 
Univ. of Mich. Press (1936), Chap. VI. 


4.2-11 
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For the same values of x and vy the asymptotic series (4.2-11) gave 
240 + O171 + 075 + 0.52 + 


If we stop just before the smallest term we get 2.57 for the sum. If we 
include the smallest term we get 2.65. This agreement indicates that 
(4.2-11) is actually the asymptotic expansion of (4.2-9). 
When the voltage is of the form 
V = O11 + kos pt) cos gt + Vy 


we may use 


ie _ eri wi 


rl 1; —y(1 + cos 0] d9 


(4.2-16 


where R is the envelope with respect to the frequency g/2m and y is given 
by (4.2-10). The integral may be evaluated by writing ,/*°) as a power 
series and integrating termwise using the result 


(1 + k cos 0)f cos m6 dé 


(—(€) — _ _ 9 
_ Om (—k)” oF k a = : +3 ;m + 1; ‘| 


1 
dr Jo 
2" m! 2 


where m is a non-negative integer, f any number, 


(a), = ala + 1)-+- (a +m — 1), (a)o = 1, and (O)o = 1. 


The integral may also be evaluated in terms of the associated Legendre 


function, 
By applying the methods of Section 3.10 to (4.2-16) we are led to 


aed ; ke 
R= (144) +m 

a Se (4.2-18) 
R~o> vAws 


he 9 
ny 1 e +» bB 

Gi of i(s — 3, $31; R°) 

s=0 . 


where the asymptotic series holds when y is very large and & is not too close 
to unity. These expressions give 


i ~ "(0% + v2 - 0-8 + +) (4.2-19 
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The reader might be tempted to associate the coefficient of Yo in (4.2-19) 


with the continuous portion of the output power spectrum. However, this 


would not be correct. It appears that the principal contribution of the 
continuous portion of the power spectrum to /7; is a Wo/m, just as in (4.2-7) 
when k is zero, The difference between this and the corresponding term 
in (4.2-19) seems to arise from the fact that the amplitude of the recovered 
signal is not exactly aQk/w but is modified by the presence of the noise. 


This general type of behavior might be expected on physical grounds since 


changing P, say doubling it, in (4.2-7) does not appreciably affect the 7, 
in (4.2-7) (which is due entirely to the continuous portion of the noise 
spectrum). The modulating wave may be regarded as slowly making 


- changes of this sort in P. 


4.3. Some STATISTICAL PROPERTIES OF THE OUTPUT OF A GENERAL 
Non-LINEAR DEVICE 


Our general problem is this: Given a non-linear device whose output J is 
related to its input V’ by the relation 


f= = / F(iu)e’’“ du (4A-1) 
£7 Jo 


which is discussed in Appendix 4A. Let the input V contain noise in addi- 
tion to the signal. Choose some frequency band in the output for study. 
What are the statistical properties of the current flowing in this band? 

It seems to be difficult to handle this general problem. However, it 


| appears that the two following results are true. 


1. As the output band is chosen narrower and narrower the statistical 


- properties of the corresponding current approach those of the random noise 
- current discussed in Part III (provided no signal harmonic lies within the 
> band). In particular, the instantaneous current values are distributed 
» normally. 


2. When the input V is confined to a relatively narrow band the power 
spectrum of the output / is clustered around the 0 (d.c.), Ist, 2nd, ete. 
harmonics of the midband frequency of V. The low frequency output in- 
cluding the d.c. is 


Ii¢@ = Ao(R) = = | F(iu)Jo(uR) du (4.3-11) 
Tv ¢ 


where R is the envelope of V. 
The envelope of the xth harmonic of the output, when x > 0, is 


A,(R) = S| F(iu)J,(uR) du (4.3-1) 
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The mathematical statement is 


co] 


I = > A,(R) cos (twnt + n6) (4.3~9 


n= 


where fn = @m/(2r) is the representative mid-band frequency of Vand 4 


is a relatively slowly varying phase angle. The results of Sections 4. 
and 4.2 are special cases of this. 


Middleton’s result that the noise power in each of the output bands 


the entire band corresponding to a given harmonic) depends only on Vy 
Yoand not on the spectrum of Vy , where Vy is the noise voltage componen! 


of V, may also be obtained from (4.3-9). We note that the total power 


. h i 

in the x“* band depends only on the mean square value of its envelop: 
A,(R), and that the probability density of the envelope R of the input i: 
volves Vy only through yo. 


The argument we shall use in discussing the first result is not very satis- 


factory. Itrunsas follows. The output current J may be divided into tw: 


parts. One consists of sinusoidal terms due to the signal. The other con- 


sists of noise. We shall be concerned only with the latter which we sha! 
call Jy. The correlation between two values of Jy separated by an interva! 
of time approaches zero as the interval becomes large. Let 7 be an interva 
long enough to ensure that the two values of Jy are substantia!|\ 


independent. Choose an interval of time T long enough to contain many 


intervals of length 7. Expand Jy as a Fourier series over this interval! 
We have 


dp < drut . 2rnt 
ly =~ + > 1a, cos —— + 6, sin — 
2 need 1 7 
) r -s 
a, — 16, = = ” eile *Tx(t) dt 
T Jo 
e 


B 
Let the band chosen for study be fo — : to fo+ , and let 


~ 


L (i ni | Pas .: (i > ‘) = te (4.3-3 


where 2, and #2 are integers. The number of components in the band is 
(tg — m,). We suppose 2 is such that this is small in comparison with 7) 7. 
The output of the band is 


“c darn Oar 
Jn = 2 E cos = t + 6, sin | (4.3-4 


$e: aan aaah Sire ater 21d ES, 
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ynere 
) r 
a, = ib, = ea [ 7 atte . mai? ee Fit) dt 
T to 
(4.3-5) 
Ny + No Ny + Ne =e ee 
n=— +n—" = foT + (n — fT) 
2 2 P 
We choose the band so narrow that 
NMg—nyK T/r or BrK1 } 3-6) 


This enables us to write approximately 


rt - - 
tla ((n/T7 foirr = —ilrfot 7 ; 1 
¢ T é L[y(l) dl 


T=] 


Gy — 1b, = 
an integer. we do this for 
Then [y(t 


1 gives us an ar 


r, = T, 7, T being chosen to make 1 Suppose 


a large number of intervals of length 7. will differ from interval 


tointerval. The set of integrals for r 


rav of values which 


we regard as defining the distribution of a complex random variable, say 


r;. Similarly the set of integrals for r 2 defines the distribution of a 


second random variable xs, and so on to x, Because we have chosen 7+ 


! 


so large that Jy(¢) in any one integral is practically independent of its values 


in the other integrals we may say that x1, ¥2 are independent. 
We have 
TI 
. ‘ —ile((n/T)—Jfo)rr .. 
d,, — 16,, = ¢ Lr 
r=] 
ri 
‘I ys SY mtr ((my + h/7 Sort 
Qni41 — Wasi = ¢ X; 
r=] 
3 Z : 
. —42r((n 7 fo)r 
an, — ib,, = Dre ” Be 
r=] 


and if m2 — my <“<9r,,as was assumed in (4.3-6), we may apply the central 


limit theorem to show that d,, , bn, , @nj41. 77° Ans, 4, tend to become in 


ae 
dependent and normally distributed about zero as we let the band width 
3— Oand T — = (and hence r; — ~) in such a way a 
In this work we make use of the fact that /y(¢ 


and imaginary parts of 41,2, °° > 


s to keep He — 
h that the real 


is Sue 


fixed. 
x, all have the same average and standard 
deviation. It is convenient to assume fo7’ is an integer. 

Thus as the band width 3 approaches zero the band output Jy given by 
(4.3-4) may be represénted in the same way, namely as (2.8 1), as was the 


random noise current studied in Part III]. Hence J, tends to have the 
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same properties as the random noise current studied there. For example 
the distribution of Jy tends towards a normal law. In our discussion 
had to assume that 67 <1. If the voltage V applied to the non-linear 
device is confined to a relatively narrow frequency band, say f, — f, . i: 
appears that the interval 7 (chosen above so that /(¢) and I(t + 7) are sub- 
stantially independent) may be taken to be of the order of 1/(f, — 

In this case Jy tends to behave like a random noise current if 8/ (fy — / 
much smaller than unity. 

We now turn our attention to the second statement made at the begin 
ning of this section. Let the applied voltage be confined to a relativel) 
narrow band so that it may be represented by equation (4.1-8) of Sectio: 
4.1, 


V = R cos (wmt + 6), R > OQ, (4.1-8 
where fm = Wm/(2m) is some representative frequency within the band 
and R and @ are functions of time which vary slowly in comparison wit) 
COS Wnt. We call R the envelope of V. 

From equation (4A-1) 


1 sit al (wmt+9) — 
I= 4 F(iue*® °°? “m** dy (4.3-7 
2dr Je 
We expand the integrand by means of 
eo) 
hoe — FS 6. i” cos neJ a(x) (4.3-8 
n=0 


where ¢€ is 1 and ¢, is 2 when » > 0 and J,(x) is a Bessel function. 
Thus 


IT = > A,(R) cos (name + 6) (4.3-9 
n=O 
where 
An(R) = 65 5 [ F(iu)Jn(uR) du (4.3-10 
Tc 


Since R is a relatively slowly varying function of time we expect the 
same to be true of 1,(R), at least for moderately small values of 7. Thus 
from (4.3-9) we see that the power spectrum of J will consist of a suc- 
cession of bands, the 2‘” band being clustered around the frequency n/,.. 
If we eliminate all of the bands except the 2‘ by means of a filter we 
see that the output will have the envelope 4,(R) when » = 1. Taking 
n to be zero, shows that the low frequency output is simply 


A(R) = = | -P(iu)Jo(uR) du (4.3-11 
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Taking to be one shows that the band around f,, is given by 


A(R) ,, 


R (4.3-12) 


The statistical properties of the low frequency output and of the en- 


velopes of the output bands may be obtained from those of R. For ex- 
ample, the probability density of 1,(.R) is of the form 


' _ J dA,(R) 
) 3-45 
pir) / 7: (4.3-13) 


| 
}where p(R) is the probability density of R. In this expression R is con- 
a 





sidered as a function of A, . 

+ It should be noted that we have been assuming that all of the band 
/ surrounding the harmonic frequency xf, is taken. When we take only a 
| portion of it, presumably the statistical properties will tend to approach 
. those of a random noise current in accordance with the first statement made 
/ at the beginning of this section. 

| When we apply (4.3-11) to the square law device we have 





da 
F(a) = — 
(4) (iu)® 
B) (0+) J >) 
Ao(R) = - ae 1 du 
2m ui 
a 2 
= GR 
When we apply (4.3-11) to the linear rectifier: 
F(iu) = — ~ 
: ie : ) 
; Ag(R) = -2[ OO ip 
é 2r Le u" 7 


' where the path of integration passes under the origin. These two results 


Ss 
* agree with those obtained in Section 4.1 and 4.2 from simple considerations. 


As a final example we find the low frequency output of a biased linear 


_ rectifier in terms of the envelope R of the applied voltage. From the table 
+ of F(iu) given in Appendix 4A we see that / (iu) corresponding to 


I =0, V<B 


I=V-B, V>B 
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—iuB 


F(iu) = —“— 


u? 
Consequently, the low frequency output is 


+00 
Ao(R) = ae [ e "8 T\(uR)u du 

2m Le 

where the path of integration is indented downwards at the origin. Whe: 
B > R the value of the integral is zero since then the path of integratio; 
may be closed in the lower half plane by an infinite semi-circle 
also follows at once from the physics of the problem. 
we may integrate by parts and get 


1 a x 


19(R) = ; | e"lGBIo(uR) + RJ (uR) \u ‘du 
ot J—a 
B aoe 
= . | iB sin uBJo(uR) + R cos uBsJ,(uR)\u "du 
& T 0 
(4.3-14 
B B 


= oe ar sin B ot. : VR? fer b 
) +r R 1 


raw 
ae 


4 


his hypergeometric function turns up again in equation (4.7-6), 
in the range —-R< B< R, 


dAg n 1 ‘ ta bB 
dR ot k? 


Als: 


When B is negative and R < —B, the path of integration may be closed 
by an infinite semicircle in the upper half plane and the value of the integra! 


is proportional to the residue of the pole at the origin: 


= —B 


Ao(R) 


Thus, to summarize, the low frequency output for our linear rectifier is, 


for B > 0, (R is always positive) 


A(R) = 0, R<B 

} 3 > a (4.3-15 
Ao(R) = we + : are sin : + VR — B, B<R 

es T R 7 


This value 


When —R < B < k 
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and for B < O it is 
f A(R) = |BI, R</|B| 
4 : 
{ |B| , |B beet 4 : (4.3-16) 
§ A(R) = +- ly! arcsin’—'+- VR? — BY |B < R 
| 2 Tv R T 
where the are sines lie between 0 and w/2.  Ao(R) and its first derivative 
with respect to R are continuous. 
When fF) ~~ From (4.3-15), the d.c. output current is, for B > 0, 
Tat i 2 
Tt ee Pe | ; 
Value I. =| —>~ + —arcsin—= + - VR? — B*| p(R) dR (4.3-15) 
7 B 2 us R us 
>< K 
where p(R) is the probability density of the envelope of the input V, e.g., 
pR) is of the form (3.7-10) for noise alone, and of the form (3.10-11) for 
noise plus a sine wave. Similarly, the rms value of the low frequency 
current Ty; , excluding d.c., may be computed from 
3-14 ea 
where, if B > 0, 
~. . SZ - ——_._ | 
lve = { —— + —arcsin— + - VR? — B*| p(R) dR (4.3-16) 
B p. 7 R us ° 
If V consists of a sine wave of amplitude P plus noise Vy , so it may be 
I I . 
As represented as (4.1-13), and if P > rms Vy, the distribution of R is 
approximately normal. If, in addition, P — B > rms Vy > 0, (4.3-15), 
(4.3-16), and (3.10-19) Jead to the approximations 
ae oe eee J 
la = -— +—-ansing +-VF-—- B+ — ¥ 
: v4 T L Tv 2Tr Pp? — B 
losed é FE F 
-_ ie — 
eCoTd ze =) + : +. Yo (4.3-17 ) 
oO F&F 2 1 2nP 
—s P? — B 
I, = 0 
i ap ¥ 
The second expression for Jae assumes P > B. When B = 0, these re- 
(} duce to the first terms of (4.2-5) and (4.2-7). By using a different 
1p is. |, method Middleton has obtained a more precise form of this result. 
Incidentally, for a given applied voltage, /a.(+-) for a positive bias | B 
is related to Zae(—) for a negative bias —  B | by 
- Ta(—) = |B\ + Tacl+) (4.3-18) 
Also r.m.s. T¢s(+) is equal to r.m.s. [y(/—). Equation (4.3-18) follows 
t 1 é | 
from a physical argument based on the areas underneath a curve of J for 
‘\ 
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the two cases. Both of the above relations follow from formulas given |) 


Middleton when V is the sum of a sine wave plus noise. They may also he 
derived from (4.3-15) and (4.3-16). 


4.4 Ovurpur POWER SPECTRUM 


The remainder of Part IV will be concerned with methods of solving the 
following problem: Given a non-linear device and an input voltage con- 
sisting of noise alone or of a signal plus noise. What is the power spectrum 
of the output? 

In some ways the answer to this problem gives us less information than 
the methods discussed in the first three sections. For example, beyond 
giving the rms value, it tells us very little about the probability density o/ 
the current corresponding to a given frequency band of the output. On 
the other hand, this rms value may be found (by integrating the power 
spectrum) for any band we choose to study. The methods described earlier 
depended on the input being confined to a relatively narrow band and gave 
information regarding only the entire band corresponding to a given har- 
monic (Oth, Ist, 2nd, etc.) of the input. There was no way to study the 
output when part of a band was eliminated by filters except by obtaining 
the power spectrum of some function of the envelope. 

At present there appear to be two general methods available for the 
determination of the output power spectrum each with its own advantages 
and disadvantages. First there is the direct method which has been used 
by W. R. Bennett*, F. C. Williams**, J. R. Ragazzini® and others. The 
noise is represented as the sum of a finite number of sinusoidal components 
The typical modulation product is computed and the output power spectrum 
is obtained by considering the density and amplitude of these products. 
The chief advantage of this method lies in its close relation to the known 
theory of modulation in non-linear circuits. Generally, the lower order 
modulation products are the only ones which contribute significantly to the 
output power and when they are known, the problem is well along towards 
solution. The main disadvantage is the labor of counting the modulation 
products falling in a given interval. However, Bennett has developed a 
method for doing this.” 

The fundamental idea of the second method is to obtain the correlation 
function for the output current. From this the output power spectrum may 
be obtained by Fourier’s transform. The correlation function method and 
its variations are of more recent origin than the direct method. They have 

* Cited in Section 4.0. Also much of this writer’s work on interference in broad band 
communication systems may be carried over to noise theory without any change in the 
methods used. 

** Cited in Section 4.1. 

%® Proc. I.R.E. Vol. 30, pp. 277-288 (June 1942), ‘The Effect of Fluctuation Voltages 


on the Linear Detector.” 
7 B.S.T.J., Vol. 19 (1940), pp. 587-610, Appendix B. 
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been discovered independently and at about the same time, by several 
workers. In a paper read before the I.R.E., Jan. 28, 1944, D. O. North 
described results obtained by using the correlation function. J. H. Van 
Vieck and D. Middleton have been using the two variations of the method 
which we shall describe in Sections 4.7 and 4.8, since early in 1943. A 
primitive form of the method of Section 4.8 had been used by A. D. Fowler 
and the writer in some unpublished material written in 1942. Recently, 
[ have learned that a method similar to the one used by Fowler and myself 
had already been used by Kurt Friinz in 1941." 

The correlation function method avoids the problem of counting the 
modulation products. However, in some cases it becomes rather unwieldy. 
Probably it is best to have both methods in mind when investigating any 
particular problem. The direct method will be illustrated by applying it 
to the square law detector. Two approaches to the correlation function 
method will then be described and applied to examples. 


4.5 Noise THROUGH SQUARE LAW-DEVICE 


Probably the most direct method of obtaining the power spectrum H(f) 
of J, where 
(4.1-1) 


V being a noise voltage, is to square the expression 


M 
V = Vy = > cn COS (wmt — ¢m) (2.8-6) 
1 


in which cm is 27( fm)Af, @m = 2tfmy fm = mAfand gi ,g2, +++ ¢ware random 
phase angles. 
Considerable simplification of the algebra results when we replace the 
representation (2.8-6) by 
l< . 
Ve= = Lae (4.5-1) 


- —2 


Here we have added a term co, 2 so as to not have any gaps in the summation 
and have introduced the definitions 


= Cy 
= —m (4.5-2) 
2rAf 


“Die Ubertragung von Rauschspannung iiber den linearen Gleichrichter,” Hocifr. 
u. Elektroakust., June 1941. Other articles by Franz are (I am indebted to Dr. North 
for the following references) “Beitrage zur Berechnung des Verhaltnisses von Signal 
spannung zu Rauschspannung am Ausgang von Empfangern”, E.N.7., 17, 215, 1940 and 
19, 285, 1942. “Die Amplituden von Gerauschspannungen”, E.V.T., 19, 166, 1942. 
The May 1944 (p. 237), issue of the Wireless Engineer contains an abstract of “The In- 
fluence of Carrier Waves on the Noise on the Far Side of Amplitude-Limiters and Linear 
Rectifiers” by Friinz and Vellat, E.V.T., Vol. 20, pp. 183-189 (Aug. 1943). 
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Squaring (4.5-1) gives the double series 
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Suppose we wish to consider the component of Vy of frequency fi, = £3) 


It is seen to be 


1 
A; cos (at — 7) om 2 Cran€, COS (Rat — Gren — Cn 
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The power spectrum IV(f) of J at frequency f; is a@ times the coefficient «/ 
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a - - % ptvpe. 
Afin the mean square value of (4.5-3) where the average is taken over t - 


y’s. Thus 


[=a 
2 +0 +0 
. \ a Bis 
VW (fi Af = < Ck—n Cn Ch—-m Cm ¥ 
tr —2 —o =. 
a: 
X ave. cos (kat — geen — ¢n) CoS (Rat — gem — 
where the summations extend over mand». Let » be fixed and consider 
those values of m which give an average different from zero. We see tha: The 
m = nand m = k — mare two such values. The only other possibilities JB of the 
arem = —nandm = —k + n, but these lead to terms containing (excep: JB produc 
when » or k equal zero) three different angles, gn, gin, ANd geyn Which Rand in 
average to zero. Using the fact that the average of cosine squared is one- Con 
half and that for a given # there are two such terms, we get to bet 
t=) f j 
onlv Vv 


2 +00 
W(f)as = > pi 


n=—0O 


+2 
=a Af > wf — fw ff 
n=—% 
where in the last step we have used 
fren = (Rk — WAf = fi — fn 
and have implied, from c_, = c,, that 
w( fin) = w(—ndf) = w(—fr) 
is equal to w(f,). 
Thus, from (4.5-4), we get for the power spectrum of / 
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with the understanding that f is not zero and 
w(—x) = w(x). (4.5-6 


[he result which is obtained by using (2.8-6), involving the cosines and 
only positive values of m, is 
x.) 


f 
W(f) = a | wx)w(f — x) dx + 2a° | wix)wif +a) dx (4.5-7 
; ! Jy j 


[his contains only positive values of frequency. (4.5-5) and (4.5-7) are 


equivalent and may readily be transformed into each other. 





The first integral in (4.5-7) arises from second order modulation products 
‘of the sum type and the second integral from products of the difference 
type. This may be seen by writing the current as 


‘ x ~~ 
I[=aV =a i z. CmCn COS (amt — Ym) COS (Wat — ¢n) 


m=1 n=1 


+ cos [ (wm i W, )t + Sm = ey |} 


The power in the range /; , f, + Afis the power due to modulation products 
of the difference type, w..f — wf, plus the power due to the modulation 
products of the sum type, w¢+ ef. In the first type fruns from 1 to = 
and in the second type f runs from 1 to k — 1. 

Consider the difference type first, and for the moment take both & and ¢ 

Fto be fixed. The twosets m= k+ (,n = fandm= (,n=k-+ fare the 

only values of m and # in (4.5-8) leading to a.¢ — we. The two corre- 
} sponding terms in (4.5-8) are equal because cos (—.x) is equal to cos x. The 
average power contributed by these two terms is 


9 


(; Creal «) xX {Average of (2 cos [(axit — wt)t — cert + etl)”} 


(4.5-9) 
= }(ackyt ct)” 


The power contributed to f; , f; + Af by the difference modulation products 
is obtained by summing f from 1 to «: 


~ S cipece = 2a” Dd) wl fertw(fo(af)? 
=l 


f 
+ 20? Af I w( fe + fw(f) af 


This leads to the second term in (4.5-7). 
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Now consider the modulation products of the sum type. The terms ¢' 
this type in (4.5-8) which give rise to the frequency a, are those for whi: 
m+ nisequaltok. Letn be 1lthenm =k — 1. The phase of this tern 
is random with respect to all the other terms except the one given by » = 
k — 1, m = 1 which has the same phase. The average power contribute 
by these two terms in (4.5-8) is, as in (4.5-9), 


\2 
5 (acycx-1) 
This disposes of two terms for which m + nis equal tok. Taking 2 to be) 
and going through the same process gives two more. Thus, assuming fo; 
the moment that & is an odd number, the power contributed to the intervs 
fi, fx + Af by the sum modulation products is 


1 &eR? ; 1 f Sie 
: 2 (ac, C.-»,) = : (QCnCe—n) a's | w(f)wif. — f) df 
4 n=1 4 n=] 0 


and this leads to the second term in (4.5-7). 
When the voltage V applied to the square law device is the sum of a nois: 
voltage Vy and a sine wave: 


V = P cos pt + Vy, (4.1-13 
we have 
1° = P* cos’ pt + 2PVy cos pt + Vi (4.5-10 


From the two equations 


cos’ pt = ; + 5 005 2pt 


ave. V4 


M 4 P 
Leas lf wna 


we see that J, or aV?, has a de component of 
= +a I w(f) df (4.5-11 


which agrees with (4.1-14), and a sinusoidal component 


P’ 
7 cos 2pt (4.5-12 


The continuous power spectrum W’.(f) of the remaining portion of 7 may 
be computed from 


2PVy cos pt + V;. 


Using 


2PVs 
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Using the representation (2.8-6) we see 


VW 
} 


2PVy cos pt = P>. CmlCOS (amt + pl — gm) + COS (wnt — pt — om) 


1 


For the moment, we take p = 2rrAf. The terms pertaining to frequency 
= ndAf are those for which 


Wm + p = 2rfn Wm Pp = 2afn 
mt+r=n m—-r\| =n 
m=n—r? m=rtn 


where only positive values of m are to be taken: If nm > r, then mis n — r 
orr+tu. lin <r, thenmisr—norr+n. In either case the values 
of mare na —r andu-+r. The terms of frequency f, in 2PVy cos pt 
are therefore 


Peia- COS (lafnt — Ojn—r})  Pente COS (Zafat — Cass 
and the mean square value of this expression, the average being taken over 
the y’s, is 


r 


2 


II 


(cinart + Cntr) = PAffw( fini) + w(fass)] 


P* Af (w(| Sn aig fo i) So w( fn > f,)| 


where f, denotes p/2r. 
By combining this with the expression (4.5-5) which arises from Vy 
we see that the continuous portion W.(f) of the power spectrum of J is 


Wf) =a Plw(f — f,) + wf +f,)] 
, (4.5-13) 
+a | w(x)w(f — x) dx 


where w(—/f) has the same value as w(f). 

Equation (4.5-13) has been used to compute W’.(f) as shown in Fig. &. 
The input noise is assumed to be uniform over a band of width 3 centered at 
fp, cf. Filter c, Appendix C. By noting the area under the low frequency 


portion of the spectrum we find 
8 . ” 
| Wf) df = a Bwo(P” + Bw) 
0 


Since the mean square value of the input I’y is Yo = 3wo, it is seen that 


this equation agrees with the expression (4.1-15) for the mean square value 
of It; , the low frequency current, excluding the d.c. If audio frequency 
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filters cut out part of the spectrum, HW’.(f) may be integrated over the re- 


maining portion to give the mean square value of the corresponding output 
current. This idea is mentioned in the footnote pertaining to equati: 
(4.1-6). 


If } consists of Vy plus two sinusoidal voltages of incommensurable fre- 
quencies, say 


Vv = Pcos pt + Ocos qt + Vy, 
CONTINUOUS PORTION OF OUTPUT SPECTRUM OF SQUARE LAW DEVICE 
INPUT = P COS 2TIf,t + NOISE 


OUTPUT OC = O(P2/2+ Bw) 


\ LET Bw@=c 


INPUT SPECTRUM 


K Set INPUT NOISE 
is 


2 ol A rn 


1 
' 
' 
\ w 1 
\ ! 
ae ie on ee” ee eee OR eee dt 


° 





fp 
° Be B ———— B- ot 2f,-B 
FREQUENCY 
Fig. 8 


the continuous portion W.(f) of the power spectrum of J may be shown to be 
(4.5-13) plus the additional terms 


vO fw(f — f.) + wf + fr] (4.5-14) 
where f, denotes g/ 27. 
When the voltage applied to the square law device (4.1-1) is® 


V(t) = O(1 + k cos pt) cos qt + Vy 
O O 
= (0 cos qt + a cos (p + gt + ” cos (p — g)t + Vy 
the resulting current contains the de component 
Qa 2 k = " - 
50 (1 + 4 +a | w(f) df (4.5-16) 
2 2 0 


*° A complete discussion of this problem is given by L. A. MacColl in a manuscript 
being prepared for publication. 
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The sinusoidal terms of J are obtained by squaring 
O11 + kos pt) cos gt 


and multiplying by a. The remaining portion of J has a continuous power 
spectrum given by 


Wf) =a Ol wf —f.) + wif +f) 
i 
+—w(f—f,-f)+—wftft+f) 
gl Ie J g@U th, — 
P o (4.5-17) 


k ! 4 ' ke oe f 
a 4 wf — fp T Je) i 4 wey T Ip aa fq) 


+x 
+a / w(a)w(f — x) dx 
— 


where f, denotes p/2z and f, denotes q/ 27. 


4.6 Two CorrELATION FUNCTION METHODS 


As mentioned in Section 4.4 these methods for determining the output 
power spectrum are based on finding the correlation function V(r) for the 
output current. From this the power spectrum, H(/), of the output cur- 
rent may be obtained from (2.1-5), rewritten as 

« 
Wf) = s| W(r) cos 2xfr dr (4.6-1) 
0 
It will be recalled that H(f)Af may be regarded as the average power which 
would be dissipated by those components of J in the band /, f + Af if J were 
to flow through a resistance of one ohm. 

The input of the non-linear device is taken to be a voltage Vit). It may, 
for example, consist of a noise voltage Vy(¢) plus sinusoidal components. 
The output is taken to be a current /(¢).. The non-linear device is specified 
by a relation between V(¢) and /(¢). In this work J(¢) at time ¢ is assumed 
to be completely determined by the value of V(¢) at time ¢. 

Two methods of obtaining V(r) will be described. 

(a) Integrating the two-dimensional probability density of Vi‘) and 
V(t + 7) over the values allowed by the non-linear device. This 
method, which is especially direct when applied to noise alone through 
rectifiers, was discovered independently by Van Vleck and North. 

(b) Introducing and using the characteristic function, which for the sake 
of brevity will be abbreviated to ch. f., of the two-dimensiona! prob- 
ability distribution of V(t) and V(t + 7). 
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4.7 Linear DetTecTION OF NoIsE—THE VAN VLECK-NoRTH METHOD 


The method due to Van Vleck and North will be illustrated by using jt 
to determine the output power spectrum of a linear detector when the input 
consists of noise alone. 

The linear detector is specified by 

{0, Vit) <0 


T(t) = LV), Vit) > 0, (4.7 | 


which may be obtained from (4.2-1) by setting a equal to one, and the input 
voltage is 


V(t) = Vy(t) (4.7-2 


where l'y(f) is a noise voltage whose correlation function is ¥(+) and whose 
power spectrum is w(f). 

The correlation function V(7) is the average value of 7(#)/(¢ + 7). This 
is the same as the average value of the function 


F(V,, V;) = {Vivo when wigobdl Ve>0 (47-3 
\0, all other V’s, 


where we have set 
Vi = Vit) 
Ve= Viet 7) 


The two-dimensional distribution of 1) and V9 is given by (3.2-4), and 


from this it follows that the average value of any function F(l,, V2) is 
ewer hae, FP ] 1 72 7 ess 
L d\ 1 [. di 2 on M iy ep[ 3 iy (Wo | 1 +- oI _ 2y, I v2 | 
(4.7-4) 
where 
IM|=y-¥;. 


For the linear rectifier case, where F(V1, Ve) is given by (4.7-3), the 
integral is 


| |"? a av, [ dV2V,Veexp| — WoVi + woVE — MeVaV9 | 
F Qi Jo 0 2\|M| 
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where we have used (3.5-4) to evaluate the integral. The arc cosine is 
taken to be between 0 and 7. We therefore have for the correlation func- 


tion of I(t), 
V(r) = [yi = vil "+ Y, COS 1] —¥s $.7-5 
; dr Tae Yo 


The power spectrum IV(f) may be obtained from this by use of (4.6-1 
For this purpose it is convenient to write (4.7-5) in terms of a hypergeo- 
metric function. By expanding and comparing terms it is seen that 


V(r) = e a vo F(-3, $s 13%) 


2m r Wo 
(+.7-6) 
ys 4 ‘ 
= vs + vi + ¥ + terms involving y;, ¥;, etc. 
+ 2r try 


' 


As will be discussed more fully in Section 4.8, a constant term 7 in Yr 
indicates a direct current component of 7(f) of .f amperes. Thus /(¢) has 
a de component equal to 


1/2 1 
[2] = —>= X rms value of V(t) (4.7-7 
Vv 


2r dar 


This agrees with (4.2-3) when the P of that equation is set equal to zero, 
Integrals of the form 
2 
Gilt) = | Y; cos 2xfr dr 

0 
which result when (4.7-6) is put in (4.6-1) and integrated termwise are 
discussed in Appendix 4C. From the results given there it is seen that if 
we neglect ¥} and higher powers we obtain an approximation for the con- 
tinuous portion W’.(f) of W(f): 


Wi) = Gif) + 2 
To a = 9 
_ w(f) Sl ere ee 
= 4 — inte 5 i w(x)w(f x) dx 


where w(—f) is defined as w/(f). 

When V(t) is uniform over a relatively narrow band extending from 
fa to fyso that w(f) is equal to wo in this band and is zero outside it, we may 
use the results for Filter c of Appendix 4C. The fo and 3 given there are 
related to fa and f, by 

; fy = fo + 


Ta = fo at ’ 2? 
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and the value of wo taken there is the same as here and is ¥o/ 8. The value 
of Golf) given there leads to the approximation, for low frequencies: 
en 1 Wo f 
Wf) = Yo [= 
TY 43 p a 
(4.7-9) 


coe (1 ee ) 
4dr to a ta 


when 0 < f < fp — fa, and to Wo(f) = 0 for fp — fa << f < fa. By setting 
P equal to zero in the curve given in Fig. 8 for W.(f) cerresponding to the 
square law detector, we see that the low frequency portion of the power 
spectrum is triangular in shape and is zero at f = 3. Thus, looking at 
(4.7-9), we see that to a first approximation the shape of the output power 
spectrum is the same for a linear detector as for a square law detector when 
the input consists of a relatively narrow band of noise. 

An approximate rms value of the low frequency output current may be 
obtained by integrating (4.7-9) 


— fifa 
i, l We(f) df 
_ wol fo — fa) _ Yo 


Sit Sir 


> 


; 1 
rms low freq. current = i X rms applied voltage (4.7-10) 
V or 


It is seen that this is half of the direct current. It must be kept in mind 
that (4.7-10) is an approximation because we have neglected ¥; and higher 
powers. The true value may be obtained from (4.2-8). It is seen that the 
coefficient (87)~"”” = 0.200 should be replaced by 


1 : r 1/2 
-[2 — = 0.209 
T 2 


W.(f) for other types of band pass filters may be obtained by using the 
corresponding G’s given in appendix 4C. It turns out that (4.7-10) holds 
for all three types of filters. This is a special case of Middleton’s theorem, 
mentioned several times before, that the total power in any modulation 
product (it will be shown later in Section 4.9 that the termy; in (4.7-6) 
corresponds to the n“™ order modulation products) depends only on the 
total input power of the applied noise, not on its spectral distribution, 


4.8 Tre CHARACTERISTIC FuNcTION METHOD 


As mentioned in the preceding parts, especially in connection with equa- 
tion (1.4-3), the ch. f. of a random variable x is the average value of exp 
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(jux). This is a function of wu. The ch. f. of two random variables x and 
vis the average value of exp (tux + ivy) and is a function of wand. The 
ch. f. which we shall use here is the ch. f. of the two random variables V(t 
and V(t + 7) where V(/) is the voltage applied to the non-linear device, and 
the randomness is introduced by ¢ being selected at random, 7 remaining 


fixed. We may write this characteristic function as 


we yy i gatas ra 
g(u, v, 7) = Limit — [ exp [iuV(t) + wVit + 7)) dt (4.8-1 
rx Jo 
If V(t) contains a noise voltage Vy(t), as it always does in this section, and 
if we use the representation (2.8 -1) or (2.8-6) a large number of random 
parameters (d,’s and 6,’s or ¢,’s) will appear in (4.81). In accordance 
with our use of such representations we may average over these parameters 
without changing the value of (4.8-1) and may thereby simplify the integra- 
tion. 
For example suppose 
Vié) = Vt) + Val?) (4.8~2) 


where I’,.(¢) is some regular voltage which may, e.g., consist of one or more 
sine waves. Substituting this in (4.8-1) and using the result (3.2-7) that 


the ch. f. of Vy(t) and Vy(¢ + 7) is 


(uw, 2, 7) = ave. exp [iuVy(t) + iV y(t + 7 


, 
aN 


aw 


ae oe . (4.8-: 
= exp ms. (iu +-v) — yuo 


' 


¥, = V(r) being the correlation function of Vy(¢), we obtain for the ch. f. 


of V(t) and Vit + 7), 
g(u, vt, rT) = exp| —~ (uw +0) — pur 
A ee a ee . (4.8-4) 
< Limit = exp [iuV.(t) + ioV.(¢ + 7)] di 

lx 7 “0 


= gv(u, v, 7) feu, 2, 7) 


In the last line we have used g,(u, v, 7) to denote the limit in the line above: 


WY ges ae, 
g.(u, v, 7) = Limit | exp [iat s(t) + ivV,(¢ + 7)| dt (4.8-5 
0 


Tx 


The principal reason we use the ch. f. is because quite a few non-linear 


devices may be described by the integral 


P= -f Fine *" du (44-1) 
Pas c ‘ 
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where the function F(iw) and the path of integration C are chosen to fit the 
device. Examples of such devices are given in Appendix 4A. The corre- 


lation function W(r) of Z(t) is given by 


II 


7 
V(r) Limit | (I(t + 7) dt 
T Jo 


Tx 


io 


Limit wT | dt [ F(iuye' vit au | Five! Vit+r) dv 
I 7 /0 Jo ” 


II 


= [ F(iu) du i F(iv) dv (4.86 


. 
Limit o | exp [ial (t) + ivV(t + 7)] dt 
0 


T-2 


= = [ F(iu) du , F (iv) g(u, v, 7) dv 


This is the fundamental formula of the ch. f. method. 
When I(t) is the sum of a noise voltage and a regular voltage, as in 
(4.8-2), (4.8-6) becomes 


V(r) = a F(iu)e wont dy [ F(iv)e 0! *” 
wip . (4.8-7) 
g(u, v, T) dv 


oe ¥ru? 


where g.(u, 7, 7) is the ch. f. of V(t) and V.(¢ + 7) given by (4.8-5). This 
is a definite expression for V(r). All that follows is devoted to the evalua- 
tion of this integral and to the evaluation of 


Wf) = +/ V(r) cos 2nfr dr (4.6-1) 


for the power spectrum of J. 

(Juite often /(¢) will contain dc and periodic components. It seems con- 
venient to deal with these separately since they correspond to terms in 
V(7) which cause the integral (4.6-1) for W(f) to diverge. In fact, from 
Section 2.2 it follows that a correlation function of the form 


ee 
A’ + = cos 2rfot (2.2-3) 
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corresponds to a current 


A + C cos (2xf¢ — ¢) (2. 
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where the phase angle ¢ cannot be determined from (2.2-3) since it does not 
affect the average power. 

Consider the correlation function for V(t) = Vit) + Vy(t) given by 
(4.8-2). Itis 


1 T T 
Limit =, | V2)V 606 + 7) dt + [ VA2Va(t + 7) dt 
T—2 7 0 “0 

(4.8-8) 


r 7 
+ | Vu(QVCt + 7) dt + Vue(t)Va(t + 7 at| 
0 0 


Since V.(¢) and Vy(¢t) are unrelated the contributions of the second and 
third integrals vanish leaving us with the result 


Correlation function of V(t) = Correlation function of V7. (¢ 
+ Correlation function of ly(f). 


Now as t — & the correlation function of Vy(¢) becomes zero while that of 
V(t) becomes of the type (2.2-3) given above. Hence the correlation func- 
tion of the regular voltage V(t) may be obtained from V(t) by letting 7 = 
and picking out the non-vanishing terms. Although we have been speaking 
of V(t), the same results hold for /(¢) and this process may be used to pick 
out those parts of W(7) which correspond to the de and periodic components 
of Z(t). Thus, if we look at (4.8-7) we see that as 7 «,Y,— 0, while the 
g. (uw, v, 7) Corresponding to V.(t) given by (4.8-5) remains unchanged in 
general magnitude. This last statement may be hard to see, but examina- 
tion of the cases discussed later show that it is true, at least for these cases. 
Thus the portion of W(7) corresponding to the de and periodic components 
of I(t) is, setting y, = Oin (4.8-7), 


> 


V,(7r) = is 


JC 


F(iu)e%°° du [ F(iv)e “*""'''g.(u, v, 7) dv (4.8-10 


where the subscript * indicates that V.(7) is that part of ¥(r) which does 
not vanish as tT —> ~. 
We may write (4.8-9), when applied to 7(¢), as 


We) = W,.(7) + W.(r) (4.8-11 


where V,(7) is the correlation function of the “continuous” portion of the 
power spectrum of /(¢). 

Incidentally, the separation of W(7) into the two parts shown in (4.8-11 
may be avoided if one is willing to use the 6(/) functions in order to interpret 
the integral in (4.61) as explained in Section 2.2. This method gives the 
proper de and sinusoidal components even though (4.6-1) does not con- 
verge (because of the presence of the terms leading to V,(7 
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4.9) Noise Pius SINE Wave ApplLieD TO NoNn-LINEAR DEVICE 
In order to illustrate the characteristic function method described in 
Section 4.8 we shall consider the case of a non-linear Cevice specified by 


| LY ieee a here 
[ = | F(iuje’’“ du (4A-1 


2a J, 
when VT’ consists of a noise voltage plus a sine wave: 
V(t) = P cos pt + Vy(t) (4.1-13 


As usual, Vy(t) has the power spectrum w(f) and the correlation functior 
Yi7). (7) is often written as ¥, for the sake of shortness. Comparing 
(4.1-13) with (4.8—2) gives 


V.() = P cos pl (4.9-1 
Our first task is to compute the ch. f. g,(u, 2, 7) for the pair of random 
variables 1’.(¢) and V.(¢ + 7). We do this by using the integral (4.8-5): 


el 


al ; _ 
Limit ;, | exp [iuP cos pt + ivP cos p(t + 7)] dt 
T’ Jo 


T 0 


oT 
- 
2 
9 
~ 
| 


(4.9-2 
= Jo PV w + v0 + luv cos pr) 
where Jy is a Bessel function. The integration is performed by writing 


ucos pl + vcos p(t + 7) = (uw + vos pr) cos pt — vsin pr sin pl 





= Vu + v + 2uv cos pr cos (pt + phase angle 
and using the integral 
1 7 
Jo{a) = " ‘eine 
2m Jo 
The correlation function for (4.1-13) has also been given in Section 3.10. 
The correlation function (7) for Z(t) may now be obtained by substi- 
tuting the above expressions in (4.8-7) 


1 % e Lnl?)u2 rr (Wo/2)v2 
V(r) = : | du Fiuj)e “or [ dv F(iv)e °°" 
4° ( vc 


(4.9-3) 
er To PV18 + e+ Quy cos pr). 


W.(7), the correlation function for the d.c. and periodic components of J, 
may, according to (4.8-10), be obtained from this by setting ¥, equal to zero. 

When we have a particular non-linear device in mind the appropriate 
F(in) may often be obtained from Appendix 44. For example, (im) for a 
linear rectifier is —« ~. Inserting this value in (4.9-3) gives a definite 
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double integral for W(7). If there were some easy way to evaluate this in- 
tegral then everything would be fine. Unfortunately, no simple method of 
evaluation has vet been found. However, one method is available which is 
closely related to the direct method used by Bennett. It is based on the 


expansion 


ga(u,v, 7) = Jo PV 12 + v + 2uv cos pr) 


7” 
= 2 €,(—)"Jn(Pu)Jn(Pr) cos npr (4.9-4 
7 0 
a= 1, g =2 for #2 1 


This expansion enables us to write the troublesome terms in (4.9-3) as 


ef Jo Pv 2 + 2 + 2Quv cos pr) 


al 2 k (4.9-5 
Lk (W, uv) coal 
= Zz a (—)"""e,, cos npr as JA Pat, (Pe 
nee &k 0 A. 


The virtue of this double sum is that it simplifies the integration. Thus, 
putting it in (4.9-3) and setting 


h, t= 


[ Fliu)ué J,,(Puye ‘9? ” du (4.9-6 


ci 


inte 
2a J 
gives 
7” a] 1 ‘ 
V(r) = z pe Rh! wih’ k€, COS NPT (4,9-—7) 


The correlation function V,(7) for the de and periodic components of J 
are obtained by letting 7 — « where ¥, — 0. Only the terms for which 
k = 0 remain: 

n” 


W(t) = D> enh no cos npr (4.9-8) 


n=() 


Comparing this with the known fact that the correlation function of 


A + C cos (2mft — ¢) (2.2-2 


Pn > 
Av + 5 cos 2rfor (2.2-3) 


and remembering that € is one while €, is two for nm > 1 shows that 


Amplitude of dc component of J = hoo 


4 / 
(4,9-9) 


, nn , 
Amplitude of “: component of T= 2hpo 
T 
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Incidentally, these expressions for the amplitudes follow almost at once fro; 
the direct method of solution. This will be shown in connection with equa- 
tion (4.9-17). 

Since the correlation function V,(7) for the continuous portion W,(f) of 
the power spectrum for J is given by 


V(r) = V(r) — V.(7), (4.8-11 
we also have 
V(r) = 5 ie Z — WW. €n COS npr (4.9-10 
n=) kel Rt 


When this is substituted in 


Wf) = s/ W.(r) cos 2afr dr (4.9-11 
0 


we obtain 
x 2) ) x 
Wij) eae a eee (4.9-12 
aa eet Qtr dr 


G.(f) = [ Vv cos 2nfr dr (4.9-13) 
0 


is the function studied in Appendix 4C. G,(f) isan even function of f. The 
double series (4.9-12) for HW’. looks rather formidable. However, when we 
are interested in a particular portion of the frequency spectrum often only 
a few terms of the series are needed. 

It has been mentioned above that the direct method of obtaining the out- 
put power spectrum is closely related to the equations just derived. We 
now study this relation. 

We start with the following result from modulation theory”: Let the 
voltage 

e Vo= Poycos xo + P; cos xy + +++ + Pw cos xy 
(4.9-14) 
x = pt, k=0,1,°°+N, 
where the p;’s are incommensurable, be applied to the device (4A-1). The 
output current is 


Pa x 
Tk: ee: a ; 
mo=V my=0 — 


\ 
* Emy COS MoXo COS 1X, *** COS My XN 


5° Bennett and Rice, “Note on Methods of Computing Modulation Products,” Phil. 
Mag. 8.7, V. 18, pp. 422-424, Sept. 1934, and Bennett’s paper cited in Section 4.0. 
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where €9 = land, = 2 for m > 1. When the product of the cosines is 
expressed as a sum of cosines of the angles mo xo + my x1 +++ Emyxty , it is 


seen that the coefficient of the typical term is Am .my , except when all 


the m’s are zero in which case it is 3-lo..0. Thus 


>Aoo...0 = de component of J 


Amo---my | = amplitude of component of frequency  (4.9-16 


a | Mo Po mpi H +++ + My py | 


For all values of the m’s, 
Nee hd \ 
Ans-.-my = | F(iu) II Jm,(P,u) du - 
us Cc r= (4.9-17) 
M =m +m +--+ + my 

Following Bennett’s procedure, we identify Vas given by (4.9-14), with 
V = P cos pt+ Vw (4.1-13) 
by setting Po = P, po = p, and representing the noise voltage Vy by the sum 
of the remaining terms. Since this makes P;, Py all very small, Laplace s 

process indicates that in (4.9-17) we may put 


II J0(P,- 0) = exp — nei (Pi teen Hh Px) 
oa ? (4.9-18) 


—yoqu*/2 


= €é 


We have used the fact that Yo is the mean square value of Vy. It follows 


from these equations that 


de component of I 


1 ss J 2)u2 
s— | F(t) Jol Puje "du 
27 Je 


II 


Component of frequency - i F(iu)J,(Puye °°"? du 
2x «de 
These results are identical with those of (4.9-9). 

The equations just derived show that /, is to be associated with the n" 
harmonic of p. In much the same way it may be shown that /,, is to be 
associated with the modulation products arising from the n‘” harmonic of 
p and k of the elementary sinusoidal components representing Ty. We 
consider only combinations of the form p; + po + p;, taking k = 3 for ex- 
ample, and neglect terms of the form 3p; and 2p; + pe. The former type 
is much more niMerous, there being about V" of them while there are only 
about V and \’, respectively, of the latter type. 
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We again take k = 3 and consider m, , my, m3 to be one, and m4, +++) 
to be zero, corresponding to the modulation product np + p; + pe + { 
By making the same sort of approximations as Bennett does we find 

n+3 
l P, Py» P3 s/s 3 _(—wu2/2) 
A n,i,1.1,0;00000 = F(iu\Tn(Pujui ee“??? du 
" 


T & 


P, PsP 
= — > es 
4 
When any other modulation product of the form np + p,, + pr, + pr, is 
considered we get a similar expression in which P;P2P3 is replaced by 
P,,P,,P,,. This may be done for any value of k&. The result indicates 
that /,,, and consequently also the (x, k)" terms in the double series 
(4.9-10) and (4.9-12) for V.(7) and W-(f), are to be associated with the 
modulation products of order (7, k), the referring to the signal and the & 
to the noise components. 

We now may state a theorem due to Middleton regarding the total power 
in the modulation products of a given order. For a given non-linear device 
(i.e. F(iu) is given), the total power which would be dissipated by all of the 
modulation products which are of order (n, k) if J were to flow through a 


resistance of one ohm is 


e.ly(0)|* 
k! 


772 1k 22 
2 éAV xh, 
Vn(O) = hn. = 3 = (4.9-19) 
The important feature of this expression is that it depends only on the r.m.s. 
value of Vy and on F(in). It depends not at all upon the spectral dis- 
tribution of the noise power in the input. 
Che proof of (4.9-19) is based en the relation 


Vai(0) -[ Wri f) af 


between the total power dissipated by all the (7, &) order products and the 
corresponding correlation function obtained from (4.9-7). 

This theorem has been used by Middleton to show that when the input 
is confined to a relatively narrow frequency band, so that the output spec- 
trum consists of bands, the power in each band depends only on Vy and not 
on the spectrum of Vy. 


4.10 MusceLLANEOUS RESULTS OBTAINED BY CORRELATION FUNCTION 
METHOD 


In this section a number of results which may be obtained from the theory 
given in the sections following 4.6 are given. 
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When the input to the square law device 


T= al” (4.1-1) 
consists of noise only, so that 1 = Vy , the correlation function for 7 is 


W(r) =a lye + 27) (4.10-1) 


where y,; is the correlation function of Vy. This may be compared with 


equation (3.9-7). When V is general, 
W(r) = ave. I[(t)I(t + 7) 


ave.a’ V(t)V"(t + 7) 
iu)” (iv) (4.10-2 


2 "i — - (83 v ° . r 
a’ X Coefficient of ,-~,,- in power series expansion 





of ch. f. of V(t), V(t + 7) 


where we have used a known property of the characteristic function. 
expression for the ch. f., denoted by g(w, 7, 7), is given by (4.8-4). I 
example, when IV consists of a sine wave plus noise, (4.1-13), the ch. f. is 
obtainable from (4.9-3). Hence, 


» 9 


ca eee ae 
V(r) = Coeff. of ; in expansion of 


a Io PV? + + uv cos pr) 


a ° (4.10-3) 
x exp | -¥ (uo a scatinds 


- 


iP , Poe ‘i “ae 
= a ~ vo} + - cos 2pr + 2Py, cos pr + 2y, 


« 


The first two terms give the de and second harmonic. ‘The last two terms 


may be used to compute HW’,(f) as given by (4.5-13). 

Expressions (4.10-1) and (4.10-3) are special cases of results obtained by 
Middleton who has studied the general theory of the quadratic rectifier by 
using the Van Vleck-North method, described in Section 4.7. 
theory of Section 4.9 may be applied we con- 


As an example to which the 
(4.1-13), to be applied to the v-law rectifier 


sider the sine wave plus noise, 
T r<0 
; 4.10-4) 
Ful, V>0 
From the table in Appendix 44 it is seen that 


F(iu) = Tv + 1)(11 
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and that the path of integration C runs along the real axis from — * to « 
with a downward indentation at the origin. The integral (4.9-6) for / 
becomes 
n+k—v—1 
7 k 1 f \ Yo/2)u2 
hy = 3 (vy + | e” J,(Puje *” du 
Tv c 


(@)" . n/2 yy 

x 'T(v + 1) 

) ej int 
= - ——— : Ff =. 5 
(eee 
21 aes = —— J: 

P’ 

2o 


where the integration has been performed by expanding J,(Pu) in powers 


ai n+ 1; -x) (4.105 


of u and using 
—au? | 2-1 + —hi AU: . 
| eo a du = ie "a ” sin Ax¥I'(A) 
p 


A 
= se ai—-—<-™ ro) (4.10-6 


° —hir 
ime 


at ari — ) 


it being understood that arg « = 0 on the positive portion of C, 
From (4.9-9), the de component of J is 


T(1 a v) »/2 vy 
hig =~ . iP; aos 5 1: —@ (4.10 -7 
v é rs 
2r (1 + *) 
which reduces to the expression (4.2-3) when v = 1 for the linear rectifier 


(aside from the factor a). 

When the input (sine wave plus noise) is confined to a relatively narrow 
band, and when we are interested in the low frequency output, consideration 
of the modulation products suggests that we consider the difference products 
from the products of order (0, 0), (0, 2), (O, 4), --- (1, 1), (1, 3), --+ (2, 0), 
(2, 2), -++ etc. where the typical product is of order (7, k). The orders 
(O, 0) and (2, 0) give the dc and second harmonic and hence are not con- 
sidered in the computation of W.(f). Of the remaining terms, either (0, 2) 
or (1, 1) gives the greatest contribution to the series (4.9-12) and (4.9-10) 
for W.(f)andW,(7). The remaining terms contribute less and less as 7 and 
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k increase. The low frequency portion of the continuous portion of the 
output power spectrum is then, from (4.9-12), 


Os Pe 
Wf) = 5, Mo2Galf) + gy mG) + 


+ 


4 his(G3(f = fo) (4.10-8 


3! 


ae = ar 
“bt jy Mla — fo) + Gif + fo)| + 


+GAi + fol + SRG — BO + Cd + BOI 


From Table 2 of Appendix 4C we may pick out the low frequency portions of 
the G’s. It must be remembered that G,,(.°) is an even function of x and 
thatO < f<K fo. 

As an example we take the input noise Vy to have the same w(f) and 
¥(7) as Filter a, the normal law filter, of Appendix 4C, so that 


f—fo)2/202 
w(f) = elles 
: 
and assume that the sine wave signal is at the middle of the band, giving 


p = 2nfo. Thus, from (4.10-8), for low frequencies and the normal law 
distribution of the input noise power, 


Wf) = | hispie 2/402 wn 1 Iie” 8¢2 
: tov x O40 lr 
+r an/ dg Nit hoe ial juntas Iispoe’? '” = (4.109) 
+ el 5 
doy T 


Although we have been speaking of the v-law rectifier, equation (4.10-9 
gives the low frequency portion of IV.(f), corresponding to a normal law 
noise power, for any non-linear device provided the proper /,,’s are inserted. 

When we set v equal to one in the expression (4.10-5) for /,; we may ob- 
tain the results given by Bennett. Middleton has studied the output of a 
biased linear rectifier, when the input consists of a sine wave plus noise, and 
also the special case of the unbiased linear rectifier. He has computed the 
output for a wide range of the ratios P*/Yo , B’/ Yo where Bis the bias. In 
order to cover the entire range he had to derive two series for the corre- 
sponding /,x’s, each series being suitable for its particular portion of the 
range. 
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A special case of (4.10-9) occurs when noise alone is applied to a lincar 


rectifier. The low frequency portion of the output power spectrum is 


eo ) = (—1} \im( — >) 1 2'4me2 
" zr fy = vi » Bi \ 2 : 2 e f*/4 
Z T m=1 mllm. oV dimer 
_ vor [grt yb rise (4.10-10 
= 26 UC | 644/2 é : ( 


1 ms 
‘aa 120° ee 
T 9564/3 ¢ ss ] 


we have used (4.7-6) and Table 2 of Appendix 4C. 
correlation 


where 
The function of 


V. = P cos pt + QO cos qt, 


where p and g are incommensurable, is 


Jo PV w + 0? + luv cos pr) K J(Qv 1? + v + 2uv cos qr) 
From equations (4+.9-16) and (4.9-17) it is seen immediately that 


hoo = : | Fu)Jo( Pu)To(Quwye “'? du (4.10-11 
27 JC 


is the d.c. component of J when the applied voltage is 
P cos pt + Q cos gt + Vy. (4.1-4 


J. R. Ragazzini has obtained an approximate expression for the output 
power spectrum when the voltage 


V, + Vy 
O(1 + r cos pt)cos gt 


= 
V; 


(4.10-12) 


HI 


ss i ‘ eas 46 > ‘ ‘ s 
is impressed on a linear rectifier.”. In terms of our notation his expression 


for the continuous portion of the power spectrum is (for low frequencies ) 


WA) = 1 (J) given by equation 


nr a? (QO? + 20) (4.5-17) for square law device 
The a@ is put in the denominator to cancel the a in the expression (4.5~17). 
We take the linear rectifier to be 


| (4.10-13 


(4.10-14 


and replace the index of modulation, &, in (4.5-17) by r. 


* Equation (12), “The Effect of Fluctuation Voltages on the Linear Detector,” Pr 
I.R.E., V. 30, pp. 277-288 (June 1942). 
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A 
near  -Ragazzini’s formula is quite accurate when the index of modulation r is 
i small, especially when y = O° (2Yo) is large. To show this we put r = 0 
in (4.10-13) and obtain 
W.(f) = | oe f hy Or wee: + f) 
_ wm? + 2%) L* ee ngg) od ni indy 


: (4.10-15) 
~1{) 4 rs i 
: a i wix)w(f — x) dx 
‘ 


where fg = q/(27m). 
tion of W.(f) obtained by specializing (4.10-8) to obtain the output power 


This is to be compared with the low frequency por- 


spectrum of a linear rectifier when the input consists of a sine wave plus 
The leading terms in (4.10-8) give 


noise. 


Bow.) = Rilw(fe —f) + wf, +f] 


,. rt? (4.10-16) 
+ hi» i / w(x)w(f — x) dx 


The values of the /’s appropriate to a linear rectifier are obtained by set- 
ting y = 1 in (4.10-5) and noticing that Q now plays the role of P. 


11 1/y\!2 
hy = : ( ) 1Fi(33 2; —y) 
2\r 


hon = (24o)? Fi(3; 1; —y) 
y = Q°/(2y0) 


Incidentally, the first approximation to the output of a linear rectifier 


(4.10-17) 





given by (4.10-16) is interesting in its own right. Fig. 9 shows the low fre- 
quency portion of IV’.(f) as computed from (4.10-16) when the input noise 
is uniformly distributed over a narrow frequency band of width 3, f, being 
the mid-band frequency. /y and /ig may be obtained from the curves 
shown in Fig. 10. In these figures P and x replace Q and v of (4.10-17) in 
order to keep the notation the same as in Fig. 8 for the square law device. 
These curves may also be obtained from equations (33) to (43) of Bennett’s 
paper. 
The following values are useful for our comparison. 


When x = 0 When x is large 





hy 





hoo 


The values for large x are 
Pro. J) given in Appendix 4B. 








= 0 


= (2m) 


Ityy 


hos 


1/x 


1/(xQ). 


(4.10-18) 


obtained from the asymptotic expansion (4B — 3) 
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LOW FREQUENCY OUTPUT OF LINEAR RECTIFIER 
APPROXIMATION ~SECOND ORDER PRODUCTS ONLY 


INPUT= V=FPcos anfpt + NOISE 


450, V0 
ourput=1= | % he 


OUTPUT D.C.= Ph,,+ B We hos 
2 Bw? i 2 (pws 
LET C=hy, —Z> = (Phog) om 9 
INPUT SPECTRUM 


INPUT NOISE 
oa ss 





-e—- 


' 
' 
' 
' 
' 
' 
rn 








fp 


“aia, 














FREQUENCY 
Fig. 9 
0.4,-— RE an en ee: ED = ania 
0.3 
0.2 
“0 0.5 ° 5 ~ 2.0 “2.5 30 aaa 7 ‘ 
x = AVE SINE WAVE POWER pe 
AVE NOISE POWER 2 BW 


Fig. 10-—Coetficients for linear detector output shown on Fig. 9 
is , , 
, 2 sa 
Phy, = V iF, (4; 1; —x) i = \ iF\(4; 2; —x) 
y= 


We make the first comparison between (4.10-15) and (4.10-16) by letting 
Q— x, It is seen that both reduce to 


WA) = Slo —) + 0h +I (4.10-19 
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which shows that the agreement is perfect in this case. Next we let O = 0. 
[he two expressions then give 


Wf) = aia, lf w(x)w(f — x) dx 


H 
| 
; 


‘where A = @ for Ragazzini’s formula and A = 4 for (4.10-16). Thus the 
fagreement is still quite good. The limiting value for (4.10-16) may also 
‘be obtained from (4.7-8). 


; 


' Even if the index of modulation r is not negligibly small it may be shown 
that when Q — x W(f) still approaches the value given by (4.10-19). 
Ragazzini’s formula gives a somewhat larger answer because it includes the 

additional terms, shown in (4.5-17), which contain &°/4, but this difference 
‘ does not appear to be serious. If the Q* + 2¥oin the denominator of (4.10 
13) be replaced by QO” + 10°k” + yo the agreement is improved. 


APPENDIX 4A 
TABLE OF NON-LINEAR DEVICES SPECIFIED BY INTEGRALS 
(Quite a number of non-linear devices may be specified by integrals of the 


form 


I= i. F(iu)e'’" du (4A-1) 
2r c 





- where the function /(iu) and the path of integration C are chosen to fit the 


device.* The table gives examples of such devices. Some important cases 


' cannot be simply represented in this form. An example is the limiter 





Eon erat cians 











I = — aD, V<—-D 

[= al, —-D<V<D 

I = aD, D<V (4A-2) 

’ which may be represented as 
ee gad. lu 
I = 7 sin Vu sin Du — 
7 Jo uw? 

. , (4A-3) 
2 ar u 

= —aD + at} e'" sin Du as 
2m Je uw 

where C runs from — * to + * and is indented downward at the origin. 


This is not of the form assumed in the theory of Part IV. However it 
appears that it would not be difficult to extend the theory in the particular 
case of the limiter. 


* Reference 50 cited in Section 4.9. 
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Non-Linear Devices Specified by Integrals 


[ F(iude’™™ du 
wJe 
C 


I = aV", n integer | Positive Loop | mth power devi 


around u = 0 | 
| 


T=a(V "Hm : | Positive Loop | nth power devi 
integer inva ° |} around#w =O | with bias 
| 


I = 0, Pie | Real « axis from | Linear rectifier 
I =alV, <a — | —«to+owith) cut-off at 
| downward _ in- V=0 
|  dentation at | 
u=0 


I = 0, V<B ie vth power recti 
I = a(V — B)’, iu)’ fier with bias 
¥ >-o | 

vy any positive number 


= 0, V <0 = | Linear __rectitier 
=aV, O<V<D ok | plus limiter 
= aD, D * Vv | 





Yy <0 
V>0 


APPENDIX 4B 


THe FuNctTION 1F;(a; ¢; x) 
In problems concerning a sine wave plus noise the hypergeometric func- 


tion 


az 4 ala+1)2 
cl! ' cle + 1) 2! 
arises. Here we state some of its properties which are of use in the theory 
of Part IV. Curves of (a; ¢; 2) are given for a = — 4, — 3.5-++ , 3.5, 
4.0 andc = — 1.5, — .5, + .5, 1, 1.5, 2, 3, 4 in the 1938 edition, page 275, 
of “Tables of Functions”, by Jahnke and Emde. A list of properties of the 
function and other references are also given. In addition to these refer- 
ences we mention E. T. Copson, “Functions of a Complex Variable” (Ox- 
ford, 1935), page 260. 
If cis not a negative integer or zero 


iF (a; C3 s) = 1 oe a was (4B-1 


iF (a; ¢; 2) = e)Fi(e — a3; — 5). (4B-2) 
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When R (zs) > 0 we have the asymptotic expansions 


Tice (1 — a)(c — a) 
F(a; ¢; 3) ~ 1 
u ; I'(a)2°-@ | T 1!2 
(1— a)(2— a)(e —a)(e —a +1), | 


15? 


a. 


F(a; ¢; —2) ~ 

slit I'(c — a)z? 1!z 
aatidat+t+e-o(Qt+a-—c) . 

+ + 


2!2? 


T'(c) | a(i+a-—c) 


Many of the hypergeometric functions encountered may be expressed in 
terms of Bessel functions of the first kind for imaginary argument. The 


connection may be made by means of the relation®! 


’ 1 z 
(Gs oe 5 oy ~ 1:5) = 2 T(y + 1)s ‘e 1.) (4B-4) 


together with the recurrence relations 


Foy Fax F x | oon I 
i. 2 (a — c) c—2a-2 
2. ac (c —“a)z — cla + 2) 
a: 0 l—c c-—a-—1 
4, =—=¢ at Cc 
5; a-c c-—- 1 l1—a-s3 
6. (c — a)s c(c — 1) cl — ¢ — 3) 


For example, the first recurrence relation is obtained from line 1 as follows 
aF(a + 1;¢;2) + (a — c)F(a — 15 ¢;5) 
+ (¢c — 2a —2)Fla;c;s}= 0 (4B-5) 


These six relations between the contiguous ,/; functions are analogous to 
the 15 relations, given by Gauss, between the contiguous 2/*°; hypergeometric 
functions and may be cerived from these by using 
iF \(a; c; 3) = Limit oF \l a, b; c; (4B-6) 
ho b 
A recurrence relation involving two ,/*;’s of the type (4B-4) may be ob- 
A recurrence relation involving two ,/ f the t tKh-4 y be ol 


tained by replacing a by a + 1 in the relation given by row four of the table 


51G. N. Watson, “Theory of Bessel Functions” (Cambridge, 1922), p. 191 
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and then eliminating ,/'\(a + 1; c;s) from this relation and the one obtaine 


from row 3 of the table. 


F(a; ¢; 3) = iFi(azc — 


There results 


ee ) Pla + tye + 155) 


c1l—ec 


Setting v equal to zero and one in (4B-4) and a equal to 


(4B-7) gives 


ont 
iF, ( 
ify ( 


bho ww 


a 
(55255) 


Starting with these relations the relations in the table enable us to find 
an expression for 1/*\(7 + 3; m;s) where x and m are integers. 
of these are given in Bennett’s paper. 


THE PoWEeR SPECTRUM CORRESPONDING TO y; 


I 
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:s) = o( 
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(Juite often we encounter the integral 
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G,(f) -[ [y(7)]" cos 2xfr dr 


where ¥(7) is the correlation function corresponding to the power spectrum 


w(f). From the fundamental relation between w(f) and ¥(r) given by 
(2.1-5), 
~ w + 
Gif) = be (4C-2 


The expression for the spectrum of the product of two functions enables us 
to write G,(f) in terms of w(f). 


expression: Let F,(f) be 


Hl 


/ F, ine" df, r 


[o etre a 


1, 


? 


(4B-7 


>, € equal to 2 j; 


(4B-8 


A number 
In particular, using (4B-2), 


: 1 2/2 z z 
ni(—); 1; -:) =e Ta +96 () a sl, (;)]. 


(4B-9 


(4C-1 


We shall use the following form of this 
the spectrum of the function ¢,(7) so that 


| Then 
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) Then 
i ao 
Be? / gilt) ge(r)e “ON dr = | Fy(x F.(f — x) dx (4C-3 
Hh] =_ h 20 
: ‘ie, the spectrum of the product gi(r)¢e(7) is the integral on the right. 
+ E e ° « . > ° 
If ¢i(r) and ¢ge(r) are real even functions of 7, (4C-3) may be written as 
@ 1 +x 
[ gilt) ¢e(r) cos 2xfr dr = i] Fi(x)Fe(f — x) dx (4C-4 
~0 2 — 30 
In order to obtain Ga(f) we set gi(z) and ¢e(7) equal to ¥(r). We may 
B-8 then use (4C-4) since ¥(7) is an even real function of 7. When ¢,{7) is an 
even real function of 7 we see, from the Fourier integral for F,(f), that F(f 
must be an even real function of f. We therefore set 
2F(f) = wf), pe 92 
find 
nher | and define w( f) for negative f by 
w(— f) = wif) (4C-5 
B-9 — Equation (4C-+4) then gives 
ee 
Gf) = [ w(x)w(f — x) dx 
L $1. g 
1 f 
= — w(x)w(f — x) dx (4C-6 
8 Jo 
1 a 
= ao i | w(x)w(f + x) dx 
» where in the second equation only positive values of the argument of wf 
rum 7 r 
F appear. 
Vv ’ she 
- — In order to get G3(f) we set ¢i(r) equal to (7), 2Fi(f) equal to w(f), and 
) ¢o(7) equal to ¥(7). Then 
—? ; wo 
i FAf) = 2 | go(r) cos 2xfr dr 
0 
5 us s 
his & = 2G,(/) 


and from (4C-4) we obtain 


; vx 

q G3(f) = 7 | w(x)Gelf — x) dx 
Rigid | (4C-7) 

1 +2 


= w(x) dx / wiy)w(f — v) dy 
16 — 20 — 30 


pone achat ASR Sa. 
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Equation (4C-7) suggests that we may write the expression for Ga(f) as 


1 +20 
G.({f) = 5 : w(x)G,(f — x) dx (4C 


This is seen to be true from (4C-2) and (4C-6). In fact it appears that 
J gy , : 
Gi(f) = 5 wif — x)G,-1(«) dx (4C-9 
-~ ws 


might be used for a step by step computation of G,(/). 

We now consider G,(f) for the case of relatively narrow band pass filters, 
As examples we take filters whose characteristics give the following w(/)’s 
and ¥(r)’s 


TABLE 1 








Filter w(f) for f > 0 V(r) 








/ 
yo f 2so”q2 ond 2 
eS fo) */2e Yoe 2(ror 


ov an 


cos 2nfor 
vo 1 y ent airl os 2 f 

ra + (f— fio? eee 

wf) = wo = ¥o/B for 


pk 3 ee 
a te oar task 





w(f) = 0 elsewhere 





We shall refer to these filters as Filter a, Filter b, and Filter c, respectively. 
All have fo as the mid-frequency of the pass band. The constants have 
been chosen so that they all pass the same average power when a wide band 


voltage is applied: 
yo = | w(f) df = mean square value of I(t) or V(é) 
0 


and it is assumed that fo > o, fo > a, fo > 8 so that the pass bands are 
relatively narrow. 

Expressions for G,(f) corresponding to several values of ” are given in 
Table 2. When » = 1, G,(/) is simply w(f)/4. G2(f) is obtained by set- 
ting 2 = 2 in the definition (4C-1) for G,(f), squaring the ¥(7)’s of Table 1, 
and using 


cos” 2rfor = $ + 3 cos 4xfor 


Filter b 
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The expression for Go(f) given in Table 2 corresponding to Filter 
exact. The expressions for Filters a and b give good approximations aro 


f = Vand f = 2fo where G2(f) is large. However, they are not exact becaus: 
terms involving f + 2fo have been omitted. It is seen that all three G,’s 
behave in the same manner. Each has a peak symmetrical about 2/9 whos: 
width is twice that of the original w(f), is almost zero between 0 and 2/, 
and rises to a peak at 0 whose height is twice that at 2/o. 

G3(f) is obtained by cubing the ¥(7) given in Table 1 and using 

cos® 2rfor = } cos 2rfor + } cos Orfor. 

From the way in which the cosine terms combine with cos 2afr in (4C-1) we 
see that G;(/), for our relatively narrow band pass filters, has peaks at 
and 3fo, the first peak being three times as high as the second. The ey- 
pressions given for G3(f) and G4(f) are approximate in the same sense as are 
those for Go(f). It will be observed that the coefficients within the brackets, 
for Filters a and b, are the binomial coefficients for the value of 2 concerned 
Thus for x = 2, they are 2 and 1, form = 3 they are 3.and 1, and for n = 4 
they are 6, 4, and 1. 

The higher G,(f)’s for Filters a and b may be computed in the same way. 
The integrals to be used are 


—f2/2ne2 
e 


a 
[ e'""" cos 2afr dr = 
0 2o0vV lnur 


ao 
re 1 na 
[ e-'"* cos Iafrdr = se meee ~ 
0 Lee =} 


In many of our examples we are interested only in the values G,(f) for 
f near zero, i.e., only in that peak which is at zero. It is seen that G,/ 
has such a peak only when » is even, this peak arising from the constant 
term in the expansion 


? 7 me 
cos*x = pa cos 2kv + 2k cos 2(k — 1)x + os 1) cos 2(k — 2)x 
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Abstracts of Technical Articles by Bell System Authors 


Historical Background of Electron Optics.’ C. J. Cavsicn. The discov- 
ery of electron optics resulted from studies of the action, upon electrons or 


other charged particles, of electric and magnetic tields employed for the 
purpose of obtaining sharply defined beams. The original Braun tube 
1896) employed gas-focusing, as did the low-voltage cathode-ray oscil- 
loscope developed by Johnson in 1920. It was early discovered that an 


} axial magnetic field could be used to concentrate the electrons into a beam, 


and this method came into wide use in the field of high-voltage cathode-ray 
oscillography. In 1927 Busch published a theoretical study of the action 
of an axially-symmetric magnetic field upon paraxial electrons, showing 
that the equation of the trajectories of the electrons was similar to that of 
the paths of light rays through an axially symmetric optical system. He 
concluded that such magnetic fields constituted lenses for electrons and pre- 
sented experimental confirmation. In 1931 Knoll and Ruska presented a 
large amount of additional experimental material and used the words “‘elec- 
tron optics” to describe the analogy. In 1932 Bruche and Johannson pub- 
lished the first electron micrographs. 

The Davisson and Germer electron diffraction experiments (1927) em- 
ployed electron beams formed by electron guns consisting of a thermionic 
cathode emitting electrons which were accelerated by potentials applied to a 
series of plates containing aligned apertures. The resultant beam was 
quite divergent. Davisson and Calbick made a theoretical and experimental 
study of the forms of such beams. They concluded that the distorted elec- 
tric field in the vicinity of an aperture in a charged plate constituted a lens 
for charged particles (1931). The optical analogy was either a cylindrical 
or a spherical lens, according as the aperture was a slit or a circular hole. 
The theory was confirmed by photographing the forms of electron beams, 
and by construction of an electrostatic electron microscope whose experi- 
mental magnification agreed with the theoretical. 

Coaxial Cables and Associated Facilities? J. J. Pittiop. (Summary of 
Talk before St. Louis Electrical Board of Trade, October 17, 1944.) Coaxial 
cables provide means of transmitting frequency bands several million cycles 
in width over a metal tube a little larger than a lead pencil, with a copper 
wire extending along its axis. Several of these tubes can be placed in a 
lead sheath. 

The frequency band transmitted over coaxial cables may be split up so as 
to provide several hundred telephone circuits or, without such division, 

1 Jour. Applied Physics, October 1944, 

2FM and Television, November 1944. 
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coaxial cables will provide for broad-band transmission service such 
required for television. 
A cable is now being installed between Terre Haute and St. Louis \ 


7 


contains six coaxial tubes to provide telephone circuits, and which may. 
in the future, find use in connection with the provision of intercity televisio; 


networks. WIL 
The structure of the tubes used with coaxial cables consists of a central B Univer 
copper conductor within a copper tube about } in. in diameter, made from § Assist 
flat copper strip which is formed around the insulating discs. Around each Fof Te 
copper tube are two steel tapes which supplement the shielding of the copper § carrier 
tube in preventing interference between tubes in close proximity. The cen- F time « 
tral conductor is separated from the outer conductor by slotted insulating 
disks which are forced onto the wire. The cables are formed with an appro Ric 
priate number of these tubes along with some small gauge pairs used for FU. S. 
control and operating purposes. Lieute 
In the case of underground cables buried directly in the earth, jute or plas- — Ameri 
tic protective coverings are used to assist in reducing sheath corrosion. — Labor 


In some parts of the country it is essential to add a metal covering outside — lems 1 
the lead sheath and the plastic or jute to protect the cables against the f in the 
operations of ground squirrels or pocket gophers. In certain areas thes 


animals have been found to carry away long sections of the jute covering S. 
and will chew holes in the lead sheath unless other metal protection is pro- § Calife 
vided. Copper is sometimes used for this metal covering to assist in light- J} Labor 
ning protection. inves! 


Repeaters in the coaxial system are now located at intervals of about five 
miles. Power for repeaters in the auxiliary stations is supplied from the 





adjacent main stations located at something over 50 miles at 60 cycles over 
the coaxial conductors themselves. 

Coaxial cables are in regular operation between New York and Philadel- 
phia and between Minneapolis and Stevens Point, Wisconsin, a total dis- 
tance of nearly 300 miles. A network of such cables totaling about 7,000 
route miles and including a second transcontinental cable route is being 
planned over additional routes. The requirements of the armed forces, 
general business conditions, the volume and distribution of long distance 
telephone messages, the availability of the necessary manufactured cable 
and equipment, and other factors may modify the extent of this construc- 
tion, the time of starting, and the routes which will be undertaken. 

Western Electric Recording System—U. S. Naval Photographic Science 
Laboratory KR. O, Strock AND E, A. Dickinson. This paper describes 
the complete 35-mm film and 333 or 78 rpm. disk recording and re-recording 
equipment installed for the U. S. Navy at the Photographic Science Labora- 
tory, Anacostia, D. C. Modern design, excellent performance, and ease of 
operation are features of the installation. 


5 Jour. Soc. Motion Picture Engineers, December 1944. 
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